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ports  a  model  of  diffraction  and  interference  near  the  critical  scattering  angle 
provided  the  product  of  the  wave  number  k  and  the  radius  a  lies  in  the  range 
25  <ka<100n.  When  the  scattering  angle  (}><  83° ,  there  is  a  coarse  structure  to 
the  scattering  with  a  quasi-period  ri  (A/a)”^  radians  where  A  is  the  wavelength. 
Tfiere  is  also  a  fine  structure  to  the  scattering  with  an  angular  spacing  propor¬ 
tional  to  (A/a).  Photographs  of  the  scattering  show  some  of  these  predicted 
features.  A  reduction  in  the  scattering  of  parallel  polarized  light  near  the 
Brewster  scattering  angle,  which  is  near  106°  for  an  air  bubble  in  ater,  is 
predicted.  Application  of  these  structures  to  bubble  sizing  and  detection  are 
summarized  and  the  theoretical  extinction  coefficient  in  water  is  plotted. 

Mie  coniputations  for  bubbles  in  water  also  reveal  backward  and  forward 
glory  effects.  These  ere  partially  manifested  as  cross  polarized  scattering. 
Observed  from  bubbles  in  a  viscous  silicone  oil  to  the  near  backward  direction 
is  found  to  have  a  strong  cross-polarized  component. 
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PREFACE 


This  technical  report  consists  of  publications,  manuscripts,  and 

supplemental  information  related  to  the  scattering  of  light  by  bubbles 

in  liquids  and  in  glass.  The  emphasis  of  the  report  is  on  the  scattering 

by  air  bubbles  in  water  into  the  angular  regions  where  diffraction  makes 

essential  corrections  to  elementary  ray  optics.  For  bubbles  these  regions 

are  backward  scattering,  forward  scattering,  and  critical  angle  scattering. 

Most  of  the  previous  literature  on  light  scattering  is  concerned  with 

drop-like  objects  where  the  refractive  index  of  the  '^catterer  exceeds 

that  of  the  surroundings.  The  present  results  should  be  useful  for  the 

1  2 

optical  characterization  of  microbubbles  at  sea  and  in  water  tunnels  , 
or  for  the  prediction  of  the  optical  properties  of  bubbly  media.  They 
should  also  be  useful  for  laser-Doppler  anemometry  with  microbubbles  and 
the  detection  of  microbubbles  in  glass. 

Certain  phenomena  described  here  should  be  present  for  other  spherical 
or  nearly  spherical  scatterers  where  the  refractive  index  is  less  than  that 
of  the  surroundings,  for  example:  air  bubbles  in  ice,  water  inclusions  in 
fused  quartz  (such  as  those  present  in  varieties  of  opal  glass),  ice  spheres 
in  water,  and  superheated  drops  in  a  host  liquid.  The  coarse  structure  in 
the  critical  scattering  region  can  be  present  even  if  the  scatterer  is  not 
spherical . 

The  first  paper  summarizes  experiments,  models,  and  results  of  Mie 
computations  for  both  backward  and  critical-angle  scattering.  Theoretical 
results  presented  indicate  tFat  axial  focussing  enhances  the  scattering  by 
an  air  bubble  in  water  into  the  near  forward  and  near  backward  regions.  This 
enhancement  for  the  forward  region,  (forward  optical  glory),  is  more  signifi- 
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cant  for  depolarized  than  for  polarized  scattering  since  the  latter  is 
dominated  by  ordinary  forward  diffraction.  At  present  the  experimental 
observations  of  axial  focussing  from  single  bubbles  are  limited  to  the 
backward  glory  of  air  bubbles  "'n  a  viscous  polymer  liquid.  (These  are 
described  in  detail  in  the  seventh  paper).  The  first  paper  also  discusses 
the  scattering  efficiency  of  bubbles  and  the  relevance  of  far-field  compu¬ 
tations  to  near-field  observations  of  bubbles. 

The  second  paper  Describes  the  first  observations  of,  and  model  for, 
the  scattering  in  the  critical  region  which  is  near  83  degrees  for  an  air 
bubble  in  water.  This  paper  is  included  here  for  completeness  though  it 
was  not  supported  by  the  present  ONR  contract.  It  was  supported  in  part 
by  ONR  Contract  N00014-76-C-0527  (R.  E.  Apfel  principal  investigator). 

The  third  paper  describes  an  improved  model  (which  includes  both 
diffraction  and  interference)  for  the  coarse  structure  near  the  critical 
scattering  region. 

The  fourth  paper  is  the  first  to  be  published  which  correctly  gives 
the  results  of  Mie  theory  for  angular  structure  in  the  scattering  by  bubbles. 

This  paper  substantiates  the  model  of  the  coarse  structure  developed  in  the 
third  paper.  The  paper  was  excerpted  from  the  M;  ster's  Degree  Thesis  of 
D.  L.  Kingsbury.  A  supplement  attached  here  gives  examples  which  were  not 
published  for  reasons  of  brevity. 

The  fifth  paper  gives  Mie  theory  and  model  results  for  air  bubbles  in 
fused  silica  glass.  The  paper  is  the  first  to  discuss  the  Brewster  scattering 
angle  for  bubbles  and  the  first  to  give  valid  Mie  results  for  bubbles  over 
the  entire  (0  to  180  degrees)  range  of  scattering  angles. 

The  sixth  paper  describes  computer  codes  used  for  some  of  the  Mie  and 
model  computations.  The  codes  shown  here  have  been  adapted  for  use  on  an  HP  1000 
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minicomputer  system.  This  paper  is  excerpted  from  the  Master's  Degree 

Thesis  of  D.  L.  Kingsbury.  For  a  more  complete  description  of  the  Mie 

3 

scattering  algorithm,  consult  the  paper  by  Wiscombe  which  is  the  basis 
of  these  codes. 

The  seventh  paper  describes  the  first  observations  of  backscatteri ng 
from  air  bubbles  in  liquids.  It  also  describes  a  physical -optics  approx¬ 
imation  for  the  backward  axial  focussing  of  scattered  light.  This  model 
shows  that  the  intensity  of  individual  axially  focussed  rays  are  propor- 

3 

tional  to  a  ,  while  the  simply  reflected  scattering  is  proportional  to 
2 

a  where  a  is  the  bubble  radius.  (The  total  backscattering  is  not 

3 

simply  proportional  to  a  due  to  the  interference  of  various  paths; 

see  paper  1).  Polarization  and  quasi-periodic  properties  of  the  glory  are 

discussed. 

The  approximations  and  physical  models  described  hero  for  bubbles 
should  also  be  useful  for  certain  cases  in  the  scattering  of  sound  from 

4 

fluid  or  elastic  spheres.  This  has  been  verified  by  Marston  and  Kingsbury 

for  scattering  in  tlie  critical  region  from  fluid  spiieres.  Backward  axial 

focussing  or  "acoustic  glory"  has  also  been  proposed  for  elastic  spheres 
5 

in  water. 

Care  should  be  taken  when  applying  the  results  of  this  report  to 
light  scattering  by  mircobubbles  in  seawater.  Our  computations  assume 
that  the  bubbles  are  spherical  and  that  the  refractive  index  of  the 
scatterer  is  homogeneous.  There  is  a  natural  tendanry  for  the  gas  in 
microbubbles  to  disolve  i nto  the  surroundi ng  liquid.^  There  is  some 
recent  evidence  that  microbubbles  in  seawater  can  be  stabilized  by 
substances  sorbed  onto  their  surface.^  The  presences  of  a  sorbed  film 
or  small  deviations  from  sphericity  may  alter  the  details  of  the 
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scattering  patterns  describee  here,  but  it  should  not  obliterate  the 

existence  of  axial  focusing  and  the  coarse  and  fine  structures. 

Philip  L.  Marston 
Principal  Investigator 
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Paper  No.  1 

Light  scattering  by  bubbles  in  liquids:  Mie  theory,  physical -optics 
approximations,  and  experiments  (P.  L.  Marston,  D.  S.  Langley, 
and  D.  L.  Kingsbury)  to  be  published  in  the  journal  Applied 
Scientific  Research  in  the  proceedims  of  the  lUTAM  Sumposium 
on  the  Mechanics  and  Physics  of  But)bles  in  Liquids.  This 
paper  was  presented  at  the  Symposium  (June,  1981). 


Abstract 

Angular  structures  in  the  far-field  scattering  from  bubbles  are 
observed  and  modeled.  Mie  theory  supports  a  model  of  diffraction  and 
interference  near  the  critical  scattering  angle.  A  new  expression  for 
the  angular  spacing  of  fine  structure  is  derived.  Photographs  of  scat¬ 
tering  show  some  of  the  predicted  features.  Application  of  these  struc¬ 
tures  to  bubble  sizing  and  detection  are  summarized  and  the  theoretical 
extinction  coefficient  in  water  is  plotted. 

Mie  computations  for  bubbles  in  water  also  reveal  backward  and  for¬ 
ward  glory  effects.  These  are  partially  manifested  as  cross-polarized 
scattering,  Observed  scattering  from  bubbles  in  ^he  near  backward  di¬ 
rection  is  found  to  have  a  strong  cross-polarized  component. 
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1.  Introduction 

Equations  for  the  scattering  of  plane  electromagnetic  waves  by  a  di¬ 
electric  sphere  were  given  by  Mie  [1]  in  1908  and  the  resulting  features 
of  the  angular  scattering  pattern  of  drops  are  well  known  [2-4].  The  Mie 
solution,  though  exact,  does  not  give  insight  into  cither  the  scattering 
process,  or  changes  in  the  pattern  resulting  froii)  changes  in  shape,  re¬ 
fractive  index,  c.  profile  of  the  incident  wave  r'ronts.  Reviews  of  the 
literature  on  light  scattering  [2-4]  reveal  a  paucity  of  information  about 
the  scattering  pattern  of  bubbles  in  liquids  where  the  refractive  index 

of  the  scatterer  n.  is  less  than  that  of  the  surroundings  n  .  Conse- 

1  0 

quently  we  have  begun  [5-10]  a  systematic  study  of  the  scattering  of  light 
by  gas  bubbles  in  liquids.  Aspects  of  the  study  are:  (1)  the  computation 
of  Mio  scattering;  (2)  the  development  of  simple  pliysical  models  wliich  give 
insight;  and  (3)  observations  of  features  in  the  scattering  which  differ 
significantly  from  both  the  scattering  by  drops  and  the  scattering  pre- 
oicted  by  geometrical  optics  [11].  In  this  paper  wo  summarize  the  main 
features  of  the  scattering  with  an  emphasis  on  the  critical  [5-7]  and  back- 
scatter  [10]  regions.  New  experiments  and  aoplications  will  be  described. 

These  are  the  first  detailed  observations  of  scattering  by  bubbles. 

Mie's  solution  [1-3]  to  the  problem  of  the  scattering  efficiency  and 
pattern  of  a  dielectric  sphere  is  a  Function  of  the  ratio  m  ^  ‘ 

is  usually  expressed  as  a  function  of  x  =  ka  -  Zva/X^  where  <i  is  the  sphere 
radius,  k  and  are  the  wavenumber  and  wavelength  in  the  outer  dielectric. 
Optical  sources  are  typically  characterized  by  the  wavelength  in  a  vacuum  A^; 
we  note  that  k  =  2irn^/A^.  For  an  air  filled  bublile,  n^  =  1.00029 
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and  it  is  usually  a  good  approximation  to  take  in  -  l/n^.  Most  of  this 
paper  will  deal  with  scattering  by  spherical  bubbles  with  plane  incident 
waves.  We  will  gain  some  insight  however,  into  which  features  of  the 
scattering  pattern  should  be  sensitive  to  deviations  from  sphericity. 

2.  Scattering  Efficiency 

The  efficiency  factor  Q  is  the  ratio  of  the  total  scattering 

SCd 

2 

cross  section  to  the  geometric  cross  section  na  .  When  n.  is  real 

valuea,  Q  also  gives  the  extinction  efficienc'/  [2-4,12].  Figure  1 
sea 

compares  for  drops  of  water  in  air  with  m  ■-  4/3  with  that  for  a 

bubble  in  water  with  m  -  3/4.  The  computations  were  performed  using  a 
slightly  modified  version  [9]  of  Wiscombe's  MIEVO  Mie  scattering  algorithm 
[13].  A  table  of  for  bubbles  exists  [14]  which  is  consistent  with 

Fig.  1.  For  a  fixed  A^,  and  a  given  value  of  x  ,  the  di'op's  radius  is 
larger  than  the  bubble's  radius  by  a  factor  of  4/3.  For  the  case  of 
light  from  a  He-No  laser,  =  0.6328  ym;  for  bubbles  in  water,  x  of  10, 

100,  1000,  and  10  000  give  it  of  0.‘^5,  7.55,  75.5,  and  755  ym,  respectively. 
The  salient  feature  of  Fig.  1  is  that  for  drops  Q  exhibits  a 

SCcl 

fine  "ripple"  structure  [3,12]  but  that  our  calculations  of  Q  for 

sea 

bubbles  do  not  reveal  a  ripple  structure.  For  both  drops  and  bubbles,  Q 

SCd 

has  a  broad  undulation  with  a  quasi-period  Ax  ~  ir/lm  -  Ij.  This  approxi¬ 
mation,  which  has  been  derived  from  the  theory  of  "anomalous  diffraction" 
[2,3],  appears  to  be  useful  for  both  bubbles  and  drops.  The  ripple  struc¬ 
ture  present  for  drops  is  due  to  optical  resonances  [12]  which  are  attributed, 
in  part,  to  internal  surface  waves.  The  absence  of  such  structure  for 
bubbles  is  probably  because  m <  1  does  not  favor  the  entrapment  of  internal 
surface  waves. 
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3.  Critical  and  Brewster  Anyie  Scattering 


For  scattering  by  drops,  diffraction  is  important  for  the  description 
of  the  forward,  backward ,  and  rainbow  regions  [2,3].  For  bubbles,  there  is 
no  longer  a  rainbow;  however,  a  new  region  appears  known  as  the  critical 
scattering  region  [5].  Diffraction  is  important  in  this  region  because  of 
an  abrupt  change  in  the  amplitude  of  the  reflected  vave  as  the  local  angle 
of  incidence  6  changes  from  Q  <  0^  for  small  impact  parameters  to  9  >  0^, 
for  large  ones.  Here  0^  =  arcsin(iir^)  which  is  the  critical  angle  for  a 
plane  surface.  Figure  2  illustrates  several  ray  paths  which  lead  to  a 
scattering  angle  4)  (the  deviation  from  the  forward  direction)  of  50°. 


The  number  on  the  left  specifies  the  number  of  internal  chords 
denote  the  angle  of  incidence  and  refraction  of  the  pth  ray. 


p;  0^  and 
The  reflected 


ray  (which  has  p  =  0  )  has  a  scattering  angle  4)  =  tt  -  2e„.  For  an  air 

u 

bubble  in  water  with  m  =  3/4,  the  critical  scattering  angle  is  4;  =  tt  -  20  = 

c  c 

82.82°.  Geometric  optics  [5,  6,  11]  predicts  that  |dl./d4)|  «  as  *  ap- 

J 

proachcs  '4  from  above  4  , .  Here  I.  is  the  normalized  scattered  intensity 
c  c  j 

defined  as  follows:  the  actual  j-polarized  intensity  at  a  distance  R  >>  a 


from  the  bubble's  center  is  the  incident  j-polarized  intensity  multiplied  by 
2 

I.  (a/R)  /4.  For  the  electric  vector  perpendicular  to  the  scattering  plane, 

J 

j  =  1,  for  the  parallel  case,  j  =  2.  This  normalization  is  appropriate  for 


bubbles  since  geometric  optics  predicts  tnat  if  the  intensity  of  the  p  =  0 

reflection  could  be  considered  by  itself  I .  (4)  ^  4)  )  =  1  due  to  total  reflection. 

J  c 

Our  observations  [5],  model  [5,  6],  and  Mie  computations  [7-9]  demonstrate 

that  instead  of  a  divergence  of  ldl./d4)|  ,  tlie  rise  in  I.  is  spread  out  over 

J  J 

the  region  14)-4>^1  =  ^  where  Q.  •=  arcsin  [0.8  (l-m^)"''^(X^/a)'"]  ^  144x"^deg 
for  water.  The  model  makes  use  of  a  physical  optics  approximation  which 
is  to  (a)  use  ray  optics 
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along  with  the  reflection  coefficients  of  a  plane  surface  to  compute  the 
complex  amplitudes  of  virtual  waves  which  simulate  the  scattering,  and  (b) 
use  Fraunhofer's  approximation  to  compute  the  diffraction  of  the  virtual 
wave  to  the  far  field  where  R  »  xa.  Step  (b)  yields  the  improvements 
over  the  geometrical  predictions  [11].  This  procedure  is  analogous  to 
Airy's  model  of  diffraction  near  the  rainbow  [2]  .  Additional  approximations 
used  in  the  model  are:  only  the  two  most  intense  virtual  waves  for  (f.  -  cj)^ 
are  included  (p  =  0  and  1);  and  the  divergence  of  the  derivative  of  the 
reflectivity  is  simulated  by  truncating  the  p  -  0  reflection  when  dg  <  0^  . 

To  compare  the  model  results  with  Mie  theory,*  we  consider  examples  of 
bubbles  in  water  which  complement  those  previously  published  [7].  Figure  3 
shows  both  theories  for  x  =  100  and  j  =  1.  It  shows  that  there  is  a  coarse 
structure  in  the  Mie  scattering  for  (j)  <  which  fins  a  quasi-period  ^  --  14.4''. 

This  structure  is  described  by  the  model  except  in  tlio  near  forward  direction 
(4>  2^). 

The  Mie  result  also  has  a  superposed  fine  structure  where  the  magnitude 

of  the  quasi -period  <  A  fa  rad.  =  360"x"^  .  Near  ,  this  fine  structure 

0  c 

arises  primarily  from  tlie  interference  of  the  p  =  0  wave  with  the  wave  due 
to  the  p  “  2'  ray  in  Fig.  2.  Its  quasi-period  may  be  estimated  froii,  the 
lateral  separation  bg  +  b^,  of  virtual  point  sources  which  would  simulate  the 
scattering  at  a  given  0  .  As  shown  in  fig.  2,  this  separation  is  the  sum  of 
tlie  impact  [larameters  b^  for  these  rays.  Standard  relations  for  tlie  far- 
field  interference  applied  to  these  sources  gives: 

arcsin[A^/(bQ  +  b^,)]  (1) 


*The  conversion  from  the  Mie  amplitudes  S.  [2,;-!  13]  to  the 

(2|S.|/x)^ 


I  .  is 
J 
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where  b^/a  ^  ~  't')/^,  (^  =  rr  +  2(02,  *  ^^2''^ 

law,  nisinpp  "  (U  is  only  an  approximation  because  it  fails 

to  include  other  (such  as  p  -  1)  virtual  vmves  and  it  does  not  completely 
account  for  the  lontj  itud  inal  spacing  of  virtual  sources  and  for  how  those 
sources  vary  with  '}).  With  ()■  ~  4)^  and  m  =  3/4,  b^  +  b^,  ~  a/0.825  and 
=  arcsin  (5.18  x  ').  This  gives  2.97”  at  x  •--■  100  which  is  in  reasonable 
agreement  wi th  Mg .  3 . 

Figure  4  shows  the  Mie  and  phys ica  1 -opti f.s  results  for  x  ^  10  000  and 
j  =  1.  The  fine  structure  quasi -period  is  greatly  reduced;  Fq.  (1)  gives 
Q^((};^)  0.0297°  which  agrees  with  the  Mie  result  of  0.03°.  The  amplitudes 

of  the  coarse  undulations  (Q  1.4°)  decrease  slightly  with  decreasing  T  until 
4)  -  75°.  A  graph  of  for  |i  •:  75°  shows  coarse  undulations  increasing 
in  amplitude  with  decreasing  T-  The  explanation  is  that  tor  x  =  10  COO, 
with  4)  75°  the  coarse  structure  is  primarily  due  to  diffraction  of  the 

p  =  0  wave  [5]  but  for  4.  75"  it  is  primarily  from  the  interference  of  the 

p  =  0  wave  with  p  >  0  waves  [6].  Plots  [7]  of  1^  for  this  x  also  show 
thi s.  transi  tion  and  that  fine  structure  is  significantly  weaker  when  j  2 . 

Figure  3  does  not  show  the  transition  because  the  diffraction  region  (f  -•  d  ,<  Cl) 
overlaps  the  region  where  interference  with  the  p  =  1  wave  is  significant, 

As  Gg  approaches  the  Brewster  angle,  0^^  -  arctan(m),  the  reflectivity 
of  the  p  =  0,  j  =  2  polarized  ray  vanishes  [2].  Consequently,  according  to 
geometric  optics,  there  is  no  contribution  to  fnom  this  ray  at  the  Brewster 

scattering  angle  4.,,  =  n  -  2Up  -  2  arctan(m  ^).  Because  the  reflectivity  varies 
slowly  near  (p^^,  diffraction  is  less  important  than  near  4^.  A  coarse  minin,um 
is  evident  in  the  Mie  near  4^  for  bubbles  in  glass  [8]  and  for  bubbles 
in  water  where  4)^  ~  10G.3”.  When  x  drops  below  5,  its  location  shifts  toward 
90°;  the  scattering  pattern  approaches  t'^it  of  j  dipole  radiator-  predicted  by 
Rayleigh  scattering  theory  [2,  3,  8]. 
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4.  Observations  of  Critical  AtKjle  Scatteriruj 

Previous  observations  [5]  were  limited  to  the  4)  region  where  the 
coarse  undulations  decreased  with  decreasing  -ji  .  We  have  a  new  apparatus 
similar  to  that  in  [5]  except  that  it  permits  observations  with  4)  down  to 
70°.  As  in  [5],  bubbles  were  attached  to  a  vertical  needle  in  distilled 
water.  They  were  illuminated  with  a  plane  wave  from  a  He-fle  laser  with 
>1^  =  632.8  nm.  They  were  pliotogra[)hed  via  a  windc'V)  with  a  camera  focused 
on  infinity  to  place  the  film  in  the  far  field.  Because  of  a  technical  problem 
we  have  not  obtainetl  precise  direct  measurement^  of  e  in  the  (horizontal) 
scattering  plane  along  with  the  photographs,  liowevor,  we  have  demonstrated 
that  )  from  Eq.  (1)  gives  an  a  which  is  consistent  with  the  observed 

T  C 

and  modeled  coarse  structure.  Thus  the  '■  values  (lunted  in  F-g,  5  were 
determined  from  the  coarse  and  fine  structures  and  not  from  direct  observations 
Figure  5  servos  more  to  illustrate  phenoiicna  des(;rii)CHj  in  Sec.  3  than  to  rig¬ 
orously  test  the  models  v<ith  real  bubbles.  The  angle  scale  is  for  a  4) 
change  in  water.  The  horizontal  displacement  is  not  exactly  linear  in  4)  due 
to  refractive  corrections  at  the  water-wiridow-ait  interfaces.  In  Fig.  5(a), 
the  broad  bright  region  on  the  left  is  the  first  ccarsn  maximum  with  4-  <  4)^. 

On  the  right,  the  coarse  undulotions  iricrease  in  amplitude  as  4)  decreases, 
which  is  the  modeled  behavior.  Figure  5(b)  is  for  a  smaller  bubble  (verified 
v/ith  direct  observations).  The  coarse  undulations  are  spread  out  and  the  fine 
structure  is  clearly  visible  with  •  0.07°. 

5.  Backv;ard  and  Forward  Glory 

Backscattering  from  drops  is  known  to  be  enhanced  by  the  axial  focusing 
of  certain  rays  [2,  3],  A  complete  descrij)f'on  of  this  "glory"  for  water  droi-s 
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is  complicated  due  to  the  ner  ssity  of  including  surface  waves  [15].  We 

have  observed  and  modeled  backscatter’ing  fnom  air  bubbles  [10]  and  find  that 

some  aspects  can  be  explained  with  a  physical -optics  approximation  that  does 

not  include  surface  waves,  for  the  following  discussion,  it  is  convenient 

to  define  y  ”  -  (I*;  y  -  0  corresponds  to  exact  backscattering . 

Computed  Mie  intensities  sliow  that  liackscattering  froiu  bubbles  may  be 

significantly  larger  tiian  tliat  expected  from  re-'^lection  of  the  p  '  0  ray. 

For  y  =  0,  the  plane  of  scattering  is  no  longer  defined  and  I.j  =  =  I. 

'  igifc  6  shows  several  broad  peaks  in  !  whicli  arj  significantly  larger  than 

the  normalized  intensity  of  the  p  -•  0  reflection  which,  by  itself,  is  (m  -1)'7 
2 

(m  +  1 )  =  0.02.  Other  axial  (e.g.  p  2)  rays  are  too  weak  to  explain  the 

magni tude  of  I . 

For  polarized  incident  light,  scattf'riu"  near  y  -  0  may  be  described  in 

part  v.'ith  (non-cross)  polarized  and  cross-polarized  normalized  in- 

(2) 

tensities.  For  I  ,  the  plane  of  polarization  is  rotated  by  90"  from  incideiU 

polarization;  symmetry  gives  0)  0  wliile  l^^^(y  0)  =  I.(y  =  0). 

in 

For  y  ^  0,  the  I  de[)(’nd  on  both  y  and  the  angle  f,  which  the  scattering 

plane  make,  with  the  incident  electric  vector.*  Figure  7(a)  illustrates  the 
computed  scattering.  Geometric  models  of  l^'^(y  0)  predict  that  the  con¬ 

tributions  of  off-axis  (p  -  3  and  4)  rays  diverge  as  y  ->■  0  because  of  a  factor 
which  accounts  for  focusing.  This  divergence  is  present  in  previous  ray  models 
[11],  but  was  not  discussed;  it  is  also  present  for  drops  [2]  for  certain  m- 
Because  of  this  divergence,  diffraction  provides  an  essential  correction  to 
ray  optics. 


*Th0  conversions  from  the  compie/  Mie  amnlitudes  S.  to  the  I^'^^(y<<  10°,  4 

^  .  .(2)  _  ,2..  2 


are  I 


(1 


S  -  S  cos2'. i  /x  and  I 


S  sin2fl  /x  where  S' 


S.|  i  S^.  Evidently  the  (ippendeiu.e  of  Eg.  (2)  is  of  ttie  correct  type. 
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The  scattering  due  to  the  p  -  3  ray  in  fig.  /(a)  (insert)  is  a  back- 
facing  toroidal  wave  front  which  appears  to  originate  at  a  virtual  ring¬ 
like  source  known  as  the  focal  circle  in  the  analogous  {?.]  p  -  7.  scattering 
from  drops.  This  source  at  point  F  is  rinqlike  because  the  figure  may  be 
rotated  around  the  optic  axis  through  tlie  center  C.  for  largo  x,  the 
stationary-phase  approximation  of  far-field  diffraction  integrals  gives  the 
following  proportionalities  for  off-axis  (p  >  2)  contributions  to  the 
taken  separately: 

/  1  ^  9 

I(p)  «  x[(c.|  +  c^)  dQ(u)  4  (c.|  -  c^)  d,,(u)  cos2;;j'  ( 


ijpj  sin2.-,]^ 


where  the  c.  are  combined  Fresnel  reflectivities  and  transii)issivities  of 


the  j-polarixed  fields,  u  ^  kb^^^sin  Y  ,  and  bp 


is  the  iiipact  [jarameter 


of  the  exactly  backscotterod  pth  ray.  Because  l  |pj  «  >.  ,  the  unnorinal  i  zed 

intensity  is  proportional  to  ka“’  while  that  for  the  p  =  0  ray  is  proportional 
2 

to  ;/  .  Consequently,  off-oxis  waves  are  dominant  when  'i  is  sufficiently 

large.  It  is  evident  from  Fig.  b  and  7(a)  that  the  modeling  of  ii'ay  re¬ 

quire  summation  of  electric  fields  froai  several  virtual  sources.  (In  physical 
optics,  intensities  do  not  add  but  fields  do.)  Interference  of  various  ring¬ 
like  and  axial  virtual  sources  depends  on  x  and  could  lead  to  the  broad  struc¬ 
tures  in  Fig.  6.  Detailed  computation  of  (y  =  0)  give  magnitudes  suf¬ 

ficient  to  explain  the  flic  scattering.  Due  to  the  Bessel  functions  in  Fq.  (2), 
the  actual  scattering  should  be  peaked  at  or  near  y  '  0. 

We  have  observed  backscattering  from  sin.ijle  air  bubbles  in  a  dimethyl- 
si  loxane-polyrior  liquid  which  were  nearly  immobilized  by  viscosity.  In  these 
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experiments,  m"^  =  1.403,  =  6.3;^. 8  nm,  a  -  0. 3  -  0.8  mm  and-  x  -  4000  - 

11  000,  The  far-field  cross-polari._'ed  intensity  had  a  dependence  on  y  and 
similar  to  that  predicted  by  Eq.  (l^b)  with  p  -  3.  Equation  {2b)  predicts 
that  when  both  siny  y  3nd  u  ;■>  1,  the  minima  in  should  be  spaced 

by  Ay  ~  IT  (kb^)“^  radians  where  b'^  -  0.44/j.  The  prediction  correctly 
described  the  observed  <  -  Z  scattering  which  is  ap()arently  dominated  here  by 
the  p  =  3  virtual  source.  The  focal  circle  was  viewed  by  focusing  the  camera  on 
the  bubble.  TheE.=  1  scattering  vws  not  dominated  by  a  single  class  of  rays. 
Cross-polarized  scattering  from  a  polydispersicn  of  bubbles  in  water  fias  also 
been  seen. 

Figure  7(b)  shows  an  enhancement  of  near- forward  cross-polarized  scattering 
due  to  axial  focusing.  Focal  circles  due  to  the  p  =  ?  and  3  rays  la»'gely 
contribute  to  the  C  ■  2  scattering.  The  forward  p  =  1  scattering  is  don)i- 
nated  by  ordinary  diffraction  [2]  when  x  is  large.  Hie  theory  gives  10^ 

when  (f'  =  0  for  the  x  and  m  of  Fig,  7.  Forward  optical  glory  has  also  been 
displayed  in  Mi"’  scattering  by  water  dro|)s  [16].  Its  description  is  complicated 
by  surface  waves . 

5.  Discussion 

Figures  3-7  demonstrate  that  for  several  angular  regions,  the  intensity 
exceeds  that  of  geometric  scattering  from  a  perfectly  reflecting  sphere  of  the 
same  a  (for  which  K  =  1).  This  information  should  be  useful  in  optical  de¬ 
vices  which  size  or  detect  bubbles  [7,  17,  18j.  Due  to  diffra.tion  near  4>^  , 
it  is  preferable  to  detect  scattering  with  i-  4)^  -  w  than  to  observe  it 

with  ()j  =  90'’  which  is  the  usual  practice  [1.'  1  for  bubbles  in  'vater  For  an 
unpolorized  source,  the  normalized  total  sea i tiM'ing  is  il^  ‘  1^)/?. 
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Though  this  paper  has  emphasized  the  far-fielci  i)'  dependence,  the 
results  are  also  applical)le  to  the  imagitig  of  bubbles.  For  example,  if  it 
desired  to  resolve  the  virtual  sources  of  the  p  =  0  and  1  rays  in  Fig.  2, 
the  aperture  of  the  imaging  system  should  have  an  angular  width  >  a  .  It 

'''j 

is  easier  to  resolve  the  p  -  0  and  2'  virtual  sources  since  the  angular 
width  requirements  are  reduced  to  Resolution  of  identifiable  virtual 

sources  reveals  the  size  of  a  bubble. 

Fine  structure  in  the  far-field  scattering,  S'lch  as  that  shown  in  Fig.  5(b), 
arises  from  the  interference  of  widely  spaced  rays.  C  )nseqiiefi tly  the  positions 
of  the  maxima  will  be  quite  sensitive  to  changes  in  the  bubble's  shape.  That 
part  of  the  coarse  structure  due  to  the  interference  o\  p  -  0  and  1  rays 
will  also  depend  on  the  shape  but  more  weaKly.  As  is  the  case  for  drops  [19], 
tl'iese  shape  dependences  riay  be  useful  for  detecting  meciianical  resonances. 

Scattered  intensities  for  bubbles  in  water  p'rujrted  to  be  from  y.io  tlieory 
have  beer,  used  in  a  study  of  cavitation  nuclei  Comparison  of  our  Hie 

results  and  model  of  tiie  coarse  structure  (which  are  consistent  for  the  a  in 
question)  with  those  in  [17],  show  that  the  latter  err  siqnificantly.  Plots 
in  [17]  predict  coarse  maxima  at  =  100"'  for  r,-  =  5  tun  arnd  >]j  =  120"'  ■for 
a  '  7.5  iim  which  are  not  present  in  our  Mio  computations.  We  cannot  find  any 
physical  justification  for  coai se  maxima  at  these  p  for  x  near  100. 

We  are  grateful  to  W.  J,  Wiscombe  for  providing  the  initial  comouter 
program  from  which  the  program  used  hero  was  derived.  This  vesearcli  was 
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Figure  Captions 


Fig.  1  Scattering  efficiency  from  Mie  theory  as  a  function  of  the  size 
parameter  x  .  The  solid  curve  is  for  drops  and  the  dashed  curve  is  for 
bubbles . 

Fig.  2  Ray  paths  for  a  bubole  with  rn  =  3/4  .  The  number  of  internal 
chords  p  is  given  a  prime  if  the  ray  enters  the  bubble  below  the  center¬ 
line. 

Fig.  3  Logarithm  (base  10)  of  the  normalized  scattered  intensity  pre¬ 
dicted  by  flie  theory  (solid  curve)  for  x  =  100.  The  dashed  curve  is  the 
physical -optics  approximation  [6]  of  the  coarse  structure. 

Fig.  4  Like  Figure  3  but  with  a  linear  scale  and  x  -  10  000. 

Fig.  5  Photograpns  of  far-field  scatterinc)  for  (a)  the  bubble  r'ldiur, 

'I  400  i,m  (x  '  6360,  ■  1.8”)  and  j  =  2;  (b)  330  i.ni  (,\  ■  43/0,  ■  2.2 

and  0  ■  1.  The  critical  scattering  angle  is  near  the  left  edge  of  ihc 
photograph;  decreases  from  left  to  right.  Coarse  structure  is  manifest 
as  broad  vertical  bands.  The  rings  are  artifacts. 

Fig.  6  formalized  Mie  backscatteri ng  froiii  a  spherical  air  bubble  in  water 
(m  -  3/4)  plotted  as  a  function  of  the  size  para'iieter  x;  for  green  light, 
a  200 

Fig.  7  Normalized  Mie  neac-backscatterinq  (a)  and  near-forward  scattering 
(b)  for  m  =  3/4  and  x  =  3040.  For  clarity,  the  computed  was  multiplied 

by  two  in  (a)  before  plotting.  The  inserts  illustrate  some  of  the  significant 
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Critical  anrj!0  scattering  by  a  bubaia:  phycical-cplics 
apprtixirnatlon  and  observations^ 

Philip  L.  Miirston 

Defxirln.rnt  of  l‘h\sr“.,  Woshinulon  Slolo  (  nii  rrsity,  ['iillrnon,  Woxhington 
DcparlrJirrU  of  f’fipir.fi”  inp  A/i/ilcrl  Si  n'ni  f.  Yolo  ( 'ui v.  .Vric  //tji  i-n,  C'onni'cdrul  OfsiW 

(Ki'ceivca  a  rrtiiimrv  IdV'tl 

The  intensity  of  liijht  soMttereii  fiv  nn  air  liuf'l'le  'O  i*  iter  m  pn  ilirlrwf  by  iIhi  peorietric-optirs 
calrulHtion  'f  Ihivii  iKt'ifo  to  a  .ini'i.b'f  derivative  a.H  tfie  i  ritual  aMitterutf:  un^lu' 

<J),  ifl  npptoiii  hed.  I'-ll'  a  t-i  (,f  diO'rarlion  I  t  ibn  antpilar  r-.'iriiin  near '!>,  are  (lem  ril.ed  here.  Ttie 
Fraunhofer  dilTr.action  for  aeatt' rinir  an.;'.-  a  eslimati  d  uMini;  a  einipiilied  |diyai<'aU)ptics 

Bpproxinnilion  A  rirpinK  and  nec.iv  of  'be  nr  field  intenailv  \  >  i.redirte.i  that  i.-i  fai’iiallv  iintila,'  to 
the  nenr-fleld  diflracUon  at  a  ■'iraitht  eof  c  Ob'-  .  v  iiion  ef  n  ilbmeier  radiim  ti.ibble'  in  water  lanth 
collimated  mor.orbratnn''C'  liluniinatiiei  rntifn no  rti"  ‘•xisieta-e  at  ibia  rininn,-  ,  !i  ha"  a  qua  a  pet '.id 
■r  25  mrad.  The  difIVartinn  rali  iilati'ai  ;n\ej,  an  apprexi.-nale  di  arription  of  tea  relative  il>  of  the 
ohaervod  maxima  and  iiiinima  I’l  no'".-'  'A  iih  a  lower  romra.m  ami  apaci'ii;  ■  I)  1  mrnd  were  also 
observed;  they  appear  to  be  caused  by  tb  ii.o  rfereme  of  tavs  witb  distiiu’t  patbst.  Impliealions  for 
the  critical  niiftlo  scattering  of  w  hile  liKht  are  iliwu».seil. 


INTRODUCTION 

T'he  piir|)nse  nl  this  p;iper  is  to  prevail  ine.i  .  iM  im  tils 
anti  till  jipproviin.ili'  I  heorel  a.il  tli'i-i  ri| 'i  "ii  "I  till  I  r,.i  I  H'II 
pheiinnieli.i  nliservcil  w  lieii  bi'.bl  I  .  si  .'ili f  d  near  llm  erilii  al 
ant’.le  irem  a  laiMile  wit  h  an  iii'i  rmr  I  'drai  live  index  li  than 
that  of  the  taller  tin  di.i.  'Idle  i  ipi  i'  al  wax  id.  npih  is  a'snmed 
111  he  very  iinieh  less  lhaii  I  hi  i  adiiis  nl  the  nearlv  pin  i  a  a  I 
htlhhle.  Appi  oxiinal  intis  t'.ix cn  here  I  Ota  ernie.;;  I  In  'lehtix mr 


111  the  ret  bad  inn  i  ml  I'ii  lent  a  ml  the  exit  ill  nf  i  e-.ei  .■  in;;  wave 
Irnnt  \  II  Id  an  I  xpn  'Sinii  tnr  the  seal  lei  ili'd  in  I  .ic  l.ir  1  lelii  Ui.il 
m  sinul.ir  in  I  nr  in  In  I  he  I  i  esiiel  lunar  In  Id  I  d  il  I  'M'  I  inn  h\  .1 
I  rue;  hi  iai,-,e  Oh  .et  \  ai  mos  ri  purl  id  here  nl  I  hi‘  ■■  an  1  'ini' 
In.  .1  nil  hid  I  in  x\  al  1  r  Mippni  1  sniiie  aipei  Is  nl  I  'ns  ^iiiiiihl  ail 
Ihentv.  The  IriMlnieii'  mav  he  ri  .iililv  exlentliai  In  ineliide 
se.illerini;  I  rnm  nl  her  i  iirv  lai  'ill  l.n  i"  'i  iidi  a'  I  h  .1  n|  a  ex  din 
ill  r,  pi  nv.  dial  I  ha  I  the  ni  11  in  pie  latlii  n|  1  all  xai'  iiif  are  nineh 
Uieili  I  'h.in  the  xvax  eleni;l  !i  nl  lipht. 

I  1 '.r.  'M  )pl  a  il  .-siu  let .  el  nnierii  a  1 . 0.''i 


l'20a  J.  Opt.  Soe.  Am  .  Viil  li'.l,  No.  It.  bi.pieadier  IliVIt  ■  .  .IP  1 1 '.M  o'.'l  ■  iv,  n,,,  .q 


MG.  )  loMI  ii  tli'rhon  in 
plan**  ol  sc  ani'i  hkj  ol  a  »  ly  with  an 
inipaii  papanintiM  h  laKjnr  ih.iii 
thn  i:nUi.ai  inp.ai  t  paianw??#*!  f* 

( )nly  i)n»  qna  ltant  (*f  tiu*  ufnss 
snciton  nl  M«*  l)iit»t»l«?  i;.  shewn 


Oavi.s'  rali'iilalcd  the  aniMilar  disl  i  iliut imi  nl  llif  mlfusiiv 
of  ii^ht  sciitli'rcd  liv  an  air  Indililr  in  w  it  ri  l  roni  ,i  inl  Inna  led 
ilU'lvCent  lieatii.  I  lie  Cidi  illal  inn  w  l^  ti.i  ei  I  on  la  nniel  ni  n|il  tea 
and  the  ('’resnel  rel  lia  l  am  i m  1 1  n  a  nl m i|  In  In  .it  I'l.me  miei 
faet";;  ealiailat  nain  nl  I  In'.  I  \  | le  nun  i  a  1 1 1  >  an  .i|>|irn .  nn.i'  mn 
that  is  uselul  nldv  when  the  railin',  nl  the  laihlih-  .i  i  m  iv 
tmii  h  ►treater  lli.ni  the  w  n  eh  n;,l  h  in  I  he  enter  nieil,.i  Let 
( he  seal  t  ern  1 a njt  le  i/i  nl  ,|  i  a dellnt  e  I  I  le  de  V  i.n  n  a  i ii;; le  nl 
f  he  seat  I  ered  r.it  I  mm  I  he  d  i  re.  l  mn  nl  t  ii.  .  nee  ii  nl  ii.ii  allel 
hea  ni  nl  I  avs  I M  e  I  i.;  I  I.  .'V  si,;  i  n  1 1.  ,i  nl  h  ,i  I  in  i  •  4  I  i.i\  ■  .  i  .d 
cnhitinli  is  |h  II  ihere  slimild  hi  .m  nniipl  dei  ma  e  III  the 
seal  ter  ed  iiiletisilv  when  the  m  at  i  e'  n";  an,;le  .  esi  ■  i  n  ,  .s, 
where 


ip,  '  rr  ~  '211  .  Ill 

sinh,  ri  (■.!) 

and  /I  =  >  I,  I'l,,  .nid  'i,  heinr,,  i'  .pi'ei u •  h. .  the  re 

frat'live  iiiihi'e.s  nl  the  nnler  and  the  niin  i  inednt. 

Fhtnre  1  ilhisiraies  the  phvsieal  nnein  nl  tins  predieted 
ileerease,  Cnn'.’ider  Imw  '  he  ,ni;;le  nl  im  idem  e  .it  'he  hiihl'h 
siirfare  9  varies  a-  the  iini'.e  t  paraineli  i  nl  the  r.iv  /■  ni<  re.ise  ., 
Uavs  with  an  inipael  pai.iiiiert'i  in  e\ress  nl  .i  erilu  ,il  value  />, 
have  1/  ?  I/,  vs  here  Ii, ,  ttiven  hv  I'ii]  r,').  IS  1  he  ent  teal  anitle  Inr 
a  p/rifie  inlerl  lee,  .'Vei'.irdmi!  in  the  c.eninel  rie  uptii-v  ap 
prnxillialinii.  ravs  with  /■  ■  h,  n  e.  'I  will  he  Iniallv 

re'leeted,  UaVs  with  li  h,  will  In  p.nliallv  lelleelid  .le 
ia..din)t  In  the  inaanilllde  nl  the  I  r'n..nel  relleetinii  |■nellielents 
aiul  I  his  |',.irl  lal  lellei  I  inn  le.id  .  In  lie  .ilirnpi  de''re:>  -e  in  i  la' 
seallered  inteii  atv  Inr  ■:>  '  '!•  i  he,  i e  ileiil.ati  d  .il  n  rnit; 
inlinsitv  alsii  me  hides  I  he  mnll  iple  t  el  h  .  iinli  and  I  r  iiisini  • 

Sinn  nl  the  rav  s  w  ilh  impaet  p.n  rUiieli  I .,  /  ■  ■  /i,  The  r.ivs  that 

penel  rale  I  he  sphere  .ind  aie  p.irl  i.illv  rellei  leii  i  niil  ■  daite  In 
seal  It  14;  Inr  lie  enl  1  re  ra njte  n|  ,,,,  Imw'i  v a  r,  I  lie\  dn  iint  nh 
seiire  I  lie  d  i>l  nil  I  deerease  ill  tile  se.il  1 1  I  "d  inteiv.n  '.  Inr  '/■ 

•A,  This  IS  ev  ideal  'rnin  1  nil  in;;  piedi' '  i.>n  I'e'i'  !S|,))<|.| 
II  Inr  the  i nil' ns 4  \  1  il  ,111  ,n  1  hnl 'Me  1  n  \  r  ■  1  1  h  n  h  e.  i,'*. 
8'.!.^*’ ,  it  i  .  .il.'i j  e  \  idi  Ml  in  1 1  le  I  nl  r'  I 'i ■  1  in:;  p: '  ■  le  t  inie  Inr 

hnhiiles  in  nlla  r  liipnds  ll'i",,  1  n|  l<ei  111  Tiie  iinensilv 

enni  t  ihul  inn  nl'  I  lie  mnll  iplv  1 1 4  leei  rd  .mil  rel  1  .n  lid  rav-  will 
he  ninitled  in  Itie  In  alment  presenli  d  here  sniee  wa  are  i  .ni 
siderin^;  inilv  the  elleels  nl  dillrm  lam  nn  the  seallered  in 
tensity  (nr  small  i',v  -  :•  ! 

'I'lu*  ealeul.iliniis  nl  I  hi  vis  predn  I  I  h.il  i'l  di-annl  in 

Uniis  at  4i,  where /,  I '/4  IS  I  he  Ini  .1  .iln  ied  inii  iisitv.  1  hi 
Iriu'l inn  nl  lij;ht  will  imi  pern,  11  ill,  'h;>  1  n  he  d '  ,1  iml miiniis. 
The  met  hi  >d  pri. I'll  I  I'd  here  is  .1  .iia(4i  nn  1  le  d  Ini  1-4  nn.il  in;; 
the  id/eets  nl  ddlf.ielinii  n  llie  reitinn '.eh,  u  n  ■  '■  |  1 ,  .m.ill 

I'  irst,  (he  apprnximale  shape  nl  I  he  It  I  n  il  w  ive  1  mill  e.  1  .il 
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cidllted  ttenim  I  re  illv  hn  tlin se  r.iv.  with  h  h,  liays  wil  h 
b  <  l>,  will  lie  mini  led  11  e  Si  e  II),  The  virl  iial  wave  I  ri  ml 
is  then  inndil  n  il  h\  .1  ph.e  e  1  m  reel  mn  line  (11 1  he  phase  slid  1 
nl  the  'mimi>dn).n  i  n  ni  l  .ni.pliinde  retlei'tinil  enellieienis 
Inr  b  •  b,  rile  iin " I il  n  d  vs .1  v  r  I  rmil  ile'-en lies  the  apparenl 
niirei  viev.eil  lev  m  'ii  i  'vi'i  lal  Irni'.i  I  he  hnhhie,  I  he 
Kr.nmhnlei  iliUi.i,  i-.m  nl  lln,  wave  eives  the  desired  ap 
prnx  jiii.n  inn  .  im  iln  ..  iiPim:;  amplilnde  and  mlensiU 

The  vv.e  e  Ir  an  1  .  .  .  imed  I  n  extend  nver  a  hall  plane  and 

is  hot  hmiled  li\  I  III  Imil  I  si  re  III  t  he  linhhie.  Tim  I  real  men  I 
resenihles  Airv  ■  n  llr.iilmn  lairreetinn  to  liesearli-'  a''" 
melrii  iipln-.  Ue.ilnieiil  ni  1  he  ramhiiw  e.mslie  '  '  I'.ven 
(hi until  .'Xirv ' .  .ippinxmi.il h III  1  niisiders  the  virlnal  wave  1  runt 
In  he  n!  ml  ..nil  e\i  1  nl  1  .mil  .1  1 1 1 1  ip  imist  he  linile),  il  is  iisel  nl 
lot  the  dniinn  ml  pnlan  .-al  inn  when  the  radius  ni  I  he  spiiere 
rx<  eeil  the  w  .IV  i  h  cr',1  li  iiv  1  1.1.  Inr  nl  apprnxililalelv  Mill,'  ' 
lie  -ml  a  ipaled  I '  n.ininn  s|  n.  i;'ein  erilerinn  in.iv  deinieale 
till- r. lime  nl  1.  ,'i',l',  ,  'I  lln  ,ippin\imahiins  preseiiled  here 

Inr  dill  I  .III  mil  'n.  .1  Is  hi  lie  m  1  he  re;;  I  nil  nl  i  ril  ieill  .untie 

M  e.i  III einei  '  4  •  .  , M  e;  m;;  i  4  a  /| 1 1 e  (  I iv  ah  ,m 

In  ihi  lie  in  'V  .il  e.  n  i  ■  ■  -i  i  ed  1 1\  1  i  I  ■  -i  i  I  I'l;'.  I  nl  i(e(  7)  I  lav 

-III  ■  '  ,r  ll  ml  e  I  ’  pi  ei '  I,  I  :  .1  1  nl  ,  e  ' .  1  11  r  I  ipl  i  IS  1 1 1  11  t  a  p  il  I  d  e 

.  re.i  ...  in  /.  :  n  ■  He  ise  .  4  1 1 1 1  e  1 1'.'.hl  a  lid  I  he  linn  I  ed 

anr.nlaiii'  ilninnilihi  ’m  i-m  i  meni  make  I  I’.eni  nnsinlahle 
tiireniin-'''  vi>  n  ■  1  "  I  I  el  e  1 1 1  1 1  a  1  It  V .  T  111' 1 4  ise  r  V  a  1 1  HI 

with  la  sei  . . nn  1 1  pnr  I  ei  I  lie  re  ni  Si  .  .  1  V  mil  ii  al  i  1 1  i.ii 

the  -impii  lilli.i  I .  Sieniv  l  ives  a  l.nr  desiTipIlnil  ni  the 

lelai  ive  pn  .ii n  ;i  ■  la.t  ■  1111,1  ,mil  niminia  nl  I  he  mteii'.n  V 
III  the  le  .Inn  '  .  ill.  ..1  .Ih.  :e  -.i  al  lenint 


I.  CURVATUni  OF  IHE  VIRTUAL  WAVE  FRONT 

I'll','  11 1  ’  ih'  I  1 1  ,ih  .  ,1  I  ',11 1 1  shin  ei  Mild  I  Hale  s\  slem  iisi-d  in 
I  he  dc'.i  1 1  pie  Ml  I'l  l  Im  v  n  i  n.d  w  rn  )  rnni  I  he  nrnim  is  eln  I'-en 
111  he  I  he  p.  mil  vv  In  I  e  1  lie  1 .1',  vv  ll  h  /|  b.  I'nll  I .  Il  4  s  I  lie  spile  I  e 
and  I  he  I  .ixe.  i .  p.n  dli  i  '  n  l  he  im  a  h  nl  rav.  I  ’nml'-  ( r .  v  i  mi 
the  virln.il  .'..ivi  h  ml  ,in  nhlamed  liv  extendi  ni  'hi  ):"n 
Miel  rn  .ill  V  'h  i ,  i  me  .ed  ,  iM  ei  i  d  rav  liaekward  .i  de  l  am  '  h 
I  mill  I  111  pi  nni  I  ;  mini  .  -  a  n  1 .  I "  I  Im  i  nurd  m.il  e  sv  i ,  m . 

I  lie  sill  till  e  I  ll  111!  1 1,1  III  ill  "1  I  III  plain  ,  4  si  al  le  r  Ini'  i'.  di'-i  nlied 
hv  I  he  pniiil  1  '  I  .\  hi  1 1'  /i  in  ■  h,  I '  •  In  b  I '  ' 

b  h,  .  .md  i;  I '  I  hi  i.idins  .4  I  he  Imlihle,  ( 'ah  nl.il  mnni  I  In 
slnpe  III  III  the  '.i  .n  I  ri  I'l  1  III  ,;i  Ves  the  I  nllnvvme  p.n  ,  i  iml  ra 


lor  t  h*  V'lrt  Mill 

".  .IV 1  1  rout : 
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♦  K'l  4^  IlUj  *'Tlcnn  «)i  Ih*'  i)<»fMm'|iii  lo»  l(>t:;Hin<)  CtXTftliaitFVi 

il  p(*'Ol  O  'M'  tjiv*  ')  viM'i  i'  vv  w*'  l(  Dill  Thi*  uppi'f  hnr  of  liMKjth 
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Wilft  Uio  tmi'f’b'  ■  •!  >  *'  ••  |!.i'  I  'l  1  ;h  iciil  ihv  point  ()  TIid  fot.it*'  1 

tlcOat*'  >'  ‘  »•*  t'vl  !u  lilt  iht.iU*  ll  Itiylui  riip. insKw*  t)|  thi' 

wav«  Iro”* 
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m  =-■  liim/i  -  (  us  (f)) 

'The  ('ulfulal  inn  nt  Ihc  dil  1 1  .n  l  nm  i  .  l.n  llll.llnt  l(V  lif  t  rilnii;; 
virtual  \vavi‘  Imnt'.  iisi:ii;  i  (inrlmaliv.  in.r)  wimli  have  llir 
sniJM'  ori|;in  Init  nri'  ml  alfd  <  III,  k\\ ,  n-  i)v  nti  .iiiiilt'  ( ;,■  ''d(  ■  •!<, 
S(i  that  the  |Kiail  isi-  r  axis  is  tin'  p.il  li  ul  l  he  si  al  li-rcil  rril  u'al 
rav  pridii  It'd  In  iii'imu'l  t  ic  ii|ni.  I'innnl  Imi  (I)  and  llin 
rotaliiin  (ranslnrniatiun  i;i\c  lla  lulInxMiii;  t'Xiircssinn  Inr  llir 
now  (■oiircli'ialos  nl  llto  viruial  xx.nr  rnnit 

u  =  X  sin'.’d,  (  \  (  o  '111, ,  (Cl) 

(■  -  — .T  cns'id,  )  X'  (7) 

i'o  obtain  a  I'ayli'r  oxiiansiun  ul  r(«)  alMiiii  u  o  i|„. 
lo  viii);  ilorivalivos  wore  (.nnipntod  tnr  II  -  H,  ilrulu  -  ;i  i; 
and  (/ 'V/i/n  '  7  (//I'l/C  (  /s/  d/i/  ' 'df(  wlnac  ;>((()  i/r/i/C 

and  (,'((()  =  till /till  Tin  laliulal  inn  ul  /■  and  ./  mod  l  ho  r<  !  i 
tnin  c  -  (I  Sind  xvl.h  h  IS  ovidoiiM rum  1 'i;!  1  I  ho  loini  ul  Iho 
Taylor  sorios  whioh  is  liiioar  in  ii  xaiiishos  hoian.-.o  ul  iho 
I'ourdiiiato  rulaliun,  tliis  is  ox  iili  111  an.  o  ,0''.  I  o  .mil  i/P',  1 
--  a  ousd,  Oniis-iun  ul  iiiliio  .n.d  hp'hur  uidi  r  li  ino,  in  Iho 
Taylor  sorio.s  );ivos 

LI  ■■  im'^  (HI 

(H  =  -  l/tl  city  I’,..  C.l) 

For  till'  case  ol'an  air  Inilihlo  in  walor,  nuiiiorio.d  1  ah  nl.iliuns 
show  that  h'.i).  (iSI  is  wii.hm  Hfi  ol  llio  vahio  ul  0(0  )  proda  iod 
by  l''.(|s.  ( I )  |7 )  pnn  idl'd  I  hill  0  •  n  *  iHIVn  In  1  ho  <  ah  ula 

lion  ol  the  diltrai  lion  that  lulluvxs,  h.ii  i.si  xoili  bo  lo.od  lor  tho 

rojtion  II  0.  (Airy's  tri'alinoni  ul  tho  raiiihow  also  roliiins 
only  tho  leadnin  torin  in  iippruxiniaiinn  tho  virtu. il  wave  truni 
but  in  that  caso  it  is  a  i  nhio  loriii. )  I  ho  sio.ii  ul'  n  mdioiiios  Hi  ii 
the  ronoctod  w  ixo  is  dnoryont 

II.  AMPLITUDE  OF  THE  VIRTUAL  WAVE  FRONT 

The  physical  optics  .ippmxiiiiatiun'  involvos:  ui  ilo 
specirication  ol  tho  pulari/aliun  ul  iho  inoidont  wavo.  and,  lii) 
the  uso  of  the  ri'lloctiun  ouoHo  ion's  lur  ji  p/mio  snrl.no  lu 
apiiroximalo  the  roriooiod  liclds  at  i.n  h  point  on  ji  oiirvod 
snrfiicc  In  the  applicalinn  ul  the  loohiiiijno  Hull  lulhav'-,  Iho 
inciilont  w/ivo  has  an  c loot  r  a  liold  pul.n  i/.il  mn  oil  hot  onl  irolv 
perpondicnlar  to  llio  pi. 1110  ul  ■.oallonni'  or  oiitirolv  par.dlol 
to  it;  a  subsoripl  j  will  ho  o'"  .ii;nod  ;  ho  v.ilno  t  in  1  Im  In  1  0.0  o 
and  2  in  tho  socmul  (.A  lino. 11  1  uiidnnaiiuii  u!  lollooiod  Holds 
may  bo  nsod  lur  uihor  nn  idoin  pul.in/aliuiis  )  I'u  u  ili  i.int 
obsorvor,  t  ho  rot  looted  w  a  \  e  appo.n  -.  t  u  i  uiuo  1 1  .an  .1  \  11 1  n.d 
wavo  ileserihod  In  Hie  lulluw  iio;  .  uni  pi  ox  .101(1111  ide  ai  (luinls 
on  the  iJ  axis; 

I  ’Hi.M  ~  rp"'  . .  MID 

whoro;  1  =  \  —  1  ,  lo' is  tho  anpuhir  I  roipiein  '  .  I  o;  Hio  i  nno,  i. 
=  2jr/X,  X  is  the  w,i\  olon;:l  11  ill  1  ho  on  1 1  1  'iii  dinni.  .nid  /  .  e.  Hie 
am|)litude  ronectiun  cuoll  u  nnl  ul  a  |  il  a  no  si  n  '.no  tor  Hit  .int;le 
of  inciilotn  0  d'’; 

ciisd  -  in  -  sin  ’ll) ' 

ci  ,  ,  ,  111* 

cusd  -I  1/1  -  .10  d)i 

-n ■  ' Cusd  -f  in  ’  ■111  I u  ' 

r  j  =  -  . —  '  ,  1 1 I 

n  ' cusd  )  I n  '  -  sill  d) I 

Konoctiun  lur  d  '  d.  whore /i,  -  /i  • .  n  o.m  In- dosoriboil  hv  I  In' 
t ulliwv ill);  .ipiiriixiiii.il lull  ■  I  h.it  six  O  1 1  .■  h  .idiii,;  di  p.  n  ii  I.-  O 


of  n  and  r  i  on  t  ~  0,  ~  II  >  0; 

ri  1  —  (8(/o  cusd,  )•''-,  (111) 

r.i  ~  -  1  +  ri-’(H</n  cosd,  ( 14) 

and  in  each  caso  l<ic/(/dl  ■  <  'i-i  Whon  d  >  d,  ,  Iho  ronoeliuii 
cuolfuionis  are  iininiuihihir  and  coni|'lex;  they  are  ~ 
ox|)(  -Ml])  .Old  rj  -  -oxi)! -m');.)  whor  "'': 

tan(Xi/'2)  =  (sill'd  -  ri  ')'iVoosd,  (l.'i) 

tun(X.i/2)  =  O'*  lati(i'>|/2).  (lb) 

Kollocli  111  lor  b  ’  /),  whoro  b  -  h,  «  11  can  bo  doscribod  by  tho 
lur'.wini;  aiiproxiniatiuns  that  doscrihe  the  loadiii);  do|)on 
don*  os  ul  ^  ]  and  <’)  ■  on  1'  ~  H  ~  II,  >0: 

Xi  "  (H(7)i  cusd,  )'  "•’,  (17) 

b't  ~  ,’j'''(H< 7ri  cosd,  )'  -,  (!Hi 

.\p|iruxiiiiiitiuiis  fur  r  ill  lorins  of  11  may  ho  ohtainod  iisni;; 
II  -  Il  '  ti/iiiii,  )  For  u  -f  I),  K(|s.  (17)  and  (18)  nivo 
li/ii '  '  '.vhoro 

til  ^  (;;//i(i)'''’/eusi), ,  (19) 

li'i-n^ti.  i.:i)i 

The  •-(le*  il  leat  mn  ul  '  '(1*  M  is  I  iirl  llor  smiiilifiod  by  takilo;  e  1 
/■  =  II  lur  d  ■,  0  _  iiinl  oipiivalenllv  lor  11  <  (';  tho  use  of  tin . 
ii|>pi  ■  ixiiii.il  mn  Is  nml  ivaled  hy  t  lie  di\  cri'oluos  |  dr /(III  \  and 
.■7,/'/..'' as  I  •  I  0  ,ind  ►  fO,  respectively 


111.  FRAUNHOFER  DIFFRACTION  PATTERN 

M  ho  seatterini;  aiii|)litude  at  lari;e  dislancos  from  the  hiiliblo 
is  proporiinii.il  (o  the  Fourier  transform  of.-]  ox|)( •-i/.'ik*  ■’): 

/,()))-■  r  e/K),"'"  ""'-'"'•(in,  ('21) 

4./  •  '• 

where  //  -  ...  'an  ancle  of  (he  obserx  .ilion  point  rolaiivo 

lu  )  Is  a  .slimed  (u  he  small  so  Iha.  ,^in))  ~  )/;  and  conslani 
mii;;iiiltido  mnlliph.alive  f.ietnrs  tii.il  procode  the  intecral 
have  heon  ninilled.  The  ai)|)ru\iimit ion  doscribod  in  llio 
prei  I'.lm;’  .ei  lions  yield  the  lulhiwili);  e\|)  ossiun  for  /, 

f;  -  i  (  e'  ’Md,],  pej) 

./Il 

wtioio,  hero  ,nid  in  the  eipialiuns  that  fmlmc,  Hio  U|ipor  si,;n 
is  lor  pul.iri/.il  lull  /  -  1  and  I  ho  lowor  mio 's  tor  j  --  2. 

Fv  ihi.iliuii  ul  I'!.!  ('.',’)  o  f  leilltalod  hv  a  ehaiice  of  v.iriahh's 
ihal  luiipilolo..  llio  (larli.il  .ipiaro  (urinoil  hy  Hio’’  (  i;u  1 
I  )(  lino  ,■  I  -  to  ,\  I  '  ’ll  .  /(  w  hole 

I*'  ■■  I) |((i /X  I  (  os(f  P  '  f '.!;!) 

SV*.  will  ai'pruximale  the  nilopral  onlv  for  0  <  ./'  •'  .A,  win*  h 
ctirrespoiiils  to  (he  re.'iuii  (iredietod  lo  ;;oumelric  uplles  |o  ho 
si  rune. I v  dhimiii.il ed;  the-  rest  net  inn  civos  ir  >0.  Kipial  iuii 
('2'2)  may  '  boll  ho  writ  ten  as  lullows:  /'^  -  i  (  Xo  cusd,  I ' ' 'c'/  '2. 

whoro 


-  ,,d  <1,111  III  I  *  t'/  i,.| 

(2.1) 

V'  .I.Mi  1  -  d/,’  '  ll,,  ,  ■ 

1 '  '  -  le  '  ’  ')(X  (1  onsd,  ) '  * 

t2.ai 
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In  the  cxiuiiK  iiliiil  tli.it  pri  iftlrs  tin-  iiiti  itriil,  i/'^(u’,(l)  is 
evaliiiltcd  liv  t lie  ii|)|iiii|inat('  varialilr  siilist it iit luii  in  lsi(.  t'd.^), 
(0  replaci's  n  iind  ir  ri  plai  cs  .•  I;  plivsn  alK  ,  v’ .  tfi'.l't  is  the 
appriixiniati'  iiliasc  shit  |  assuci.iti'd  wit  li  r,  i it  I  he  pisinictrir 
ray  scattered  te  tiu'  uIim'i  vat  am  aii!'li'  'I'he  seennd  term  in 
this  exponent  lal  arises  1 1  uni  l  he  ei  nnplel  ion  ol  t  lie  sipiare. 

VVe  now  ex.onine  tlie  eondiiion  for  oinittmp  the 
ex()lii/  ,  ( .'.a  )|  lai  lor  w  il  hni  the  ml  e;:  i  ,d;  when  tin.  siinphri 
cation  is  pernnssihle,  tlie  inlei;ral  inav  he  eslim.ited  iismi; 
Kresnel  mlepraU.  ||  I  is  a  sol  I  leient  Iv  slowlv  varvin;; 

fmu'lion  ol  .  1  he  prim  iple  eoiit  riliul  loo  of  t  he  mteprand  oe 
curs  when  |.’|  T  heeanse  the  phase  o|  i'xp|i  I  .a  j  issta 
t  iotiary  at  j  -  h  t  hor  a  hieid  diseussion  ol  stat  lonarv  phase 
aii|)roxini.it  ions,  see  Seael  'h  Koi  a  •  0,  c  ^  a  ,  and 
it  is  evident  lhal  |v'  ,1  .'.a  l]  i/  ,1'  I'h  eonseipieiil Iv  tlie  lai  tor 
in  ipiest  ion  inav  lie  omit  ted  il  the  i  .idnis  ol  I  he  I  in  hide  is  snt 
(iciently  l.irpe  1 1 'll  v' .  ( '.!,m  «  I,  l''or  ev.iniple'  waler  has  n 
1  d.'i.'t  and  \  isede  li|;hl  Iromalh  \e  ia'>i  r  ha-,  a  waveienpth 
in  waler  of  n  '  0  i;:t'es  (nn.  I''or  pol.n  i/al  ion  1  t he.  eondii  ion 
on  i/,(-."l  piles  ia,'ri',l /imd  '  >s  1  and  hn  polari'/.il nm it  pives 
l(l/,a‘IO  pnid  '  ■')  1  whole  I'.ns.  I  I'.n,  r.'tl),  and  (d.'il  have  :ilso 
heell  Used 

Alter  some  ro.n  t.nu'.ini;,  the  re  -ult  ol  this  apiiroxiniiil ion 
may  he  wrilleii  a-  lollnw- 


A’  ('d  . . ' 'i|/''la  )  -  J'l  ■  ■  »1,  Cdi;) 

where  /■'{ii  'l  is  hresiiei's  iiitepral'"' 


andf'l-'o)  =  .h'l  ■  e)  =•  Km  any  stale  of  iiu  ideiit  po 

lari/alion.  ihe  lollowinp  fai  tor  pi-es  the  dependeiiee  of  the 
intensity  on  >/  where  i;  ‘  e.  the  eoinplex  ooniiijmle  ol  f,': 

HI!*  =  |('(a  I  d  'i:\-  I  |.K(a  I  (  ‘/if.  (.’{(I) 

'The  iiorinali/at  ion  ii'.ed  in  l''i|s,  I  '.'  1 )  and  i  .1th  is  siieli  lhal  ):)i* 
•‘dasfi/i.  -  (,'d  lii'eoines  l,ii'',e  and  posit  IV  e  As  expeoled. 

*  0  as  i/i  •  i/i,  heeoines  l.iri'.e  ami  posii  iv  e,  however,  I  lie  rep  ion 
le  <  II  IS  mil  o|  Ihe  taupe  of  tin  s  al  id  il  v  I  esi  tor  'he  approx  i 
Illation  used  ni  deriv  inp  K.(|  (.im 

I  )avis'  aiiaivsis'  ol  l  he  peolnel  nr  .i  il  let  inp  inav  he  used  to 
iililaiii  the  iioi  in.ih/.il  a  111  ii'ipniod  to  expri- the  ale.ohile 
iiileilsily  /  as,oei,iled  with  the  o  llei  led  i.iv  at  a  distanee  />’ 
»  (I  (roin  tho  venter  ol  the  Inihhir; 


/  --  /„(a//i’ I  S.  Cll) 

where  hi  is  I  h"  mndelil  mleiisit  v  ami  pp'* /h  >  for  liirpe  n 
SO  1 1 -.at  l.l  I  )  c  mm  ides  wil  h  K.i|.  t  :i  I  ol  lie  I  1  III  t  he  apliro 
priale  hnul  Korl  herinm  e.  il  l!ir  mleili  leme  helweeii  Ihe 
tolallv  reflevleil  i.ivs  ,ind  Ihe  iiiolliply  relleiled  ravs  are 
oinil'ed,  the  Intal  iiileiisilv  /,  lievome., 

I,  /nhi'h'l  I  (, 

will' re  < f"h,'i)  IS  ihe  total  emit  rihiil  imi  hv  !  he  nnill  ipiv  tenevled 
rays.  l''(ir  an  .ir  hiilihle  in  water,  lihimin  ited  hv  tmpol.iri.'ed 
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HCj  'I  View  ol  Sir  .ipp.ii.ilie,  in  the  pi, me  of  svanoiinr)  wIikPi  w.is  hrio- 
eontal  in  llin  lat-  .i.iii  ry  Tlie  tioMorn  ol  the  tinhtilr..  w.is  in  i  onl.u.l  with  .i 
ii.aerlle  Iti.il  pievr'n'o'l  ihr  h-mhli  Ifiini  llij.iiiiu')  nut  ol  the  '.e.inef iii'i  pl.tne 
The  tori  ol  llie  hiitihlr  (whti  i.  o*  the  r.ule  rli--, played  tiere)  was  viewed  try  a 
miciosr.opii  Pie  w.iyelemiih  m  an  ol  the  illumination  was  bJ2  H  nni 


I  ip  hi ,  'I 'a  hie  III  ol  lie  I  1  shows  lhal  fi"  is  a  siiiool  h  liitiel  ion 
III  pi  lor  1, 'I  m-ai  .  lOil  ih.il  =  II.IIP'.’T  and  (i 

(7.h“ )  -  n  no  ti  ( 'onseip lent ly  ( i'  for  an  air  Inihlile  in  Water 
IS  apiiroxiiihiieh,  1',  lo  :  i|  I  he  asvnipliit  ie  value  nf  V  S 

Kipi.ilimi  '  .0)1  .omi.ii  in  1. 1|  HI  lo  1  he  well-kiiowo"' I'di'siiel 
(neiir  /le/r/l  dilli.n  '  "ii  hv  a  si  i  op, hi  -edpe  svhere  I  he  parijiii 
elei  a  i  .,  lo  lli.ii  .  i  .r,  | iro| im I u nial  lo  ihe  lateral  disidaieineol 
ol  the  oh- er\ .0  en  |ion.i  lo  in  ihe  edp.e  of  the  p.eoinetfn 
shadow  Kor  hot  h  I  a-i  in  I  la' repimi  ii  >0th',  ioleiisil  v  p. 
predieled  (o  o  .i  dial  i  with  an  me  ilia  I  Ion  aiiipl.llmle  1  hat  lie 
creases  r.ipidiv  with  nu  rr,e  inp  ii  .  Tins  hehavior  is  predn'led 
fnr  Ihe  liv  h,  In  ;,iii  m.i'v  ot  a  hiihhh'  hee  nise  nf  the  alirupl 
variiilimi  in  the  .la  phindr  of  ih(>  virinal  wave  Itoiil  ne.ir  the 
erilical  .'iiij’le  .oid  hei  ,nme  ol  t  he  .  .'/r-rnli/M'  of  t  his  vvave  I  root 
indoeed  hv  I  hr  rrll...  t  loo  I  rmn  I  lie  rorv  ed  surfave. 

Krnin  llir  aforrinrni  lo,  ril  hound  o|i  f- 1  for  llir  prinviplr 
ennt  III  111  I  loir.  Ill  I  hr  iiilr,'  I  ,il  10  r,M  I,  wr  i  an  est  imal  r  h  n 
th"  II'  ri'pion  of  ml  r  rest ,  t  he  errors  o'l  rod  need  in  I  he  ea  ho  la 
lit. 11  III  the  'i,  Irmn  Ki|  ,  (  10)  and  t.'O),  |n  the  interval  0  ■  n 
<  It.  K.(|  (hto  pirdii  I  .  two  inlensily  niaxima.  In  tliis  rePimi 
the  pi  ineiple  roiil  liho' Ions  oernr  wit  h  II  ■'  O  "i  l/|„,o  -I'  I 
ie)l-  Aid)'  'with  a'  -  It.  Wil  11  t  he  all  iieinent  iom  d 
and  wllh  II  1  le  n.  the.  eslini.ile  'lives  Om.in  *on  ,\ 

coin  pa  I  is(  III  111  K.i|  i  i  i  s\  it  h  'Z  i  ealeiil.il  ml  iisinp  l''.(|.  1 1  Oi  and 
this  shows  lhal  I  hi  in  I  'r  ini  i  odoi  ed  in  r>i  is  less  ihan  0'  r  A 
similar  emnp.o  1  mi  o|  !■  .  |  .  '  I  .S|  and  t  ./d  I  pive  a  hound  Im  i  he 
erroi  in  '  ■  "I  .s' .  < ''iie-riin.  ni  Iv  I  he  use  o|  a,  i  p,  jns 

I  died  Koi  I  In  1  ir.  .1 1'  on  r  (i ,  the  0|)per  lino  I  ol  ill 

tepralimi  III  Ko  I'.'i  1  e  n  iM  e  w  r  .i  ii;  el  ol  appn  ixniial  imi,  I'he 
m.iiii  |o  1 1!  n  ,1  n  .|  i  h 'hi  ippi  ox  mi  il  n.'  .s  desi  r  i  hrd  ili  Sri-., 

II  III  e.  III!  .mno.riiin  yvilh  rx|ii.|  iiiiellts  drsirdied  ill  I  lie 
seel  Ion  I  h.O  1 .  -It.  i’.v 

IV.  OBSrRVLD  fiCATTrRING 

I  he.  s( ,  I  ion  .  h  .|du  phol .  .p  i  ,iph  ii'  ol  iserv  :ll  imis  ol  ihe 
er  1 1  leal  aoph  v  .oi.il ..  ais  m  l  lie  '-1  ,it  1 1  reil  ml  rie-il  V  pi  edn  ted 
in  I  he  pri  V  I.  lie  i  i  l  ;mi  1  i pm  e  :  1  -hi  iw  s  ,mI  i.ipr.mi  ol  I  lie  ,(|i 
pn.itus  111  lilleil  wall  r  was  plaeed  m  an  ahmmimn  eon 
laiiier  1  he  mm  i  w.ill  .  ol  Ihe  l  oiil.iiner  li.iil  heeii  hlaeki  iieil 
So  that  I h.  \  .ill  m  III  li  mill  h  ol  the  si  r.iv  lip, lit ,  A  steel  needle 
[lassed  vei'i  1.  dh.  I  hrme'li  .i  riihher  sepluni  in  the  Iml  imo  ol 
the  emit  111  II  I  I  he  ol  In  r  end  i  >|  the  in  I'dle  was  ei  ii  im  s  led  In 
a  sv  ’  inp,e  t li.il  i  mO ,oi inl  air  1 1 1 ir i/m\t al  holes  1.0  eni  m  d i 
.iiiiet er  pa-  I  I i  I  n I  r m.h  I  hr  w.dl .  ol  I  he  em 0 .iiiii  r  al  riphi  .m 


.l.t)|)l  Sue.  Am,,  \  of  tpl.  Nil  0,  Septeinher  tf)7*,l 


I  liilip  1..  Marstmi 
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30  millr  .  <in  If, 


FIG  4,  F^iritoyrtiph  ijf  IIhj  rntii.,)!  .nuiUW.tt  dilh.-c.  Inhi  ,i  •■•n-  .• 
baildji  hiivft  ;i  chor.irlPtisln  sf}p;iriili(Mi  nl  ;ippt(txitri.iPtr,  'tit.w*  *l'-  v  n.  ii- 
polnri/ocl  wilh  its  ckicltii;  (uil.l  p.ii.ilk'l  I'l  llif  S'  .iM.'iiiu!  pl.i'U 


iMi.j  iii  i  KMKd;,  (loiii  Idtl  to  fiytil  The  verllril 

.|■ll  ..ii-  . . .  liy  tlilldK  lion  Tltfl  incident  beiint  w;is 


'I'ln'sc  ltdic.s  wcri'  St 'll  I  I'd  ys'il  h  ( )  nitt;  •  to  in'i  .  cnl  n.it  > 
Iciiks  and  wcri'  tus  I'li  il  liy  oiil  ii'tdli'i  i  ii  ni  i '  i.ilfd  iiid  nv 
liillililc  (iiilld  III' l■l'('llH'd  III  lilt'  Hal  li|'"l  I  111-  HI  I'llli'  l'\  di 
iTi'iisinp'  the  v'liliiiiii'  III  till'  ''>'riii|',i'.  lliis  inil'ldi'  w.i  ■  staid' 
priividi'il  its  viiliiiiic  was  l<t'|il  snlTiririillv  mall  TIh  liiild'l.- 
was  a|)|)riixiiiiali'l\'  sidii'iiial  in  it-  ii'iiii  i'  I'laop  'I'mi  .■ 

(lisciissioii  III'  I  hi'  stall!  Ill  V  iiT  I  mill  ill",  ai- I'l 'lid  I 'll  1 1 . .  '  d'l  ■ 

s(‘l'  Hi' I’.  II.)  rill'  I  nisi'  III  t  III'  I  HI  I  dill'  wa  ■  in  i  mil  .i  i  il  li  I  h- 
iii'cdli'  diirin).:  Ilir  nliscrs  al  imi 

'I'lii'  |)iilafi/i'd  lii'ain  Iriiiii  a  li  iiiW  I  In  Ni  I  lo' i  v-.i-  .  > 
paniU'd  In  prndiu'c  an  S  inin  diaiii  paralli  I  Imai  i  riii-  I'X 
pandi'd  licam  passed  tinrnialls  ami  Inn  i/niilalK  '■  liiinii'Ji  uin 
I  if  I  he  windows  in  snrh  a  svav  that  I  hr  i  niit  r.il  r.iv  s  ol  i  hr  In  am 
pa.ssed  cliise  In  I  he  i  nnl  er  nf  I  hn  laihhln,  A  i  ainni ii wed  I  hn 
hlihlile  lliriinph  the  sci  inal  wiiidnw  I'hi'  nneiiii'  .nnl  ihi-  nnii'i 
wall  III'  the  altitiiinnia  i  niilainer  vsei'n  lilanki'iii'd  in  ii  'hii'i'  iln 
intensil  S' III  si  |■a>  lipjil  lhai  wmdd  ha'.niill'  rw  .  ,n  nh  i  iinil 
!  he  seal  lerinp  iiT  I  he  I  in  hide,  I  ,ii:hl  seal  I  in  ed  liv  I  In-  liiihh.’i 
near  I  he  eril  ieal  seal  leriiip  aiipde  had  an  .nii; In  ol  im  aimn  •  ai 

the  water  inwdiiss' iiilei'laee  III  ap|iinMiii,ili'l\ 'll!" 

’I'liis  li^'hl  is  ri'I'raeled  al  the  ssiiidnss  nr  niii  il.ien  I  hi 
eaniera  lelephnln  lens  iNikkor  inanulaet  ni  I'l  hs  N'lkmi 
I  lie.  I  had  an  el  leel  i  ve  h  nal  li'ni;l  h  /  i  d  'i  ill  m  m  I  hi-  leii.  u  ,i 
placed  wilh  its  axis  hiiri/niilal  and  ap!  II  iiMiiiaii'K  p.ii  .dll  I  1 1  ■ 

I  he  rerr.'ieli'd  i  ril  ieal  I  ;i\'s  I'hi  . . '  'In  im  idin' 

iniiniiehrmiialie  liphi  was  leaialls'  i  hm.rii  isiiii  llii-  i  lei  lie,  lii  ld 
parallel  In  I  he  seal  lennp  plane  l  h.il  1 1  a  n  spinel  i  ■ .  ■  '  m  lie 

preei'dini'  analsse.  I'hi-.  i  hniee  h  n  l  he  pi  l.n  i/.n  imi  ip-e...  i  In 
t;reale.sl  inapiiiliidi  Inr  the  rnel  I  n  mnl  nl  i'  oi  I'.p  'lli 
( 'nnseipii'lll  Is  ll  pises  the  limsl  lapld  deeiease  tin  |i,|  with 
inereasitii;  '  and  il  shmild  la  ihr  pn!ari/at  inn  tin  vslneh  the 
iiinissinii  III'  I  he  ssave  In  ml  wit  hi  i  ■■'Uin  l'',i|  i  I  svill  pinilni  e 
the  sni.dlest  errnr  in  I  he  pin  lii  i  "d  se.-ii  lerinp  I'he  diaiiiei  ei 
III  the  linhhle  in  ihi  I  rnm:  /i/im.  e,,i  n  ".iiied  wiih  a 
inieriisenpe  Inealed  in  the  air  diieells'  ahnse  the  I'ldiiili 

I’liiilntiraphs  were  lakeii  svilh  the  eaniera  leie.  Im  ii.eil  i" 


i  l.nii-  ■  'h  a  'h,  niin  plane  ss.e.  al  the  nplieal  Iransliirni 

j  I  a  I  in  II  in.  I, ell)  '  i,i  si  here  p' and  i/i,  '  deiliile 

I  a'  lie  .'  ,  Mie  I.'.', II  II  side  the  I'l'll  1  hat  enrrcspilllll  III  i/i  HIlH 

■a  '  he  .in  1  I  'l  I'l  -I  ,|"i'ni  !s  niei'enieius  in  I  he  dislanee  / 
l.n.  ip'.el  1.1. ','11..  (Mlise  .lie  I'l'laled  In  ihe  iliereinellts 

■  ■  h  '  .\  /'iiii  \/ "  /.  Nnnieneid  caleulal  inn  III  I  he 

.1'  ..I.  I  i;l  I  •  .ind  M  e  a  iler  air  hilerl'iiees  shiiw 

II  Ill  111  a  III  e.  hi  relaii  d  In  I  he  I'lirrespiindint’ Al)' 

h'.  'll  ".Nil  S',  nil  a  ni'i;  I  II!  ihle  eii'in.  Tile  sileeess  id  this 

1.,.  I  h... i  ,ri  Ii'""h  lie  ehnii  e  nl  ihe  plane  nf  ilie  s'iewed  windosv 

•  .  1 1  I  h.ii  1 1  I ,  nr.ii  Is  per  I  lend  nil  I  111'  In  I  he  eril  ieal  rny.  'I'liest* 

1 1  ..I  ill  s  i.'i  IS  III  - 1  annl  n  Ill'll  II  1  ips'e  A  i)  -  A//ii  /'  [Ol  eourse,  if 
ihi  1  nli:,  II  I'll. I'  esieinal  In  ihe  sphere  had  hefn  iinifnrm  sn 
'  1 1.1 1  111 .  s^  "idi  ns  S',  .e  needed,  i  I  seniiii  I  lint  lie  heei'ssary  Inin 
\..i  ll"  ■  I  ,i|.|ii.i\iinali"n.  I'hen  .\l  nh'  eiaild  lie  nsi'd  as  the 

ii  I'l  l "11  p.ii.innier  in  l'a|,  r.’li  lieeanse  id'  the  Kiiiirier 

I  I  1 1  I  "I . I !  1  n  I  n  M'l  I  le  ,  1  ll  a  lens, '  I 

li  111  ll  .1  I  I  "ll  is  e  I  I'pi'i  iihii  I  ll  111  III  a  phiiinpraiih  III' I  he 

I  ll  11  ll  .ineli  -I  ll  ll  iiin;  ss  ll  h  a  hidihle  dianieler  Ilf  I  f)  4.  ().l).''i 
■lO'i  I  he  ■  .ll"  mihi  .III  s  I  Al)  lelai  is'e  In  an  arhil  rarv  refer- 
...  e  M  I  h..  him  n  a  d  ss  ,is  '|'n  \  svil  h  an  expnsnre  time  of 
1  I.. I  I  he  III . I  III  s  ei  t  ll  al  d.ii  k  hands  ss'il  h  a  eharaelerislie 
p,i.  in:  I  ipm  .  sinialels  nirad  1  ^  | ,.)“  |  are  examples  id 
III.  ini'  le  O',  ii.inini  1  inedii  led  hs' !''.i|.  (  r.!)  'f'he  ei  ilU'Cnl  ric 
I  111  'll  I '  1 1  in;  I"  I  all  .1 1 1 '  i  eiit  ered  i  li  ise  I  ii  I  he  iiiosl  eeni  rill  dark 
! I  1'  1' i  .1 1 .  .ipp  n  eiii  Is  ,in  .11 '  1 1  al  l  assiM  ialedss’ilhspiirinnsre- 

I  h-,  I '.  n  I  al  liei  I  rnm  I  he  I'S  It  Ss  iiidnss'  III'  In  I  111  I  III'  lelepliolii 

II  II  Ihi  link  rie.iiin  nn  ihe  lei  I  nl  the  liinire  IS  i  11(1  ieal  is'e  of 

lie  I  e.  1 1 1,  I M  III  in  I  in  riisi  I  \  I  hal  i  n  en  rs  sshen  le  <  I). 

Ml  1  nil  leaiiiini  11  1  eaie.  ssere  I'hlained  frntii  Iss'ii  ncpatlS'es 
1  .  i.i.  il'i.ili  ,1  ipi.ml  ilal  IS  e  eiiiiipariMm  wilh  llie  Iheiiry. 
I  1,  .  ■  1  onip.ui  s  I  111  relal  ise  inlensil  y  inferred  from  a  scan 

■  I  I  111  I'.il  I  "  Pir  .1  Inilddi  d l.ilnel er  1 1(  I  .ll.'i  ,1  (I.O.a  min.  ’f'wii 

liei  p.i'  iiiieiei.,  ssi  .(■  mi  ll  Im  I'lii  ll  nepatis'e  In  make  the 

|i  .1  ■  .  .11  ll  h  I  hr  I  hem  \  pnssihle:  I  i  I  The  liieal  inn  ef  the 

. . I"  .  ii  d  '  emild  nut  he  di'lerniined  exiierimentally 
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I'hihi'  I.  .M.irslmi  liiini 


^  milhro4ion» 


(fomn  mi<;ro()unsiHiiiwlt't  sc, in  nl  ihn  (itinti  |(.i|iIik  iu'imIiui'  III,'  iMililili- 
diamelei  W>is  I  60  f'.n  Di,,  solid  I  imvc  i  .  mvii  i.v  I  d  i  'ii|  which  r-  ihi- 
fesull  ol  the  physii  .il  '■  ;,tiL'i  .ipproonviiioii  I  wi,  (  .u.iitiidnrs  ,id|usUHl 
to  (iicilil.ilo  Itiis  (,oii'|!  iii'.nn 


SO  ihilt  an  iirliitriii  V  I'liiist  dill  dni  ji'  was  .uiifiHi  to  I  tih  dlisi  issa 
of  ihr  mi'd  an  I'd  iiilcnsit  s  in  ol  Hr  i, ,  n|i(  ino/r  i  he  iip.mMiiriil 
with  thi’  llu'orv.  1  h  I  ’I  hr  ali-.ulnlr  inlrnsil  s  rinild  not  hr  dr 
termined  pxprrinu'r.tdllv  >n  ih.it  th  rx|irriniriii.il  iniriisitv 
/,.  was  olilaiiK'd  liy  innil  i|il\  ini',  inrasuirinrnK  hv  :i  roiisl.int 
factor  K  that  was  scirrird  lo  o|itnin,'r  ai'.imiirnl  with  Kii 
(dO).  'I'hf  inaundndr  uf  (!'  was  sul fii  iriil Iv  sinall  llial  il  was 
neulcctcd  in  tin'  cnnipulalirn  ol  ihr  Ihroirtn  ,d  prediction  ai 
that  cinild  he  used  as  the  prrdielioii  I'he  relationship 
between  tlie  measured  (ransiinllaiu c  of  the  nepaiive  <  and  /, 
is''': 

/,  =  K(i  -  r,,''') 

where  to  was  inea.snred  at  an  iinexposed  region  ol  the  nei’ative 
and  the  eontrast  index  >  was  eHlimaled  Ironi  niainilaetiirer’s 
literature  to  he  ('.70  for  tlie  conditions  of  Idin  (levelopincnl. 
'I’he  r(r;)  ineasurenienl.s  were  sniootlieil  hv  li.nid  lo  redme  the 
mnunilndc  of  the  eirculai  Irinc.e  artifaet  in  Ihv  iiilened 
'I'he  scale  of  the  smooth nij'  onK'  si e; In  Iv  exceeded  I  lie  pa'  nip 
of  those  fringes.  ’I'lie  valid il  s  ol  I  he  [iroeedine  was  suhsian 
liated  by  eomp.nni;',  paiaiiri  hul  ihsliiu  se.'Uis  Imm  ihesame 
negative. 

Kip  lire  ■')  illusir.iles  .1  I  lose  .e'lrnaeol  ol  ihr  rrlalivr  po-i 
I  ions  of  the  maxima  a lul  mini  iiia  rl  l  hr  inlnisii  \  hands,  w  il  h 
those  prediiled  hv  l'',i|,  lain.  I  Ic'  n  vahirs  wrir  raleiilaled 
usinp  l''.i|s.  I'.')  and  'I'l'ilh  II  I  mil  .''i  I  o-i  inin; 

('(le)  and  .S'i'i'l  were  idil, lined  liom  lal'lr-  I'he  va.Iui", 
and  loealioiis  lor  the  ni.ixlina  and  nnnan.i  wric  ohlanied  lium 
Kresnel's  oripinal  i  .ihni.il  loii  1 1  a  i  he  de  >  i  ipl  nni  ol  eder  dil 
fraelioii  las  ipioleil  hv  Nhi-'  .iir'l  'I'hi  ir  i.s  an  .i|ipairiil 
diserepain  y  ol  I  his  l  hem  ,  ■>'.  il  h  l  h,  l"llr\\  nr  Iral  i nr.-  oi  |  hr 
nieasiiremeni :  lil  there  appear;  |r  hr  a  diilrrrin  r  ol  inlni 
sities  ill  ihr  rrpiiiii  a  :  1  Ih.il  raiinm  I"'  r\|i|.inirii  hv  ihr 

omission  ol  ( i"  I  n  nil  ihr  ea  Iridal  e  a  i;  In'  lie  aiM|i|il  udr  ol  thr 
ilUelisit  V  osrill.iMm Is  exrrrd  .  I  hr  I  a  rde  i  loii  -  h'l  n  ’  and, 
( iii)  I  here  is  ,in  n  in 'X  plained  dr"|  i  in  ihr  mlrnsil  v  lir.ir  ir  a 
Allot  her  nri',.il  i  \  r  shows  ipiahl  al  lx  ''h,  ii.id,M  disaeacriarnls 
lint  also  a  siind.ir  a'.'rrrnieni  lor  ihr  h  ,■  al  eat  ol  ihr  m.i\iiii,i 
and  Miiniina  with  ihr  ili.'ax  l  ie  ';o.,aiiiv  :>iel  ipi  ililx  ol 
mra'airenirnt  s  I  h.n  i  an  hr  ol  a  aim'd  w  il  h  1 1 ,,  ix  .ol.ihlr  app.i 
rat  ns  arr  liol  adn  pial  r  I  or  del  rn  miei I  h,  ,  h  I  iriri  e  a  ,  ol  i  his 
eah  iilaliiin. 


A  conipiirismi  ol  photopraphs  taken  with  bubbles  of  dif- 
ferinp  radii  shows  that  the  anpular  spacinp  of  the  inletisii  v 
oscillation  del  leases  wiih  iiu  reasinp  n,  which  isex|R'cted  frotu 
the  lorm  ol  Kq  ('.'  !)  'I'his  decrease  is  also  evident  to  the  eye 
if  the  radius  e,  nu  reased  as  the  hubhle  is  viewed  thr()ii„h  a 
telesr  ope  loi  usi  d  on  ml  mil  V  Intensity  variations  similar  lo 
thope  (lisplaved  in  I'  ips  I  and  ft  were  also  observed  with  the 
ineideiU  [lolari.'.ilion  ot;  -  1. 

V.  FINE  STRUCTURE  IN  THE  SCATTERING 

('arelnl  exainin,il  ion  ot  the  photopraphic  nepalives  revealed, 
in  sonie  l  a  i  s,  ,m  om  dial  ion  in  the  inleiisity  pcriislie  in  r(  wlm  h 
was  .siipenmpo'.r  1  on  I  he  diffraetion  sluieture  descrihed  m 
Sec.  I\'  '  '  Kor  e\,imiilr,  oM'illations  in  iniensity  with  a  period 
ol  app;,,'  'll.  o  I  I  :  ml, id  are  visible  in  the  oripinal  jtholo 
p.raphi,  ),imi  .  ,  ■  lie  rd  ni  Kip,  I.  I  It  is  unlikely  I  hat  I  h  is 
'.I I  u,  I  MO  X  II  ■'  i'  -  III  1  hr  priiiled  xersion  appe.irinp  in  I  his 
|oiirn.d  III, I  It  I  oix  I  d  ill  I'dr, )  I'lxamii'al  ion  o|  I  hr 

nepal",!'  I,  ■.  lii.i  '  Ir  1 1 'oj'.r  spaeiiip  i;  decreases  xvil  h  in 
rrea'  II,,  ',  i  •'  ■  '  ii.  in , ,  i ,  ,sr,  i| le  XX  as  used  lo  measure  i) 

tr.,m  till  O';'  ',x,  i  n  i  ,ohhle  ■■  with  1 .7  liim  C  ?u  ‘d.'J  nun. 

d'h,'  o  'di  0,  l.n  I  iimiil  he  tilled  Ijv  the  equiition;  i;  ~ 
\i,',  ,  1.  t  '  '  I  r  IM,  however,  il  was  (lillieilll  to  ol) 

lam  loll  I  n  ■  I . .  Ilia  111'- I  roin  each  nepal  ive  h'jeaiise  ol 

ii,iia  all, I 'I  ■  '  1  1 1 .1' I  111  I  he  oseilhitions,  Ivnerpy  assi  > 

eiai  rd  X'  1 1 1;  ■ .  I '  no  i  ix  ,  iiirnlenl  with  h  <  I),  will  event  nail  v 
leaxi  if,  '  'i,!  '  ,hi,-.  1,'d  xxuh  i/i i,''. ,  Krom  the  linearily  ol 
I  he  1  qua  1 1-  I":  oel  i  hr  inapnitiide  of  A,  d  is  likely  that  the 

osedlaia,.  in  ,aii  ,d  hx  I  he  inlerfereiiee  (d‘ lhat  liphl  wil  h 
the  1 1  n  ir.ill ,  r ,  :  I,'  ' ,  q  1  r'.hi ,  liowex  er,  the  number  of  reflee- 
tune,  and  '  X  laii  i  .'  I  ,,r  p.iih  of  this  mterfeiin!;  liplu  lias  imI 

heeii  ,l'  ■■e  i',.  I  I  ;,'  i  | ■  iiat ion  rl  I h"  nuiioi  virl ual  sinii te; 

(ol  1 .  I  \ 

fiiiad.o  .  1  !  I  ’  .  r  e  '  d  osi  dial  ions  can  he  oli',erved  in  the 
riniihoxx  Jioiii  hqiiid  drops  xvilh  n,  >  rl,,  and 

inonorhroi.,  ,'  r  d  a  ru  m  ■lain,  'I'lie  uu  ill  or  li.is  round"'  I  ha  I 
I  licit  ■  p.ii  r,;;  ,1  ,1  ,  ri  ',rly  pro))orl  ioiial  to  (1  and  for  drops 

ofxyleii,  ■  I  '  n  111 '.x.iier;  U.ii.  ’I'hese  friiipes  were 

found  I'l  h,  'I  !  I  'h'leri  nip  miermioiniplitude  shape  o.-o 
eilhilioii  .  ,1  r  "ir  ■  i,  I  liameti  drops  since  ellaiipes  m  ihe 
drop's  liipr  1  i,ri  m  mpiihir  slnfl  in  ill,  posit  ioli  of  these 
lrm,',r."  ■  ,,o  ,■  I' a  111  h,  XX  hen  I  he  shape  of  I  hr  drop  USX  il 

1, Ur, 1,1,  I  111,11  ■  ai  ',1  III  a  I'rxx  I  b,  the  trinp.rs  appi  an  d 
to  lilm  1,  •  'r  I  ,  ,  'I  ■  d  I  i  •  dial  ions  111  the  shape  ol  l  he  hulihlr 
shonld  '  '  ■  ■  ,  III!  i  mr  si  r  III  I  lire  lo  hhir,  India  d, 

hiihh!'  l.ii,.  .  ai,,'i,  a  I  ,  lie  I  hv  room  X  iliral  lohs  mav  hax  r 

him  1 1  ■  I  I  i"  ' '  ■  .  I  1 1 II  ■  r  nepal  ix  es  xvhe  la  nniie  xvere  i  i , 

■le'x, 'I  I,  v|„  ,,  ,  '"  i  A,  Ir  I pii  idlv  1 ,111)  si,  I'lir  Mmil,ii  It  II  . 

•  ui'i  diM, . :  laii  li  aiiplr  sr.il  I  el  a  lep  xxilh  |■,llnlMl\\ 

SI  i!  I  m  ■  ,  ,1  .  ■  ,  o  .  ,  'ni|iarin;;  Kip.  -I  xvilh  pholnpr.qih.- 

ol  I  hi  lai  i '  a  '  1 '  n  ;  i  ,  1 1  ii  I  I  i 


VI.  DI?rUi'''SlflN 

'I 'hr  ,  I  at  n  .d  II  ',lr  si'.iMi  r  in;;  of  ripht  I  lU  hides  xvas  phoh . 

praphed.  ihe  'niliMr  ih  .  In  .  xxere  in  the  raiipe  l.li  '-Ml  iinn 

I'l.ii  h  ph"  i.'i  1 1 'll  1  •  inid.ir  III  appear  anee  lo  Ki;;.  1  in  I  h.il  Ihr 
rnipm  •,  'm  ml  'r  -i  I  x  '  n  a  hi  the  nil  n  al  anple  is  '  pii,  e  xisihle  I ,  'i 
rilhri  I, Or  rl  I'l,  I'!'  rl  | " 'l,i  I  i/.il  loii .  f'ollsrqnriil  Ix  iheir  n 
I  ,  o  1 1  il  ' " I  ■  "  I  r  n  1 , 1  I  hr  pi  l'd  u  l  ion  ol  pri iiiu'l  l  ie  I 'pl  a  - 
nr.o  |i  ,  .  ii’i  ',  '■  il  I'lr  ,m  Imhldes  ol  llips  ire  ill  xx.ih a 


I.MI) 


P.’UI 


.1  Dpi  Sim  Ami  .  \’(i1,  li'l,  N"  a, 'srpn  adirr  P.IVtl 
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I’hilip  I,.  MarsPiii 


FurtlK-rniDre,  the  hkm  iiriTiK'iUs  ol  llic  ixisitions  nl  the  in¬ 
tensity  niiuiiiia  .111(1  niiniina,  inc.isdred  relative  to  an  mule 
tenniiu'd  aMj;le,  ,iic  m  teasnnalile  agreement  with  1m|.  (dtit 
An  exiK’rimental  test  ol  icvei  .il  aapec  tsnt  I  his  ('alenlation  iiiav 
re(|nire  the  use  ol  phot  oeli'i  1 1  ii  det  is  t  ■  u  s  hei  atise  ot  limit  a 
tions  in  the  present  .ippai.il  ii  These  iiu  hide  t he  I'ollnwinj; 
features:  (i)  the  ainplitiide  ol  the  rin;;mn  "•  die  intensity,  lid 
the  ahsohilc  positions  in  the  maxima  and  minima,  (in)  the  rale 
(if  decay  ol  the  inleiisiiy  lor  le  5  1,  (iv)  the  rnnp  •  ol  radii  tor 
which  H(|  CIO)  may  he  used,  and  (v)  (|n;mlil.i"  ■•.i-ore 
menls  of  the  interlerenee  phenomena  deserihe  .  i  see.  \ 

I InforUmately,  the  l.irneiiess  of  the  hohhie  radii  make:'  it 
impractical  to  calciilale  the  Mie  seallerni); 

A  brief  diseussion  of  this  use  id  the  physical  oiilu  s  ap¬ 
proximation  IS  in  order.  Surlaee  wave  phenoiiiena'  do  not 
appear  to  he  included  in  this  eiileulation;  Kii;,  is  evidence  ih.ii 
they  are  not  needed  lor  an  approxim.ile  deseriplion  id  these 
ohservalions.  .slurlaee  w.ives  can  he  iiiehided  in  tin-  W.ilson 
transformation  method  whii  h  has  heen  used  to  lormiilale  the 
scatterinu  from  i  ((ivity.'"^ '  '  eomparison  with  the  result-. 
in  Ki)!.  has  not  heen  perlormed  hee.ius;*  ol  I  he  i  iimplii'aled 
nature  of  predieled  seallerinc  am(ililude  je  .','  Uef  l.H,  lv|s. 
(■t.9‘d)  and  (1. 101 1];  indeed,  the  present  aul lu>i  h  as  heen  uiiahle 
to  determine  il  the  rmnia;:  is  eonsisleni  with  those  amplitudes. 
A  second  shorleoininj;  id  the  approxmial ions  presented  m 
Secs.  I  III  is  that  iT.e  eurvalure  of  the  huhhie  orl houonal  to 
the  plane  of  scaitei  iiij;  is  omitled  until  1‘A|.  (dl)  where  its 
Keotnelrie  effect  is  included  in  the  o  '  factor. 

Diffraction  m.iy  have  a  major  effect  on  the  far  field  .ip 
pearanee  of  the  erit  ie.il  anple  s.al  terinj;  for  till  hides  in  water 
illuminated  liy  i  ollini.ded  while  liuhl  I.el  i/i, '  aiid  i/c  "  denote 
respectively  i/i,  tor  violet  and  lor  red  liuht.  The  dispersion  of 
water  and  l'',(|s.  1 1 )  and  1 1!)  j;ive;  <i\  '  -  <;>,  "  ~  17  mrad;  con¬ 
sequently  the  t;e(imetrical  sciitleriii)'  (Fiy.  IHof  Hef.  I ),  as  seen 
through  a  telescope,  should  appear  lihdsh  in  the  ret;ion  i/v " 
<0  <  0,.'.  It  is  evident,  however,  from  Fig,  .'i  that  diffraction 
smearing  of  the  othersvise  ahriiiit  change  in  intensity  (or 
monochromalie  illuniinalioii  will  lend  l>.  reduce  signifieantiv 
th'  coloration  predieled  iiv  geometric  opt  ics 

'This  work  was  motivated  hv  an  inx'  siigation  id  optical 
methods  for  measuring  the  refrai  live  index  id  superheated 
liquid  drops  suspended  in  immiscihle  host  liipiiils.  At  siiffi 
cicntly  high  lemperal iires,  such  drops  frequentiv  have  n,  < 
n,i  and  they  should  exhihit  critical  angle  scalieriiig  phenomena 
similar  to  that  deserihed  for  huldiles  Tixeii  though  the  dis 
tinetive  critical  .ingle  u  allering  predirled  liv  geometric  optics 
is  altered  hy  diffraction,  rt  may  he  iiilerred  from  the  ahsolule 


angles  of  the  inlensily  maxima  and  minima,  provided  the 
present  theory  is  louiid  to  give  I  he  (ih.xo/idc  angles  wlnm  n  is 
known. 
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Scattering  by  bubble  in  water  near  the  critical  angle:  inter¬ 
ference  effects  (P.  L.  Marston  and  D.  L.  Kingsbury)  J.  Opt. 
Soc.  Am.  71,  192-196  (1081)  [An  Erratum,  reproduced  below, 
was  published  in  0.  Opt.  Soc.  Am.  71,  917  (1981).  This 
notes  errors  of  transcription  whiclT'are  not  corrected  in 
the  reprint  given  here.] 


errata 


In  the  caption  to  Fig.  3  of  our  roper,'  »)■!/>-  (/>,.  should  be 
written  »)  =  </),  -  0,  svhich  is  the  definition  of  >j  given  in  the  text 
of  our  paper.  In  the  line  below  Kep  (If)),  ti,  »  (8/nol'^'^  cos  0^ 
should  he  written  ns I  =  (8/nn)''-'/cosW,  .  Also,  the  top  line 
of  the  right-hand  column  on  page  194  should  be  written  as 
". . .  error  from  using  Kq.  (2‘2) . . . ." 
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kiMicivcd  .AiiKUsI  :;h.  lOltll 

A  physical  aiplIcK  iipproxlinalldii  i.sdcrivrcl  tor  liKlit  hy  diclni  lrit'  .splicrcs  wil  h  nd  rai  l  i\  i'  liulit  cs  lass  I  linn 

llu'ir  siirrinmdint's,  and  it  is  appikd  (o  air  luddili>s  in  wali-r.  'I'hi'  apprnxiniaiinn  (;ivcs  llii>  cnarsi’  slrnclnri'  in  ihr 
■si  allarini;  when  Ilia  saallarinn  an(>la  '/>  is  near  ilia  ariliaal  saallariii);  an)  la  i/i, .  wliara  0,  ^  8.1°  Inr  liiilililas  in  walar. 
Dillraaliiai  lias  liaaii  nlisarvad  lo  lia  inipiirlanl  in  Ilia  ariliaal  ra);i(ai  liaaansa  nl  an  aliriipl  alian^ia  in  Ilia  anipliinda 
nl'  Ilia  rallaalad  wava  |l‘.  1..  Marslon,  .1,  Opi.  ,Soa.  Am,  119,  I  i’ll.')  I',"  1 1  1 1!IH0||.  Inlarlaranaa  llial  is  dna  In  a  ridraalad 
wavi  prodiiaas  osaillalioiis  in  Ilia  iiilaiisily  willi  an  annular  quasi  pariod  ol  nianniliida  iX/n  I'  ’  rad  inair  0, ,  wliara 
A  is  Ilia  wavalaiinlli  and  a  is  Ilia  radius.  Unlika  dilTraalion  ralalad  osi  illalinns.  Ilia  inlarlaranaa  nsaillalioiis  in- 
araasa  in  maniiilnda  as  Inrward  saallarinn  is  appmaahad.  Opiiaal  liinnalin);  llinainli  linlililas  is  also  discnssad. 


INTRODUCTION 

An  approximntedeHcription  of  the  seiittering  of  light  try  large 
(lielei'trif  spheres  may  he  obtained  from  ray  optics,  provided 
that  the  oliservation  angle  is  not  in  a  region  in  which  wave- 
front  diffraction  is  dominant,'  In  the  ordinary  case,  in  which 
the  refractive  index  of  the  sphere  exceeds  that  of  the  sur¬ 
roundings,  forward  scattering,  the  rainhow,  and  the  glory  are 
the  regions  in  which  diffraction  is  important.'  '  In  the  less 
explored  case  of  a  sphere  whose  refractive  index  is  less  than 
the  surroundings,  a  new  region  appears  in  which  diffraction 
is  important.'  This  region  may  lie  called  the  critira/  seal- 
lirinU  ri’fiion  because  it  is  associated  witli  light  retlected  from 
the  sphere's  surface  when  its  angle  of  incidence  is  close  to  I  lie 
critical  angle  normally  associated  with  total  relied  ion  from 
plane  surfaces.  Diffraction  effects  in  this  region  have  been 
observed  and  described  with  a  physical-optics  approximation. ' 
The  purpose  of  this  paper  is  to  discuss  the  effect  of  the  in¬ 
terference  of  the  reflected  wave  with  certain  rays  that  pene 
trate  the  sphere  and  are  refracted  to  angles  in  the  vicinity  of 
the  criliial  scattering  region.  This  interference  leads  to  a 
modification  of  the  previous  model'  that  becomes  significant 
in  t  he  critical  region  when  an  air  huhhie  in  water  has  a  radius 
<0.4  mm. 

I,et  the  scattering  angle  0  denote  the  angle  of  the  distant 
oliservation  point  nieasured  with  respect  lo  the  direction  of 
[iropagation  of  the  incidert  [ilatie  wave  and  the  sphere’s 
center.  'I’he  critical  scattering  angle  is  0,  =  n  —  211, ,  when 
0,  =  arcsin  (n  ')  is  the  critii  al  angle  of  incidence  and  n  = 
(q,,Ai,)  >  l,n„  and  q,  being,  reipectively,  the  refractive  indices 
of  t  'le  outer  and  the  inner  media.  Air  hulihles  in  water  have 
n  ~  4/11  and  0,  82.82°.  Figure  1  illustrates  several  rays 

that  have  a  deviation  angle  0  of  50°.  Each  ray  is characteri/ed 
by  a  parameter  /;,  where  p  -  1  is  the  number  of  rellections 
from  the  internal  surface  and  p  =  0  has  only  an  external  re- 
llection.  Kays  that  entei'  the  bubble  below  the  center  line  are 


indicated  with  a  prime.  The  angle  of  incidence  measured 
from  the  surface  normal  is  II,,,  and  the  corresponding  refracted 
angle  is  p,,.  Kays  with  p  e*  0  have  II,,  <  II, . 

The  previous  physical-optics  approximation'  described  the 
diffract  ion  of  the  virtual  wave  front  associated  with  t  he  p  = 
0  ray  only.  It  was  shown  that  diffraction  is  important  in  the 
region  I0-0,  |  "v  (X/a )'  -’ rads,  where  n  is  the  sphere  radius 
and  X  is  the  wavelength  in  the  outer  media;  consetiueiitly, 
diffraction  is  important  near0,  even  if  X/a  is  quite  small.  The 
pre.sent  paper  shows  that  the  interference  that  is  due  to  the 
p  =  1  ray  leads  to  a  modulation  of  the  scattered  intensily  as 


0  ?  I 


Kig.  1.  Itiiy  paths  in  tlie  .scatleriiig  plane  with  a  sealiering  angle  0 
--=  .'■,11'’.  The  nnniher  givi's  the  ray  paranieler  p  'The  dashed  ray  is 
forward  seatiering  lhal  is  doe  lo  Innneling  from  the  relleeted  ray  and 
is  fliseussed  in  .See  I. 
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a  lutu'tion  ot  (/).  'I'he  /)  i  2  rays  will  Ire  omitted  since  they  lead 
to  virtual  wave  fronts  much  smaller  in  amplitude  than  the  /» 
=  t)  and  p  =  1  Wave  trouts.  The  p  =  2'  ray  leads  to  a  super- 
|)osed  low-ami)litude  modulation  of  the  scattered  intensity 
with  an  annular  period  ^A(,\/o),  where  a/.\  is  the  lateral 
separation  ot  rays  0  and  2'  after  emergence  from  the  sphere. 
For  an  air  bubble  in  water,  ray  tracini'  gives  \  ~  t).H2  with  </> 
=  <t\  ,  whereas  observations  close  to  0,  (Sec.  V  of  Kef.  4)  gave 
\  =  0.7(j  ±  0.04.  'I’he  angular  period  of  this  modulation  is 
«(X/(i)'  "  when  a  »  X.  ’I’he  emphasis  of  the  present  paper 
is  on  the  coarse  structure  in  the  scattering;  consequently, 
primed  ravs  are  omitted.  Mie-theory  results''  (to  lie  published 
separately)  substantiate  these  approximations  in  the  region 
(0,  —  40®)  <  0  <  0,  for  n  -  4/0. 

1.  DESCRIPTION  OF  THE  VIRTUAL  WAVE 
FRONTS 

Diffraction  and  interference  effects  in  the  scattering  are 
modeled  by  first  using  ray  ootics  to  describe  virtual  wave 
fronts  positioned  within  the  s(  atterer.  'I'hese  wave  fronts  are 
then  allowed  to  diffract  to  th*  far  field.  Fortunately,  it  is  not 
necessary  to  express  the  exact  amplitude  and  phase  of  the 
wave  fronts  in  terms  of  ('artesian  c<s)rdinates.  'I'his  is  because 
aiiproximate  dependences,  dciived  below,  justify  the  use  of 
the  stationary-phase  approxin  a.ion  in  the  diffraction  inte¬ 
grals. 

Figure  2  illustrates  ('artesian  coordinates  used  in  the  de¬ 
scription  of  virtual  wave  fronts  in  the  scattering  plane.  The 
origin  is  chosen  to  be  the  point  where  the  reflected  ray  with 
(Xi  =  0.  contacts  the  sithere.  'I'he  positive  c  axis  is  the  path 
of  the  scattered  critical  ray  predicted  by  ray  optics.  The 
positions  of  the  virtual  wave  fronts  associated  with  rays  0  and 
1  are  denoted  hy  the  functions  and  I’ltu),  where  the  ti 
axis  is  perpendicular  to  the  i'  axis  and  lies  in  the  scattering 
plane.  IMtase  shifts  of  ray  0  that  are  due  to  relleclion  be.vond 
the  critical  angle  are  not  included  in  (n  bid  are  accounted  for 
with  a  separate  phase  factor  in  the  description  of  the  ampli¬ 
tude,  The  reference  phase  is  selected  such  that  CotO)  =  0, 
Kay  tracing  gives  Cn  =  n  at  ti  =  (),  and  the  leailing  term  in  a 
Taylor  series  for  r,i  is' 

Cii  =  (la  ( 1 ) 

where  o  =  -la  cos  d,  )  The  'I’aylor  scries  for  lu  may  he 
obtained  hy  including  effects  of  ret  ract  ion  and  t  he  change  in 
velocity  within  the  sphere;  the  leading  term  is 

I'l  =  '/.(U/'-'.  (2) 
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The  resulting  wave  fronts  are  shown  in  Fig.  2.  The  calculation 
ol  the  interference  of  Waves  0  and  1  uses  the  exac*  |)ath-lenglh 
ditlerence  that  is  not  easily  expressed  as  a  function  of  ii. 
Klfects  ot  wave-front  curvature  orthogonal  to  the  scattering 
plane  will  he  approximated  with  normalization  factors  to  he 
introduced  suhseciuently. 

In  the  description  of  the  scattering,  it  is  only  uer'cssary  to 
(■(insider  two  cases  with  orthogonal  incident  polarizations, 
since  the  general  case  may  he  ohtained  by  a  linear  combina¬ 
tion.'  The  subscript  /  is  assigned  the  value  1  when  the  inci¬ 
dent  electric  vector  is  entirely  (leriiendicular  to  the  scattering 
plane  and  the  valiu*  2  when  entirely  parallel.  Waves  scattered 
to  the  far  tield  associated  with  p  =  tl  and  p  =  1  rays  may  be 
simulated  hy  a  Witve  on  the  u  axis  with  the  following  ampli¬ 
tude: 


V 


Fi|!.  2.  Positions  of  the  virtual  wave  fronts  for  p  =  0  and  p  =  t 
s(  .iltering.  The  d.ished  ray  is  the  incident  critical  ray  that  is  letlected 
to  tiecome  the  c  axis.  A  portion  of  tlie  huhhie's  surtace  is  also 
shown. 


)  =  h,ij  exp|i(-  fee,,  -  wt  )|.  of) 

where  i  -  v  ~1  .  tc  is  the  frequency,  t  is  the  time,  atid  h  = 
2Tr/X.  'I’h;-  expt-  iuit)  factor  will  be  omitted  subsequently. 
Kay  tracing  and  the  physical-optics  approximation  give  the 
following  amplitude  factors'  '*: 


hn.j=-0(do).  u  <  0,  (4) 

h,).;  =  a  >0.  (Tl) 

/ii.;  =  0,  a  <  0,  (It) 

hi,;  =  1 1  —  CytOi)-')  V  x/)  .  a  >  0,  (7) 


where  the  T;  are  Fresnel’s  amplitude  rellection  loefficietits 
for  a  plane  surface," 


_  sintp,,  -  Ol,) 
sintp,,  -f-  0,,) 

tantp,, -I- tt;, ) 

and  /);,  is  the  refraction  angle  predicted  by  Snell's  law, 


IH) 

(9) 


p,,  =  arcsintn  sin  (>„).  ( 10) 

The  .Stokes  symmetry  relations  have  been  used  in  the  deri¬ 
vation  ol  Eq.  (7),  where  1)  is  Van  de  Hulsl’s  divergetice  factor,' 
which  may  he  written" 


__  _ _sin  l>]  Cos  ll\ 

2|1  -  (n  cos  d|/cos  p,)|sin '/' 

(/.isthe  scattering  angle  ofthep  =  1  ray  itredicted  hy  rayoii- 
lics. 


0  =  2(/)|  -  (),).  (1'2| 

and  (v  is  a  positive  constant  normali'/.ation  factor  determined 
suhse(|uentiy. 

I  he  phase  shilts”  t),  are  associated  with  the  (nearly)  total 
rellectioti  (d  ray  0  as  t/i,  exceeds  (I,  (tunneling  through  large 
hubbies  can  be  tieglected;  see  Sec.  4): 

tan((),/2)  =  (sin'-tt,,  -  a  ' -)'  ■'/cos  d,,.  (lit) 

tan((V;/2)  =  a  "  tan((')|/2).  I  PD 

I  he  (\  vatiish  when  do  =  d, ,  hut  they  become  positive  as  d,, 
exceeds  II,  ,  which  represents  an  (idi  aaccaicat  of  the  p  =  0 
virtual  wave  front.  (The  reader  is  cautioned  that  some  au¬ 
thors  write  (^1  atui  (V.  with  an  iticorrect  sign  for  their  choice  of 
the  sign  (d  u;  in  the  time-depetidence  factor.") 

Aiiproximate  dependence  of  Ip,  ,  on  small  values  of  u  is 
ohtained  hy  first  relating  d,,  to  a  fora  ~  Dandd,,  ~  d,  This 


I'M 
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leads  to  the  following  proportionalities:  1  —  ha,;  «  (—«/«)''"- 
for  u  <  0,  ,'i  i  j  (n/o) '  -  for  u  >  0,  and  ' 

(If.) 

where  =  (H/nii)'  -’  cos  d,  and  /!,.  = 

2.  FAR-FIELD  DIFFRACTiON  AND 
INTERFERENCE  OF  THE  SCAITERED  WAVES 

Let  W  denote  the  distance  to  the  observation  point  from  the 
center  of  the  bubble .  When/?  »  kn-,  the  scattering  ampli¬ 
tude  is  proportional  to  the  F'ourier  transform  of  //,);  -f 
which  may  he  written  as  /,  =  /oj  +  where 

/,,,;(sin  tj)  =  ^  expt- ihu  sin  rjldu  (Ifi) 

and  tj  ~  0,.  —  0.  An  approximation  for/a  ,  with  i)  >  0,  which 
was  introduced  in  iv'ef.  4,  is  to  let  /  /a  ,  =  0  for  u  <  0  hecau.se 
( /(I,;  decreases  abruptly  as  u  becomes  negative.  It  is  not  clear 
if  ray  optics  tan  be  used  to  model  the  virtual  wave  in  the  region 
in  which  |(9/ia,;7du  |  is  large;  however,  its  use  leads  to  the  dif¬ 
fraction  of  retlected  energy  away  from  the  critical  region  and 
supports  the  approximation  Uo.jiu  <  0)  -  ■  Most  of  the 
diffracted  p  =  0  wave  in  the  critical  region  comes  from  the 
virtual  wave,  where  u  is  small  and  positive,  and  on  approxi¬ 
mate  result  is  obtained  by  using  Eqs.  ( I )  and  (15)  in  the  entire 
positive  a  domain.  The  integral  may  be  appr(/ximated  by 
using  the  principle  of  stationary  phase'"  with  the  result  that 
fa  j  =  {\a  cos  0, where'*'" 

gii,;  ==  expdyiij)  If(it')  -  F{-  ®)|,  (17) 

w=  [(a/A)cos  (/,  1'^- sin  rj,  (18) 

7o,;  =  cosO,,)'^'  +  ,  (|9) 

4o; 

wheic  the  first  term  of  Eq.  (19)  is  the  phase  shift  of  the  pro- 
metric  ray  reflected  to  0  and  Fiu  )  is  Fresnel's  integral, 

F(w)  =  ^  exp[i'(7r/2)z“|dz,  (20) 

Wiien  the  observation  point  is  in  the  critical  region,  Eq  (2) 
may  be  used  in  the  evaluation  of  f\  j: 


wave.  Consequently,  the  error  from  using  Eq.  (2,'t)  in  the  total 
scattrrinp  should  be  acceptable  even  if  jj  is  small. 

The  scattered  intensity  normalized  to  that  of  the  geometric 
residt  for  a  perfectly  reflecting  sphere  of  the  same  size  is 

Ij  =  'I'Apo.!-^  P\,,\  ''-,  <24) 

where  the  normalization  was  obtained  by  noting  that  Igoj  |  *• 
\'2  as  ic  ••  “>.  The  actual  intensity  is  the  incident  intensity 
multiplied  by”  Ij{a/I{)-/4.  The  value  of  y  in  Plq.  (7)  is  deter¬ 
mined  to  be  2  by  noting  that  only  then  doesgij  become  the 
value  predicted  by  ray  optics."  Because  of  the  interference 
of  the  pp  j,  the  intensity  is  sensitive  to  the  computation  of  (7,)., 
—  7i,;).  Comparisons  with  Mie  theory  have  shown  that  the 
|)hnse  difference  given  by  Eqs.  (19)  and  (28)  is  limited  in  its 
usefulness  because  of  errors  introduced  by  Eqs.  (1 ),  (2),  and 
(15).  The  agreement  with  Mie  theory  is  significantly  im¬ 
proved  when  (7o,;  —  7i,;  )  is  taken  to  be  the  exact  value  pre¬ 
dicted  by  ray  optics''": 

>().;  ~  7i  .;  =  2/m  (cos  (/|  -  cos(/(i  -  n  “'  cos  pi)  -  r),  (//,,). 

(25) 

where  On  is  predicted  bv  ray  optics  to  be 

On  =  (TT  -  0)/2  ==  0,  -I-  (7?/2).  (26) 

To  evaluate  Eqs.  (22)  and  (25)  for  a  given  value  of  0,  it  is 
necessary  to  soK"  the  transcendental  Eqs,  (10)  and  (12)  for 
0\  and  pi.  The  evaluation  of  Eq.  (24)  requires  the  evaluation 
of  Eqs,  (7)-  (14),  (17),  (18),  (22),  (25),  and  (26). 


3.  DISCUSSION 

Figures  8-6  illustrate  the  result  of  three  approximations  for 
the  sc.sttered  intensity  for  an  air  bubble  in  water  normalized 
to  the  ray-optics  value  for  a  perfectly  reflecting  sphere.  In 
each  figure,  n  =  4/8.  The  ray-optics  result  for  a  bubble  is 
ecpial  to  1  ■+ 2.'i  i  j- when  tj  >  0  and  T;  ((/(i)' when  7;  <  0.  It  ex¬ 
hibits  a  cusp  at  the  critical  angle,  which  is  also  evident  in  the 
ray-optics  predictions  plotted  by  Davis"  and  by  Welford.'- 
The  predictions  for  the  diffraction-only  model  developed  in 
Kef.  4  are  given  by  |g,i ,  |  '^/2.  The  present  model  must  even¬ 
tually  break  down  in  the  region  where  <  0  because  of  the 


fij  -■  /i|j  exp|-  i(/;u  sin  17  -I-  (/o<vu'-'/2)||du.  (21) 

The  pha.se  of  the  complex  exponential  is  stationary  when  ti 
=  (7,  where  U  -  a  cos  0,  sin  t).  Since  h\  j  is  a  slowly  varying 
function  of  u,  the  stationary-phase  approximation  may  he 
used  to  evaluate  Eq.  (21)  provided  that  U  is  well  within  the 
domain  of  integration."’  The  result  is  /i ;  =  (Xa  cos  //,  )'/'-' 
p\.i/'l,  where 


Pij  -  25 1, ,(//])  exp 

'qT  +  7t;  +0(X/a)'-’. 

[  \4  1 

(22) 

7t,;  = 

kisin  1})’ 

2a 

(28) 

In  Eq.  (22)  the  stationary-phase  point  that  is  used  in  the 
evaluation  of  h  \  j  corresponds  to  the  (/i  given  by  Eqs.  (ID)  and 
( 1 2).  The  u.se  of  the  stationary -phase  approximation  breaks 
down  as  r;  -  *  0,  since  U  0;  however,  this  is  just  where  lu  j  ► 


t;  degrees 

Fig.  8.  Three  apprnxiiniitions  for  scattered  intensity  is  a  innetion 
of  !;  =  (/(  —  it>,.  with  II  =  4/8,  the  electric  field  perpetidiculnr  to  the 
scattering  plane,  and  ho  =  ID.OOt),  The  da.shed  curve  is  the  result  of 
ray  optics,  which  adds  the  intensit  ies  of /;  =  t)  and /)  =  1  scattering. 
The  dotted  curve  is  the  diffracted  p  =  t)  wave  only  and  is  from  Kq. 
1 17).  The  solid  curve  is  given  hy  Fq.  (24),  which  includes  the  dil 
trad  ion  of  t  he  p  =  0  wave  and  t  he  interference  of  t  he /)  =  1  wave. 
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Kit!.  1.  Saiiu'  iis  Kit!,  il  l)ul  with  parnllcl  itolarizalioii. 


‘  .'U  v< 

,,  I'KjRfF!; 

Kit;,  •>.  Same  as  Kit;.  tail  ilti  parallel  piilarizaliia). 

assumption  thiit  ( ',i ,  Ui  <  (.1  =  0.  lit  the  l'it!ures,  the  moilel 
has  Iteen  extended  into  that  retrion  hy  evaluatint!  Ki),  (24)  witli 
!,',,(>)<())  =  0.  Fresnel’s  integral  /•’(a  l  was  evaluated  with 
a  numerical  ai>proximation.'  ' 

'The  huhhle  sizes  in  these  litjures  are  manirest  hv  consid 
eritiK  ilKimitiittion  hy  it  He  Ne  laser  with  a  wavelentith  in  air 
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ol’tid'J.H  nm  and  a  wavelength  in  water  X  "s  474. (i  nm.  'I'hen 
k(i  =  10  000  gives  n  =  0.75r)  mm.  When  measuring  the  scat¬ 
tering,  refraction  corrections  nitty  he  required  to  relate  the 
sciittering  angle  in  thr  icn/cr  to  that  seen  hy  an  external  in¬ 
strument.* 

'I'he  diffriiction-only  theory*  predicted  oscillations  in  the 
intensity  thtit  are  homologous  with  respect  to  huhhle  size, 
polarization,  and  refractive  index.  They  decrease  mono- 
tonically  in  magnitude  as  r/  increases  and  have  an  angular 
quasi-periiKl  £(\/n  )*'''*  rad.  Figure  II  shows  that,  for  j  =  1  and 
ka  =  10,000,  interference  alters  the  scattering  significantly 
for  r]  2  10°,  After  the  fourth  maximum,  t  he  oscillations  in- 
crrnsr  in  magnitude.  It  is  noteworthy  that  the  quasi-period 
of  the  newly  predicted  oscillation  is  also  roughly  (X/n  rad 
|as  may  he  shown  hy  Kqs.  (1)  and  (2)1.  Figures  4  (5  show  that 
the  effect  of  interference  is  even  larger  for  j  -  2  and  for  ka  ~ 
loot).  The  intensity  averaged  over  a  quasi-period  is  roughly 
the  ray -optics  result.  Consequently,  for  large  spherical 
huhhles  in  white  light,  the  ray-optics  results  are  useful  except 
lehen  le  <  2.5,  where  diffraction  corrections  become  essential 
tie  =  2..5  corresponds  to  tj  =  4.4°  at  ka  =  1(),(KX)  and  rj  =  14.0° 

at/«i  =  loot).) 

The  experimental  data  plotted  in  Ref.  4  were  for  ka  It) 
OHO  and  j  =  2.  A  comparison  of  Fig.  4  with  those  data  indi¬ 
cates  that  it  is  not  surprising  that  the  diffraction-only  model 
correctly  described  the  relative  positions  of  the  first  four 
maxima  hut  that  it  underestimated  the  magnitude  of  the  in¬ 
tensity  oscillations  for  r;  >  H°.  A  new  apparatus  has  l)«en  built 
that  permits  idrservations  in  the  region  t]  >  7°,  which  was 
obscured  in  the  original*  apparattis.  Visual  observations  for 
air  Indrhles  with  a  ^  1  mm  in  monochromatic  light  suggest 
that  the  magnitude  of  the  intensity  oscillations  tends  to  in¬ 
crease  with  increasing  >)  hut.  as  of  the  time  of  this  writing, 
photographs  of  a  (piality  to  permit  quantitative  measurements 
have  not  been  obtained.  These  observations  preceded  and 
motivated  the  present  theoretical  effort. 

The  intensities  given  hy  Kq.  (24)  have  been  compared  with 
Mie  theory*  '  *  with  n  =  4/1!  and  ka  =  25,  100,  1000,  5000,  and 
10,0(K).  The  computations  were  performed  using  Wiscomhe’s 
Mie-scattering  algorithms*’’ with  minor  modifications.  'I'he 
comparision  shows  ’  that,  for  {)  <  rj  <  ;!()°,  the  coarse  structure 
is  approximately  described  by  K,q.  (24).  Frrors  in  ihe  pre¬ 
dicted  liK'ations  of  intensity  oscillations  are  typically  less  than 
one  tenth  of  a  ipiasi-period.  Furthermore,  the  phase  differ¬ 
ence  of  the  scattering  amplitudes  for  the  two  polarizations 
given  hy  Rq.  (2H)  (below)  reproduces  the  coiirse  structure  in 
the  Mie  results.  Rquation  (24)  differs  most  from  the  coarse 
structure  in  the  Mie  results  when  j  =  2.  This  may  he  due  in 
part  to  its  omission  of  the  p  =  2  wave,  which  has  a  smidler 
amplitude  for  j  =  1  than  for  j  =  2.  For  example,  at  >/  =  111)° 
the  intensities  of/)  =  2,  scattering  in  the  units  of  Figs.  11  (1,  an 
given  by  ray  optics'  to  be  0.012  and  0.015,  respectively. 

RquiOion  (24)  does  not  describe  fine  utruelure  in  the  Mie 
intensities  which  is  most  prominent  when  j  =  1.  'I’hisstruc 
ture  hits  a  quasi-period  ~.\(X/fi)rad  10.11°  for  ka  -  1000), 
which  indicates  that  it  could  be  ineluued  in  the  model  with  the 
addition  of/)  =  2'  scattering,  Ray  optics  verifies  that  the 
amplitude  of  2' sciittering  in  the  criticid  region  is  largest  when 
;  =  1.  K.xeept  for  the  line  st  rueture.  thi^  Mie  intensities  de 
crease  f,riidually  as  r;  decreases  through  zero.  I'he  error  in  Rq. 
(24),  with  ))  negiitive.  genendly  decreases  wit h  l;a  as  diflrac) ion 
becomes  more  important,  ’ 
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iMiuation  CJ  t)  iiiiiy  In-  cunvcrlfd  to  the  units coinntonly  usoil 
ill  Mil'  alttoritlims  throiitth  llic  tollowintt  Iransl'orinatioii: 

(27) 

wluTi'  N,  arc  the  scattcrinj;  amplitudes  in  (he  usual  units.'  '  ’  ''' 
The  [liiase  dilTerenee  of  .S'l  and  'S'.i  is  needed  it  the  seatterinn 
for  arhitrary  stales  of  the  ineident  polarization  is  to  he  eal- 
eulated.'"'  The  phase  differenee  is  predicted  to  he'  ‘ 

arntN'.d  -  arn(.S'i)  =  arntA’ii;,.  +  a’i,,')  -  arnU',i,i  +  aM  i), 

(■2K) 

where  arn  denotes  the  complex  arnument  of  the  (|uantit >'  in¬ 
dicated.  The  present  calculation  uses  an  expt—  iuil)  time 
factor,  whereas  the  S,  are  often,'  '  '  hut  not  always,"’ specified 
hy  usiiij'  ex|)(+  iwt ),  so  the  rittht-hand  side  ol  Kq.  (’28)  should 
he  mul)i))lied  hy  —  1  for  some  aiiiilical  ions. 

The  coarse  structure  in  the  scatterintt,  which  arises  because 
of  the  interference  of  the  rellected  and  refracted  littht,  may 
he  a  useful  experimental  tool.  'I’he  annular  iiosition  of  the 
maxima  should  shift  if  there  are  channes  in  the  size  of  a  huhhie 
or,  for  the  case  of  a  drop  with  n,  <  ri„,  if  the  refractive-index 
ratio  channes.  .Shifts  in  the  structure  of  rainbow  scatterinn 
have  heen  useful  for  detectinn  shape  oscillations  in  droiis.'” 
and  a  similar  techniipie  could  he  used  for  detectinn  small 
channes  in  the  scatterinn  from  huhhles. 

Diffraction  associated  with  the  critical  annle  and  structured 
scatterinn  that  is  due  to  interference  should  also  he  present 
in  the  sciilterinn  hy  cylindrical  huhhles  or  dielectric  cylinders 
with  n,  <  n„.  l'’or  a  circular  cylinder  with  a  symmetry  axis 
perpendicular  to  the  incident  beam,  the  critical  scatterinn 
annli’  nti<l  ray  tracinn  in  the  ncrittcrinn  plane  are  the  same  as 
those  described  here  fur  s))heres.  The  divernence  factors  used 
in  the  descritition  of  the  virtual  wave  fronts  differ  from  those 
of  a  sphere.  The  nfotut’tric  optics  of  a  larne,  perfectly  re- 
flectinn  cylinder'"  show  that  (nenU’vtinn  tutmelinn)  lh(i,(u  > 
())|  is  no  lonner  constant  hut  is  '  lsin((/i/'2)l'  The  normalized 
mannitude  and  positions  of  the  intensity  oscillations  should 
(lifler  from  those  of  a  sphere.  Critical-annle  diliraclion  and 
interference  structure  should  tdso  he  present  in  acoustic 
scatterinn  from  lluid  spheres  and  cylinders.  "’ 

4.  ELECTROMAGNETIC  TUNNELING 
THROUGH  SPHERICAL  BUBBLES 

The  present  model  htis  assumed  that  the  reflection  is  total 
when  III,  >  II,  riie  existence  of  evanescent  waves  within  the 
huhhlc  wiil,  however,  frustrate  the  total  reflection,  as  in  the 
(  asc  of  plane  dielectric  sltdts  separated  hy  an  air  nai>.  '  I'he 
cITcct  of  lunnelinn  on  nearly  totiil  reflection  at  curved  di 
elect  ric  interfaces  has  heen  studied  to  determine  losses  from 
curved  optical  wavenuides  ' however,  the  sense  of  the  cur 
\,'ilurc  differs  from  that  in  the  present  case.  The  effect  of 
tunnelinn  throunh  huhhles  can  he  estimated  hy  usinn  the 
|)lane  slab  results.  '  in  which  the  nap  i/  is  c(|ual  to  I  he  huhhie 
diameter  and  ll„  becomes  the  annle  of  incidence  at  t  he  plane 
siirlace.  The  result  of  I  his  approximal  ion  is  that  the  wave 
I  hat  has  I  untie  led  t  hrounh  the  laihhle  will  always  he  trans 
milled  in  the  forward  direction.  The  dashed  rav  in  Kin.  1 
shows  the  rav  transmitted  because  o|  the  Irusiraled  relleclion 
ol  I  he  /  I  =  11  rav.  For  each  /;  =  I)  ray  wit  h  i/,,  ^  //, .  |  here  is  a 
lorward  Iransimtted  ray  with  r  =  II,,.  The  inicnsilv  trans 
mil  lance  l  hrounh  a  |)lane  nap  i',  foil  nil  to  he  ol  order  t/.i/ 1 
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when  0„  =  II,,  and  it  deereases  rapidly  as  11,,  exceeds  II,. 
Tunnelinn  should  have  a  nenlin'I'h’  effect  on  (he  reflection 
from  a  huhhie  unless  ha  is  smtill.  This  conclusion  is  tdso 
sut>ported  hy  the  suicess  ol  the  present  model  in  re))roducinn 
Mie  results"  when  ha  =  ’2.'>. 
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Mil'  si'iitti'i'iii);  iilqiiril  hms  witi'  iisrd  In  ruin|)iili'  .M-iittnri'd  mli-nsitii’s  nnd  plijisc  dil Irri'iu'cs  Inr  iiir  htddili's  in 
Wilt  IT.  Hi'siilts  ill'!’  |)li)t  ti'd  its  it  liinrtinii  III  t  111'  M'lil  Icrin);  nn);]r  I/I  111  I  111'  qcMri'id  I'iinqc  III  tit  I  iMi“  Inr  si/t>  iinriiiiii'li'rs 
All  111' ‘.l.'i.  Mid.  1  .(Mil),  iiiiii  lO.IKHI  iriirri‘s|iiiiidiiiq  111  riidii  II  ^  l.d/iiii  In  ll.M  niiii I.  ,'\s  i.i  ili'iTi'iiscs  lii'lnw  llic  rrilirnl 
si'iiUi'riii);  iiM);li'  id  |  lie  iiili'ii>ii  v  iiirri'iisi"'  iiiid  iiiidi'iqncs  liriindly  siini  i'd  iiM'illiil  inns  I  Inil  iii'r  drsrrilird  liv 

ii  idivsii'iil  ii|il ii'.s  ii|i|iriiMiiii;it inn  di'Vi'lii|M'd  liv  llin  .'IiiiIiiH’k  in  n  M'pni'iili'  piililn  iilimi.  Mir  sriil Irrinq  nlsn  rxliiliils 
rinriy  spnrrd  iiM'illiil inns 


INTRODUCTION 


( )lis(‘l'val  inns  Ilf  lip'll!  sen  1 1 1'rcd  liy  nir  liiilililcs  in  wiilcr  in  I  hr 
yicinil  y  id  I  hr  ci'il  irni  srni  tni'iiiq  anpic'  indicalc  lhal  (lilTr.-K 

I  inn  is  sip'iiifirain ,  even  if  ( he  hnhhie  rad  ins  n  is  as  larqe  as  I 
inm.  The  enarse  s(  riienii'e  in  I  he  eril  ieal  seiMteriiif;  repiiin  has 
heen  deserihed  wilh  a  |)hysieal  n|)lies  apiiroxinial inn  thal 
aeenuiils  tnr  I  he  di  (Trail  inn  id  I  ip  lit  rid'leeliid  frnni  ( he  hnhhie's 
surface  and  the  intertereiiee  nt Certain  refracted  lipht.  ’  In  this 
Letter  we  use  Mie  (henry  ’  '  In  ealeillale  tar  fiedd  seallered 
intensities  and  phase  diflereiiees  in  (he  erilieal  repinii  id 
spherical  hnhhies  in  water  and  (innpare  these  exact  results 
wit  h  tlinse  111  I  he  physical -npl  ies  apprnxiinaliiiii. '  I’revinns 
applieatinns  nt  Mie  t  henry  In  hu hides  einphasi/.ed  I  he  eah  n 
lal  inn  I  it  the  seallerinp  etTieienev  fact  nr''  ''  and  the  radial  inn 
pressure. ''  In  eniilrasl  with  the  seallerinp  hy  spherical 
(Imps, the  exact  anpular  striieture  ir  the  seallerinp  hy 
huhhies  is  nnt  wei!  expinred. 

( 'nnsiderat inn  nt  the  apprnxiniale  seatterinp'  intrnduees 
(he  features  In  he  expected  in  I  he  exact  Mie  re.^ulls  and  ta 
eililales  t  he  ehniee  nf  t  he  ranpe  and  trecpieiicv  id  I  he  imp  Ies 
needed  Inr  a  d  esc  rip!  inn  nt  I  he  anpular  si  riirl  ure.  I''ipure  I 
i  I  lust  rales  seyeral  rays  wit  h  a  seal  teriiip  anple  i/i  nt  .ad".  In  t  his 
( ipure  and  (he  Mie  ennipul  at  inns  I  hat  Inllnw.  (he  ridat  iye  re 
t  raet  i\'e  index  ii,  /ii.,  is  taken  In  he  'M  \,  where  n,,  is  I  he  relrae 
t  ive  index  nt  the  water  and  ii,  is  (he  rc  tract  ive  index  nt  I  he 
hiihhle  ennteiits.  Hays  are  eharaeterized  hy  a  parameter  /i. 
w  here  /  i  —  I  is  t  he  numher  nt  reflect  iniis  wit  hiii  t  he  hulihle  and 

I I  =  n  has  Hilly  an  external  ret  led  inn:  K,,  is  t  he  anple  id  inei 

deuce  at  the /it  h  ray  when  it  lirsi  tniiches  I  he  siirtace  nt  I  he 
hnhhie.  When  sin  II  '  //,  hi,,,  t  he  /)  =  n  ray  is  Inlally  nd'lecled 
(till  eftects  III  tunnel i lip  '  and  siirlace  curyaliire  are  neplectedi 
Tile  cnrrespnndinp  cnndilinii  nn  the  scatlerinp  anple  id  (he 
rel  lected  ray  IS  .,'1  '  w  here  .... .  I  he  cril  ical  scat 

lerinp  anple.  is  related  In  m,/ii,,  I  hnitiph  Lps.  I  1 1  and  ( L’l  nt  Het. 
I.  (I'nrair  hnhhies  in  seawater,  yi.  may  hens  larpe  as  .s  I'M 
The  physical  npl  ics  appinxi mat  inns’  ’  lead  in  l  he  predict  inns 

III  si  met  ure  in  t  he  anpular  scat  I  erinp  wit  h  an  anpular  spread 


III  the  nrder  nt  I  A/ii  I'  rad  |  ~  I  l■|(/,■M  I  '  depl.  where  A  is  I  he 
wayelenplh  in  t/ic  mili'r  infdiiiin  and  i;  =  .Structure 

wilh  a  (|uasiperin(l  nfihat  mapnilude  will  he  releried  Inns 
ciiiirsf  slnirlurf.  in  cniil  rast  In  /iiic  ,s7n/c/ mit  in  I  he  scat  lerinp 
I  hat  has  a  cpiasi  perind  m  ( he  cril  ical  rep  inn  nt  /c.s.s  I  Inm  A/ii 
rad  Inr  :t()()°/A'(/ 1.  The  I ine  si riicl ure  nripinal es  in  part  t rnin 
Iheinlerlerencenfihe/)  =  'J  raxwilh  the  rest  nt  I  he  seat  lerinp. 
'I'n  (lescrihe  the  peneral  efi'ecis  nt  the  critical  anple  nn  the 
seal  lerinp,  Mie  cnnipmat  inns  are  needed  wit  h  a  spread  nt  |  i/i, 
-  i/i|  sniiiewhat  larper  than  ( A/ii )'  '  with  anpular  steps  snine 
what  smaller  than  \'2ii  In  prevent  sainplinp  ermrs. 


COMPUTFK-PKOGRAM  CONSIDKRATIONS 


The  Mie  snlutiniis  were  enmputed  hy  iisiiip  the  l-'nK  I'K.W 
MIKX'P  siilirniiline  cleyelnped  hy  Wiseninhe, chnseii  in 


I'ij;  I  \i:\\  I  III  1 1  is  III  t  III'  ‘'c;il  I  ciini;  pliiiic  \vj|  h  ;i  s(  ;H  t  cnin;  iiliulr 

'»o'  Till'  iiiiinl»rr  ;i(|  lilt  nit  tin'i'K  Ii  IMS  i:iM''- !  Ill'  rav  |i;ii';itiH'l  ur  /  » 
Tin-  apiattMiiLil  um  )',i\('n  in  Ki't.  iin  Indrs  the  iiili'i  h'n'iicc  and 
dittiat  t  D*n  III  w  iiM  s  a  >'''Mial  «'d  with  i  as  s  n  and  I 
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pri'lcrcu'c  In  Ms  I'nsU'i'  Mll-'A  I  null  iiic  licnmist’  nl  |  he  incmnrv 
liniiliil inns  111  llw'  HP  KHII)  i  nnipulcr  sv.^lnin  used  in  llm  vvnrk. 
\\  f  tniiiiil  ll  nci'i'ssary  In  t  nincrl  nil  Wismiinlin’s  inicrnni 
vnrinlilcs  In  HP  (Iniililc  prcrisinn  ( I  I  sinnilii’iinl  (ii>iilst.  iiml 
we  wnrc  I  hen  nl'lc  In  duplii  nlc  liis  '  nlilns"  liir  scnl  Inrinn  liy 
n  drnp  willi  /i,//i,,  =  l.n.  iWiscnml'n  's  sulirniilinrs  iissmiif 
I  I  (litlit  sinuln  pri'cisinn.  I  MnsI  i  nniplnxanlhnu'lii  wnsrr 
wi  ll  I  ('ll  as  rnal  aril  linicl  ii.'  Iim  aisc  nf  I  lie  alisciu'n  nl  cun  ;<lc\ 
dniililc  priaisinn.  'I’hc  ln|;arii tiniii  dnrixalivc  nl  llic  Hi 
latti  Bessel  Innelinn,  eninninnly  deiinled  liy  .A,,,  was  laltn 
laled  liy  dnwnward  reeursinn,  wideii  is  known  In  lie  slalile. 
'I'lie  iirnurani  was  tested  against  \\  iseoinlie's  tallies  tor  /i,/ii.. 
=  l  a  nl  the  eoniplex  sealleriiiK  ainplitud'.  a  and  elTii  i,  nev 
I'aelors  and  anainsl  the  sealtei  ini;  elTieieiu  y  I'aetors  lor  small 
si/e  parameters  when  n,lu,,  ~  iViTi  as;'  veil  m  Uel's.  7  amt  H. 

Pel  iV  denote  the  mini  her  ol'  lerins  iiu  hided  in  the  i-ompn 
I  at  ion  ol  the  Mie  series  and  /V,||,|,  dennle  the  iinniher  ol'aiuiles 
r.ir  which  the  Mie  seatterinn  is  coniputed.  Sin-e  MIKVO 
uses'’''ii  iV  =  kn  +  l.uriiA'ii  l'  ■'  +  2  and  iV.,,,,,  inusi  also  lie  in 
creased  in  proportion  to  kn  lo  resolve  the  line  siriu  c,  exe 
ciil ion  I inie  is  roughly  proportional  In  l/.n )  ’  Inr  a  lixed  spread 
111  angles,  Whereas  the  results  lor  kn  ~  ‘.2.7  re()iiired  a  Tract  ion 
111  a  mimile,  those  Tor /id  =  1  n.DIIU.  wheie  I  he  ant;le  step  si/e 
was  1 1.(111.7°,  reiiuired  'J.a  h  Tor  Htltltl  antdes.  I'I'hese  times  were 
reali/ed  alter  vecinri/al  ion  oT  I  wo  tV,,,,,  Inojis.  wtiich  speedeil 
up  execiilioii  hy  a  Tactor  oT apiiroxiiiiately  1.1 

Hiulier  Values  oT  .V  than  that  ttiven  in  Wisconihe's  eipialion 
above  were  tested  Tor  the  liuhhle  ease  with /.'d  =  HKIand  KMH' 
and  louod  to  yield  identical  results  to  six  decimal  places. 

In  contrast,  the  physical -optics  approximation  program  was 
executed  in  seconds,  regardless  oTsi/e  parameter.  K<|ualions 
used  in  that  pro),'rani  are  listed  al  the  end  oT .Sect ion  II  oTBcT. 


KBSUr.lS 

t  'alculaled  iiilelisilies  /,  as  a  Tuiu  t  ion  nl  I  he  scalteriiin  anyle 
i.'i  are  shown  in  I'd^s.  2  -1.  The  iioriuali/alion  was  chosen  so 
that  /,('/'!  =■  1  represeiils  perTecI  reTlecI ion  according  lo  (jeo- 
met  l  ie  opt  ics, '  I  '  Al  a  distance  /\  »  kn  '  Trom  the  center  oT 
the  sphere,  the  actual  ;  pnlari/ed  inleiisilv  is  the  inciileni 
;  pnlari/ed  inteiisiu  multiplied  hv  / ,(n ///)  '/ i,  where  ;  =  I 
denotes  polari/al  ion  nl  the  electric  vector  per|)eiidicular  lo 
I  he  scat  I eri lit;  plane  and  ;  =  '2  t  he  paralhd  case.  I'he  van  de 
1  lulst  iinrinali/.al  inn '  used  in  Wisconihe's  rout  ines  reipiired 
that  the  sipiarcd  niodulus  oT  the  scalleriin;  ani|)lilude  com 
piited  with  MllsX'l)  he  multiplied  hy  12/knl  '  to  he  expressed 
in  I  he  units  oT  l''it;s.  2  I.  ’ 

.\s  is  shown  in  the  Tinures.  I  he  exact  sol  ut  ion  sui't'esls  that 
a  Tine  ripiile  slriicliire  is  stiperposed  upon  the  coarse  structure 
that  is  described  hy  I  he  phvsical  o|)l  ics  approxiinalion.  The 
ripples  have  a  ipiasi  period  oT  a  lion  I  (l.H'J  (  \/n )  rad  near  i/i, . 
dcci c.isini;  slightly  with  To  avoid  severe  distortion  re 
siill  iiij;  Trom  saiupliiii;  errors,  it  was  necessary  to  use  an  anule 
step  si/i-  less  I  ban  A,'7n  rail.  This  Tine  si  ruci  u re  is  dll"  I o  the 
iuterlcreuce  ol  ravs  not  included  in  the  approximation,  [iri 
iiiariK  I  he  p  =  'J' rav.  ( 'oiiseipient  ly.  I  he  ripples  should  he 
dauipeil  out  Tor  hiihhles  cniilainintt  a  itas  that  ahsorhs  littht. 
I'  iiie  structure  from  air  hulihles  has  lieeii  observed.' 

I'he  plots  in  I'  itt.  .7  show  I  he  |ihase  diTTcreUce  oT  the  seal 
I  ered  waxes 


(b) 


T'.i;  t'lileuliiled  iiiiriiiali/.eil  sciillerinn  iiilensilies  Inr  All  =  lU.IMKI, 
The  electric  vei  Inr  is  |ii.riillel  In  I  he  srnlienni;  pliiae  The  snliil  curve 
IS  Irnm  Mil-  llienri.  'The  dash  curve  is  the  plivsieiil  nplies  a|) 
prnxiiiiiil mil  e.iveii  in  Hel 

Ti  =  ai'n'.'s  d  -  arttl.'s  1 1,  ( 1 1 

win  re  the  .S',  are  Mie  siallerini;  amplitudes  Tor  a  lime  de 
lieiidence  oT  expl  -i  u  l ).  where  w  is  I  he  Tretpielii  v  and  I  I  he 
time,  'I'his  phase  dilTereiu  e  can  he  used  alonu  with  1 1  and  /  - 
to  compute  the  scallered  inleiisity  and  polari/alion  tor  arhi 
Irarv  stales  oT  incident  polari/alion.'  The  physical  optics 
approximation  Tor  h  is  "iven  hv  l'',i|.  I'JH)  oT  Bel.  'd.  As  is  de¬ 
scribed  in  BcT.  'J.  it  was  necessary  to  reverse  the  sit;n  ot  a  Trom 
the  output  oT  I  he  MIPX'I)  program  because  Wiscomhe  assumes 
a  lime  dependence  oTexpl  (  iu'l  I. 


DISCUSSION 

^'or  ■;>  >  •;>, .  the  attreemeni  w  it  h  the  approximal  ion  is  best  at 
low  size  parameters  where  diTTraclion  aeeounis  Tor  the  scat 
lerint;  in  this  rei'ion  in  the  approximation. 

The  coarse  siruci  lire  oscillations  al  anules  i/i  <  i/i.  are  similar 
to  the  predictions  oT  diTTraclion  theorv.  which  itives  a  renecled 
intensity  proportional  to  t  he  SI  plated  modulus  oT  i''resnel  in 
lepi'als  |1'',<|.  Idtl)  ot  Bel.  l|.  The  Tirst  maxiimun  oT  the  re 
Tlei'ted  inteiisil  V  occurs  at ' 

~  -  I  ,'.:i \/n  cos  //,  l'  '  Iradi,  I'd) 

The  Cornu  spiral  behavior  oT  the  T'resucl  ilileprals  [iredicls 
a  decrease  in  I  he  amplil  ude  ol  I  lu  si-  iul  eusil  v  oscillal  ions  with 
deereasinn  o  and  .ilsn  a  decrease  in  their  -piasi  (leriod.  Tin- 
amplitude  decrease  is  most  exideiil  at  laipe  si/e  parameters 
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Kig,  II.  l.ikc  Kig.  2  l»it  witli  (ill  ha  =  UMK)  mid  thi'  tdciiric  vector 
pcrpciidiciiliir  to  nciittcring  plane  and  (li)  ha  «  UK)  and  the  electric 
vector  parallel  to  the  acattering  plane. 

(Fig.  2)  where  the  diliraetitin  pattern  la  lean  apreati  out;  at 
stnall  ai'/.e  parametera  (he  (lilTraction  decrenae  ia  lost  hecaiiae 
of  the  increaaing  amplitude  of  the  intensity  oaeillationa  aa- 
aoeiated  with  the  interference  of  the  p  =  I  ray.  'I’heae  oacil- 
laliona  are  also  of  quasi-period  rtid,'^  diminishing 

with  (/).  In  addition  to  the  j  polarizations  of  Kiga.  2  4,  (he 
physical-optics  approximation  has  been  verified  for  the)  not 
shown  and  for  holh )  at  ha  =  .'itlOO.  Both  the  approximation 
and  the  Mie  theory  predict  that  the  positions  of  the  coarse 
maxima  depend  only  weakly  on  ),  the  coarse  maxima  lie 
slightly  closer  to  0,.  when  j  ~  1,  'I'he  approximation  predicts 
that  the  intensity  variations  of  the  coarse-atruclnre  oscilia 
tions  are  greatest  when  )  =  2,  which  agrees  with  Mie 
theory. 

Mie  theory  shows  that  the  (iiw-ntrui  turv  intensity  varia¬ 
tions  are  larger  for )  =  1  than  for  )  =  2;  t  his  is  prohahly  because 
p  =  2'  scattering  ia  greatest  when  j  =  1  For  :i()°  <  i/i  <  !•(>“, 
'.hi  piisilioUK  of  the  fine-structure  intensity  oscillations  for 
)  =  2  are  typically  shifted  from  those  for  )  =  1  by  one  half  of 
a  fine-structure  quasi-period,  except  when  i/i  is  slightly  less 
than  ijif 

It  is  evident  from  Fig.  .'i  that  the  coarse  slructure  of  the 
phase  difference  for  Mie  scattering  roughly  follows  that  of  (he 
physical-optics  approximation.  'This  has  also  been  verified 
at  ha  =  2.'i.  The  gently  sloping  curve  of  the  approximate  re¬ 
sult  near  </),  may  he  undersUsKl  in  terms  of  the  phase  advances 

of  the  (totally  reflected)  p  -  0  ray;  (>■>  is  approximalely 
(ri„Aii)‘‘'Ai  near  <h,..  (hanparisoii  id'  Figs,  .h  and  H  shows  (hat 
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the  maxima  in  the  coarse  structure  ul  |<i|  are  associatetl  with 
minima  in  the  coarse  struclure  of  the  /,  An  anomaly  is 
present  iti  Fig.  ."dl))  for  i/>  ~  .17'*.  Mit‘ cotnpulalions  o(  /,  in 
dicale  that  l\  is  (piile  small  al  1  his  anomaly,  for /,(i  =  10(1, /| 
hasa  minimtnn  value  of  H  X  Id  '  in  the  unilsof  Fig.  iUh)  al  7' 
=  H(>.70.''i‘’.  The  plot  of  ("i  for />(/  =  2.7  reveals  a  similar  anomaly 
at  i/>  ~  211°. 

rite  range  of  huhhie  sizes  itssociated  wit  h  Figs.  2  .7  will  he 
illustrated  hy  ctmsidering  green  illtnninal itai  with  a  wave 
length  in  air  of  770  nm  and  it  wavelength  in  water  \  ol  I  Li  nm; 
then  Aa  =  10,000  givs  a  =  (172  pm  and /.'o  =  27  gives  a  =  1.(1 
pm. 

It  is  remarkable  (hat  the  physical oplics  approximalion 
describes  (he  general  fealtires  of  the  broad  iittensily  maxinta 
at  i/i  ~  117°  for  ha  =  27.  A  stationary  phase  ai)|)roximatiiai 
was  ttsed  in  llefs.  1  and  2  toeliminate  a  (thase  term  |i/,(.'.  a  I 
in  F,(|.  '24)  of  Kef.  1]  from  the  diffraction  integral.  The 
arguments  given  in  Kef.  I  that  justify  that  approximation 
wtaild  appear  to  l.reak  down  for  ha  -  27.  The  rough  success 
of  the  model  indicates  that  the  maximum  can  be  attributed 
in  part  to  the  first  ililTraction  (naximum  and  in  part  to  the 
interference  of  the  refracted  wav  associated  with  the  p  =  1 
ray;  however,  27  may  be  dost'  .o  Ihi'  lowest  value  of  hci  for 
which  this  model  is  useful. 

'I'he  princi|>al  inirpose  for  develo|)ing  am'  U'sting  the 
physical-o()tics  ap|>roximation-  was  to  see  if  the  coarse 
structure  in  the  scattering  could  be  relaleil  to  the  iliffraclion 
and  interference  of  rellecled  and  refracted  waves.  The  model 
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Fii:.  I,  l.fki’ Fi^.  2  liu(  wall  ),(i  =  2.'|  mid  (lie  l■l('^■(ric  vis  lor  (ul  per 
|ieli(lieul,'ir  mid  ( III  piiriillel 
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■>  I'liiiM'  tlilli'i'ciH'i'  Inr  Ihi'  MMlIcriiii:  iiiu|)lllii(li's  in  the  Iwi. 
|)(iliin/;;liHii  kind's  liii  liil /,ti  -  liinii  iiiid  ihi  All  =  'I'hrsuliil  nirvr 
l~  Iriiiii  Mic  lln'iirv  iiik'  I'.q.  l  1 1,  Tlir  riiiAc  Is  the  plivMial 

.i|itU's  M|i|ir<ixiiiiiiUoii  liv  r.’Mi  III  Hcl,  ■.! 


miiy  liiivc  |irii('liciil  liiisvcvci',  (Acm  IIupiikIi  il  miiils 

(lie  riiic  si I'lK'l lire  in  ihc  scnllci  inn.  I'"!'  nxmnplc.  in  itic 
si'iiltcrinn  liv  il  splu'l'iciil  )i(ilv(lis|i('i'si(in  uf  Itnblilcs,  the  lini' 
striiciiiri' wniild  lie  Inst  iind  Ihc  iiiinsc  structni'c  rcliiini'd  il  ilic 
sprciid  (it  linlililc  si/cs  nr  sviixclcnnllis  was  mil  Inn  liii'nc. 
I'dirl lu'i'inni'i'.  Ihc  nindcl  niiiv  lie  nsclul  in  llm  design  nl  in 
si  riirni'iil  s  I  hill  drpend  (III  I  lu'  sciil  Ici  inn  III  linlil  by  linlililcs. 
blxiiinplcs  incindi'  nil  diA'iccs  Ini  in  nlii  npliciil  iii(‘iisiii'i‘in(‘nl 
III  iniiTnbiiblili'  pnpilliil  inns  ill  sen;  lb)  liisiT  Dnpplcr  iiiic 
ninninlry  nl  billibics  in  liijilids;  iind  n  I  I  be  dnlcclinn  nl  Imlibln 
linpcidi'i’l  inns  in  idnss.  .Inlinsnii  iind  ( 'nuke'-''  liiivn  nsi'd  In) 
In  nbliiin  evidence  nl  ii  Ini  k  nl  bubbles  III  sen  wil  b  i  iidii  nl  li'ss 
llinn  nppi'iixiiiiiilely  dll  ^n\\.  In  nppnreni  disnipeeiiienl  willi 
siiniliir  ncniisllc  nieiisiirenienls.  'I'lieir  inslrimieiil  innde  use 
III  sciillei'iim  wil  II  i/i  --  Dll" ,  wbicb  is  sbiiW'ii  by  l■’i|;s.  d  nnd  I  In 
be  iniicb  sinnller  I'lir  niicinbilbbles  llinn  Inr  n  i;eiiniel  rie  re 
I'leclinn  ri'iiin  ii  perteci  lerleclnr.  Ibilb  Ibe  pliysicnl  nplics 
nppi'iixiniiil  inn  nnd  Ibe  Mie  llienry  siierresi  Ibnl  liidible  de 
led  inii  insi  rninenls  relviiii;  nil  bisinl  if  scnl  lerinn  sbiiidd  nc 
cepi  lit;l:l  willi  less  I  linn  I  bnl  nl'lbe  lirsl  dilTrniliiin  mnxi 
111  mil  i's  predicled  by  l‘',(|,  I D I ,  ( ieniiiel  ric  npl  ics'"  jdves  /  ,  > 


1  Inr  all  i/i  ill,:  C(iiise(|iienl ly,  il  sij;nit'icniil ly  nveresi ininl es 
I  lie  scnl  lered  inleiisil  V  cliise  In  i/i,  wlien /.'ii  <  11)00. 

Keller"  '  '  Ims  nllcinpled  In  use  liirlll  sciilleiini;  with  </>  ~ 
DO"  In  Micnsiire  Ibe  si/e  speciriini  nl  niiernbiibbles  in  wiiter  I'nr 
llie  piirpnse  nl  teslintr  iiindels  nl  cnvil  nl  inn.  Het'ereiices  M 
nnd  1  .'i  rrive  nii|;illiir  scnl  lerinn  pnl  lirns  (  re  purled  In  be  nb 
liiiiied  Irnin  Mie  llienryl  willi  A  ~  O.lidHH  ^nii/ 1  .ddd  nnd  ii  = 

I .  .n,  7..n,  and  III  pm  I  wliicb  enrrespnnd  InA'ii  ~  1  d.  Oil,  100,  nnd 
Id'il.  Cnnipiirisnn  nl  the  resiills  nl  niir  Mie  priigriini  nnd 
Mi.idel  slinw  lliiit  Ixeller's  pntlerns  err  sip'iiiriennl ly;  Ibey  dn 
mil  include  line  slriicliire  nnd  lliey  incnrrectly  describe  llie 
ciinrse  si riicl nre. 

We  nre  r;rnlel'iil  In  W.  .1.  Wiscnnilie  liir  priividin^  Ibe  inililnl 
cnnipiiler  pmurnm  trniii  wliicli  llie  prn);riini  used  here  was 
derived  iind  In  H.  Medwin  liir  briiininr;  1  be  wnrk  nl  .liilinsnii 
and  ( '(Mike  In  inir  nl  lent  inn.  'I'liis  resenreb  was  siippnrled  in 
pari  by  I  be  Wasbiiii;lnn  Stale  I  biis  erxily  Keseiircb  nnd  Arls 
( 'niiimil  lee  nnd  by  1  be  I  llTice  nl  Nnx'iil  Kesenreb.  I'  I ,.  Mar 
slnn  is  an  Allred  I'.  Slnnii  Kesenreb  |•'ellnw. 

‘  Am  Imr  In  wbiini  cnn  espniidence  sliniild  be  addressed. 
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SUPPLEMENT  TO  PAPER  NO.  4 


This  supplement  compares  unpublished  Mie  computations  (solid  curve) 
with  the  physical -optics  model  from  Paper  3  (dashed  curve).  The  refrac¬ 
tive  index  ratio  is  n./n  -  0.75.  The  notation  is  as  in  Papers  3  and  4 

1  0 

except  in  Figures  SI  and  S2  we  give  the  total  normalized  scattering 

=  (I^  +  I^)  /2  for  the  case  of  unpolarized  incident  light.  For  Figures 
S3  -  S5,  scattering  for  the  same  ka  but  for  the  orthogonal  polarization 
is  given  in  Paper  No.  4;  the  associated  figure  from  that  paper  is  noted 
in  parenthesis. 

FIGURE  CAPTIONS  FOR  SUPPLEMENT 


Fig. 

SI 

Log 

(base 

10)  of 

I  for  ka 
u 

=  25. 

Fig. 

S2 

Log 

(base 

10)  of 

I^  for  ka 

=  100. 

Fig. 

S3 

for  ka 

=  100  ( 

see  F ig. 

3b). 

Fig. 

S4 

‘e 

for  ka 

=  1000 

(see  Fig. 

3a ) . 

Fig. 

S5 

>1 

for  ka 

=  10000  (see  Fig 

.  Za).  See  also 

Fig. 

,  4  of  Paper 

1  for  (p 

from  75“ 

’  to  85°. 

Fig. 

S6 

Pha 

se  difference 

for  scattering  amplitudes 

(as 

in  Fig.  5a, 5b) 

for  ka  =  25.  Note  that  the  model  does  not  give  the  fine 
structure  present  in  the  Mie  result. 
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One  of  the  causes  of  light  scattering  in  glass  is  the  pres¬ 
ence  of  small  gas  bubbles.  In  critical  applications  such  as 
optical  fibers  this  scattering  is  not  negligible;  a  microbubble 
trapped  at  the  core-cladding  interface  may  scatter  a  signifi¬ 
cant  fraction  of  the  light  in  the  fiber.'  Here  we  use  Mie  theory 
and  physical  optics  to  describe  the  main  features  of  bubble 
scattering  in  glass.  We  anticipate  that  these  results  will  be 
useful  for  the  detection  of  bubbles  in  quality  control  appli¬ 
cations. 

We  are  unaware  of  prior  publication  of  the  Mie  intensity 
for  0  180°  scattering  from  bubbles  in  any  dielectric.  Struc¬ 
ture  in  the  Mie  scattering  from  bubbles  can  be  related  to  dif¬ 
fraction  and  interference  near  the  critical  scattering  angle  and 
to  the  Brewster  null  in  the  reflectivity  of  parallel  polarized 
light.  For  scattering  from  a  sphere  with  a  refractive  index 
which  exceeds  that  of  the  surroundings,  the  critical  scattering 
angle  is  not  present;  also,  the  effect  of  the  Brewster  angle  is 
less  pronounced. 

We  approximate  the  bubble  shape  as  a  perfect  sphere  and 


assume  the  glass  has  a  refractive  index  of  1.46,  as  does  the  glass 
u.sed  in  the  cladding  and  outermost  core  layer  of  a  typical 
optical  filler  preform.  (In  the  resultant  filler  a  bubble  trapped 
at  the  core  -cladding  interface  is  broken  up  into  a  sequence  of 
much  smaller  bubbles  elongated  in  the  direction  the  fiber  is 
drawn. ' )  We  take  the  bubble  gas  to  be  air,  so  that  the  relative 
refractive  index  n,/n„  cy  0.68, 'i.  A  bubble  is  described  by  its 
size  parameter /eci,  where =  2ir/X,  A  is  the  wavelength  of  the 
light  in  the  glass,  and  a  is  the  bubble  radius.  Here  we  con¬ 
sider  bubbles  with  size  parameters  2.'i,  and  100.  For  A  = 
0.49  pm  (which  corresponds  to  692.8  nm  in  air  for  visual  in¬ 
spection)  ka  =  100  gives  a  =  6.9  pm;  ka  =  5  gives  a  =  0.9,'i 
pm.  For  A  =  0..'i6  pm  (0.82  pm  in  air,  a  typical  operating 
condition  of  fibers)  ka  =  1(W)  gives  a  =  8.9  pm. 

If  the  bubble  is  illuminated  by  a  collimated  light  beam,  the 
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Fig.  1.  l.ognrithm  (bsse  lU)  of  the  normalized  scattered  inlensity 
predicted  hy  Mie  theory  (solid  curve)  for  ka  -  lUO  und  the  ele<  tric 
field  perpendicidar  to  the  scattering  plane.  The  dashed  curve  is  the 
physical  optica  approximation  given  hy  K(|.  (24)  of  Ref.  :i  which  is 
useful  when  scattering  angle,/)  is  close  to  the  critical  si-uttering  angle 
,/,,  94°.  The  approximation  fails  to  iescrihe  forward  region  (,/,  < 

20°)  and  the  backward  regi  ui  (■/,  ~  180°). 


Fig.  2.  hike  Fig.  1  hut  with  electric  field  parallel  to  the  scattering 
plane  (_/  =  2  scattering). 
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Fig.  :i.  Scattering  for /,'<1  =  2.'').  Mie  theory  for ;  =  1  (dotted  curve 
and  j  =  2  (solid  curve,.  Dasheil  curve  is  from  physical  optic 
approximation. 
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aiiKiilar  (iistrihut  ion  of  the  scatterwl  intensity  is  described  by 
Mie  theory. At  a  distance  R  »  ka-  from  the  biil)))le  center 
and  a  scattering  angle  (/>,  the  ^-polarized  scattered  intensity 
is  given  by  the  incident  y-polarized  intensity  multiplied  by 
/jti/iltu/ZO-V't,  where;  =  1  denotes  polarization  of  the  electric 
vector  perpendicular  to  the  scattering  platie,;  =  2  the  parallel 
case.  'The  normalization  of  the  functions  Ij(it))  is  chosen  so 
that  =  1  represents  perfect  reflection  according  to  geo¬ 
metric  optics.''  ''  Using  a  program  described  in  Ref.  !>  incor¬ 
porating  a  modified  version  of  the  MIKVO  subroutine  devel¬ 
oped  by  Wiscomhe,"  we  obtain  the  functions  shown  in 
Figs.  1  -4. 

'I'he  coarse  structure  in  the  vicinity  of  the  crit  ical  scattering 
angle  0,-  may  he  understcKKl  with  reference  to  a  [jhysical  optics 
approximation  developed  by  Marston."''^  Most  of  the  scat¬ 
tered  intensity  in  this  region  is  due  to  three  types  of  rays  which 
we  denote  by  the  parameters  p  =  0,  1 , 2,  where  p  is  the  numlrer 
of  chords  wit  hin  the  bubble  and  the  p  =  I)  ray  has  only  an  ex¬ 
ternal  reflection.''  When  the  IcK'al  angle  of  incidence  0  satisfies 
sind  >  sinW,.  =  a,  Ai„  the  p  -  1)  ray  is  nearly"  totally  reflected; 
the  abrupt  boundary  of  total  reflection  at  (I,,  results  in  a  far- 
field  diffraction  pattern  of  this  ray.  Rays  with  0  <  0,.  are 
mostly  transmitted  inside  the  bubble.  Rays  with  p  -  1  are 
refracted  twice  and  interfere  with  p  =  I)  rays  in  the  region  0 
<  0, ,  where  0,.  =  rr  -  2(1,.  is  the  scattering  angle  of  rays  inci¬ 
dent  at  W,..  Marston’s  physical  optics  approximation  con¬ 
siders  the  diffraction  and  interference  of  p  =  0  and  1  rays. 
Small  quasi-periodic  fine-structure  variations  ride  on  the  large 
quasi -periodic  coarse-structure  variations  and  are  due  mainly 
to  the  interference  of  certain  neglected  p  =  2  rays  {p  =  2'  rays 
in  the  notation  of  Refs.  II  and  T)).  If  the  average  intensity  over 
an  angular  range  containing  several  fine-structure  variations 
is  desired,  the  approximation  may  he  used  to  predict  /;  in  the 
range  20°  S  0  S  0,  =  911.54°.  'Phis  model  breaks  down  for 
small  bubbles,  as  can  he  seen  in  Fig.  4,  due  to  its  use  of  the 
Fresnel  reflect  ion  and  transmission  coefficients  of  planar  in¬ 
terfaces  and  the  stationary  phase  approximation  of  diffraction 
integrals. 

'Phe  physical  optics  a|)proximation  predicts  a  quasi-period 
S  rad  for  the  coarse  structure,  and  a  fine-structure 

quasi-period  of  0.95(X/o)  at  0, ,  decreasing  '  with  0.  'Phese 
features  are  clearly  evident  for  ha  =  25  and  100,  as  is  the 
predicted  broad  decline  in  intensity  for  0  >  0,.  - 
arcsinj  1 .2{X/a  cost/,  l'^"!  (identified  in  Ref,  5  as  the  approxi¬ 
mate  location  of  the  last  coarse-structure  maximum). 

'Phe  broad  minimum  in  1}  near  1 10°  corresponds  to  light 
scattered  from  the  Mrewster  angle:  parallel  polarized  light 
incident  on  a  plane  glass  air  interface  is  totally  transmitted 
when  ((  =  114.41°.  Phus,  when  ka  is  large  enough  to  use  the 
Fresnel  coefficients,  we  expect  a  negligible  contribution  to  l  i 
at  0  =  IH0°  -  2(114.4°)  =  1 1  1.2°  from  the  externally  rellected 
ip  =  0)  ray.  Phis  0  can  also  he  written  as  2  arctan(n„/n, ). 

For  ka  =  5  ( Fig.  4),  the  fine-structure  quasi-period  is  ap¬ 
proximately  (‘(pial  to  the  coarse-structure  (piasi-period.  'Phe 
minimum  in  I-.',  at  102°  and  the  slight  dij)  around  45°  are  the 
Inst  remnants  ol  this  structure.  Ask  a  is  reduced  below  5,  we 
find  that  the  scattering  pattern  approached  that  of  a  dipole 
radiator,  as  predicted  by  Rayleigh  scattering  theory:  /|(0) 
becomes  a  flat  line,  and  /■.|(0)  is  symmetrical  about  a  minimum 
at  90°. 

We  have  compared  these  plots  with  ol  hers  made  for  n,/n,, 
=  0.75  (air  bubbles  in  water)  I'or  the  same  values  of  ka  and 
found  them,  as  expected,  ipiite  similar.  'Phe  most  noticeable 
difference  is  the  slightly  larger  (pmsi-periods  for  bubbles  in 
glass.  In  both  cases,  the  best  angles  to  look  for  scattering  will 
he  those  angles  where"  0  'S  0,  —  (X/o)''^".  Corrections  are 
sometimes  needed'  when  computing  the  effective  scattering 
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Fig.  4,  lake  Fig.  11  hot  with  ka  =  5. 


angle  observed  by  a  detector  due  to  refraction  of  the  scattered 
light  as  it  leaves  (he  glass  and  enters  the  detector."  'Phe  coarse 
and  fine  structures  present  when  ka  S  25  may  he  used  to 
obtain  independent  estimates  of  a,  'Phat  diffraction  is  im- 
|)ortant  at  0,  can  he  seen  by  comparing  the  Mie  results  to  the 
geometric  iiredictiiai,'  including  only  the  p  =  0  ray,  that  ) 
=  I  and  with  the  physical  optics  prediction  that  /;(0,  )  =  0.25; 
the  Mie  results  lie  much  closer  to  the  latter  than  the 
former. 

'Phe  resuh  that  scattering  is  best  observed  when  0  S  0,.  - 
(X/d)''''--’  rad  can  also  he  used  in  certain  cases  for  scattering  by 
nonspherical  bubbles  if  their  aiiproximate  size  is  known.  For 
example,  if  the  bubble  is  spheroidal  and  either  the  major  or 
the  minor  axis  lies  in  the  scattering  plane,  the  critical  scat¬ 
tering  angle  predicted  by  ray  optics  is  the  same  as  that  pre¬ 
dicted  for  a  sphere.  'Phe  diffraction  related  shift  in  the  first 
maximum  is  roughly  (X/d)'^"  rad,  where  a  becomes  an  average 
(which  tlepends  on  the  orientation  of  the  spheroid)  of  the 
semimajor  and  semiminor  axes.  'Phe  details  of  the  coarse  and 
fine  structures  will  dei)end  on  the  bubble’s  shape. 

'Phe  approach  to  (he  scattering  described  here  is  to  allow 
the  light  to  enter  and  leave  the  glass  via  plane  surfaces  and  to 
neglect  surface  reflections."  .Scattering  near  90°  has  been 
used  to  detect  bubbles  in  glass."’  For  glasses  with  n„/n,  =  1.5, 
0,  =  90.4°  and  the  Mrewster  scattering  angle  is  1 1 2.0°.  Rre- 
vious  apirlications  of  Mie  theory  to  bubbles  are  reviewed  in 
Refs.  5  and  1 1.  Approximations  for  scattering  by  inhomo- 
geneities  within  fiber  waveguides  have  recently  been  pub¬ 
lished.'-’''  hut  these  do  not  describe  the  angular  structure  for 
scatterers  with  ka  >  I . 

'Phis  Work  was  supported  by  I  he  Office  of  Naval  Research. 
I’.  1,.  Marston  is  an  Alfred  1’.  .Sloan  Research  Fellow. 
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K.  There  may  he  an  exeeplioii  lo  this  miUleline  lor  hirne  spherieal 
hiihhles  la'eause  haeksealterint;  ean  he  enhanced  hy  the  axial 
l'nensin){  nl';)  =  II  (slnfv  rays.  Ohaervaliiins  nlThis  enhancement 
fur  hiihhles  in  liquids  are  deserilied  in  I’.  I..  Marston  and  I).  S. 
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index-matehinn  liipiid  and  illuminalinn  and  ohserving  it  via 
windows. 
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Paper  No.  6 

Description  of  computer  codes  for  Mie  and  model  computations 
(D.  L.  Kingsbury  and  P.  L.  Marston)  [This  is  a  slighlty 
modified  version  of  an  Appendix  in  D.  L.  Kingsbury, 
"Light  Scattering  Near  the  Critical  Angle  in  Air  Bubbles 
in  Water  and  Glass,'  M.  S.  Thesis,  Washington  State  Uni¬ 
versity  (1981).] 
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Append  L  x 
COMPUTER  CODES 


Three  programs  were  vised  t:o  generate  oiir  plots.  MIESC, 
incorporating  a  modified  version  of  Wi.scoinbe's  MIEVO,  produced 
the  Mie  and  CRIT  the  approximation  results.  The  versions  given 
here  store  results  in  disk  data  files  where  they  can  be  accessed 
by  the  plotting  routin  MTEPL  for  di.splay  on  the  HP  2847A  monitor 
or  for  paper  plots  on  the  HP  9872A  plotter.  All  programs  are 
written  in  Hewlett-Packard  RT-IVB  Fortran.  CKIT  tvaluate.s  Llie 
a  p  p  r  ox  ima  L  1  on  cicvolopetl  by  Marston  .iml  Klignsbury  (1981). 

MIESC 

MIESC  asks  for  tlie  followina  inputs; 

1.  REF  INDX  is  the  complex  relative  refractive  index  Uj /n^ . 
The  imaginary  part  nmst  be  non-positive. 

2.  BETA  is  the  size  pnrauio.'ter  ka  “  27Tn/A,  where  A.  is  the 
wavelength  in  the  outer  medivim. 

3.  AST  is  the  smallest  scattering  angle  for  which  scattering 

angle  for  v/hich  scattering  functions  are  computed,  DPHI 
the  angle  ;;tep  size,  and  NUHANG  -  is  the  total  num¬ 

ber  of  angles. 

If  NUMANG  is  even,  results  are  calculated  for  the  angles 

<Pj  AST  +  i(DPHI),  i  =  0  to  NUMANG/2 ,  and  their  complementary 
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angles  (180*-  )  .  (If  NUMANG  =  1,  90*  must  be  the  only  angle. 

If  NUMANG  =  0,  only  the  extinction  efficiency  factor  is 

computed . ) 

For  example,  if  iij  =  1.33  -  i(.03),  n  =  1.0,  and  scat¬ 
tering  functions  are  desired  for  <p  -  40*  to  90*  with  increments 
of  .5*,  one  should  type  in 

/. 7508, -.0282, 40. , .5,101 

MIESC  then  asks  for  the  name  and  cartridge  number  of  an  existing 
datafile  in  which  to  store  tlie  output.  What  data  is  stored 
depends  on  what  the  output  buffer  RBUE'  lias  been  equivalenced  to. 
In  the  version  included  hero,  each  reco’d  v/ill  consist  of  three 
real  numbers  giving  4'/  I<('^)/  ?iid  lj(<i>),  and  the  records  are 
entered  in  order  of  increasing  except  that  the  record  for 

(180*-  ((>)  immediately  fc.lows  that  for' 

If  the  phase  difference  DIFA  is  not  desired,  the  code 
that  computes  PHI,  PH2 ,  and  DIF.^  is  skipped. 

To  bring  subroutine  MIEVO  down  to  more  manageable  size, 
the  original  comments  were  largely  omitted.  The  commented  ver¬ 
sion  is  given  in  Wiscombe  (1979).  Most  complex  variables  have 
been  converted  to  B-dimcnsional  .irrays  or  to  separate  variables 
prefixed  with  "R"  or  "I."  For  example,  NUH  br'came  NUM(l)  anu 
MUM(2),  denoting  the  real  and  imaginary  parts,  respectively,  and 
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ZET  became  RZET  and  IZET.  (This  was  done  because  of  the  absence 
of  double  precision  Fortran  in  RT-IVB  Fortran.)  Most  complex 
arrays  were  similary  converted  to  separate  arrays  for  the  real 
and  imaginary  components. 

Subroutine  names  prefixed  witli  "DV"  refer  to  functions 
in  the  HP  Vector  Instruction  Set,  which  considerably  speeded  up 
the  execution  of  two  former  DO- loops. 

To  keep  the  loaded  program  within  a  32K  size  limit,  the 
original  MIEVO’s  complex  BIGA  array  was  discarded.  (BIGA,  which 
stored  the  A„ ,  contained  as  many  elements  as  the  number  of 
terms  used  in  the  Mie  series.  Thus  for  ka  =  10  000,  a  total  of 
2(10  088)  =  20  176  double  precision  must  be  stored,  requiring 
161  408  bytes.)  Required  A„  values  are  instead  stored  on  and 
retrieved  from  a  system  disk  using  EXEC  calls. 


The  use  of  the  Hewlett-Packard  minicomputer  instead  of  a 
regular  "pay"  computer  system  was  absolutely  necessary.  Extra¬ 
polation  from  Table  8  of  Wiscombe  (1980)  indicates  a  MIEVO  exe¬ 
cution  time  of  over  100  seconds  on  a  CRAY-1  for  ka  •-=  10  000  and 
NUMANG  =  8000.  The  several  size  parameters,  repeated  testing, 
and  dozens  of  plots  required  made  our  choice  of  computer  systems 
inevitable. 
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CRIT 

Required  CRIT  inputs  are  tlie  relative  refractive  index 
RI  (now  assumed  real);  the  size  parameter  BETA;  PTS  and  STRT, 
which  determine  the  angles  for  scattering  computations;  and  LAST 
ANGLE,  the  largest  PH  (<p)  .  N  -  no /i''i  -  1/RI  is  used  to  calcu¬ 
late  the  critical  incident  angle  0,.  ; 

THC  =  arcsin(l/N). 

(Arctangents  are  used  in  the  program  due  to  the  unavailability 
of  the  arcsin  fuixction.  )  The  critical  scatt  ering  angle  <f>  is 
then 

PHC  =  Tr  -  2  ( THC ) . 

For  a  given  incident  angle  TH-THC,  the  angle  PH  at  which 
the  refracted  ray  is  scattered  is  given  by  an  equation  derived 
using  Snell's  Law  and  elementary  geometry 

IH  -  2(arcsinl sin(TH)/Nl  -  TH). 

For  a  desired  smallest  PH,  one  chooses  STRT  such  that  TH 
=  STRT(THC)  substituted  in  the  above  equation  gives  this  PH, 

As  CRIT  executes,  TH  j s  stepped  up  in  increments  of 
E  =  (THC  -  THj„j^  ) /PTS .  As  TH  approaches  THC  the  corresponding 


step  sizes  for  PH  avow  steadily.  Consequently,  at  selected 
angles  TH2 ,  TH3 ,  and  TH4 ,  the  step  size  in  T’H  is  reduced. 


The  function  G1  i5J  the  geometrical  divoigence  factor 


5irt(TH)  cos^TH) 

G I  =  — j f - - - _ — 

21 1  -[f^c«>»{TH)/cos(R)J|si.,(rPH) 

in  the  expression  for  the  intensity  of  refracted  radiation 
(Davis,  1955;  van  de  Hulst,  1957),  where  R  --  arcsin(Nsin(TH)  )  . 

IP  and  IS  are  the  iiitensi ties  of  the  parallel  and  per¬ 
pendicular  refracted  radiation,  respectively,  computed  using 
Fresnel's  intensity  transmission  coefficients  (1  -  FP )  and 
(1  -  FS)  for  plane  surfaces.  The  coefficients  are  squared 
because  this  radiation  is  refracted  twice  --  at  entrance  and  at 
exit . 


The  variable  W  is  the  upper  limit  of  integration  on 
tlie  Frest\el  integral  \ised  in  the  reflected  radiation  calculation 
(Marston,  1979).  Subroutine  FRES  evaluates  this  integral  using 
a  numerical  approximation  (Abramowitz  and  Stegun,  1955). 

The  phase  differences  B1  and  D2  between  the  reflected 
and  refracted  ray.s  for  perpendicultir  (Bl)  and  parallel  ( B2 ) 
polarization  have  three  compjonont;'. :  the  phase  shift  the 
reflected  ray  undercjoes  at  the  bubble  sui  f ace  (DEhl  and  DEL2 )  , 
the  phase  patli  difference  ETAC,  and  a  Tt/4  phase  difference  aris¬ 
ing  from  the  curvature  of  the  scattered  wavefront.  (For  the 
reflected  ray  this  phase  shift  is  already  included  in  the  Fres¬ 
nel  integral . ) 
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CRIT  is  now  ready  to  compute  the  scattering  fvnictions  at 
angle  PH: 


St  -  *■  (]/^)(FC  +  iF3  )  exp  (  -  iBl ) 

Sj  =  VTp  +  (l/f2)(FC  +  iFS)exp(-iB2) 

where  FC  and  FS  arc  the  output  of  subroutine  FRES, 


MIEPL 

Output  was  plotted  by  protjram  ftlTEPL  using  the  HP  1000 
Graphics  package.  MIEPL  is  not  included  in  the  program  list¬ 
ing  which  follows. 
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I  f  f:  ^  ^  f  ♦  f  ♦'  f  f  ♦»  f  f  C'  I  H  T  r-  i  i  'Mj  -.'  r**-' '  /  *  ».  *  1.  *•  ' 

BLOCr  C'hTh 
C OH P L t : .'  ■-  0 )-i I'  K  ,  S  I  ,  ? 2 

D 0 ij i;: I,,  i-  P R E  i:  1 : 1  i::i r (  e E  r  H ,  ;< rtu ,  p  i 

lOMHOH,  '  IHOO  r/pl  ■  2  >,  '.-.•BwE  I-  ,  01  •>  0'1  2  ),  S2'-  -  i  2  •  .  EETm  OHU'.' 5  t  2  '  .  M2i.  LIT  , 

♦  UUrU'lHi:,  ,  ilHI"  1  ■  HH  ,  OE:-:  r  ,  LUOUT  /  tCIIHt''.  lL!..t<  ,  I  ’^E'.  T  .  ME:0,  1  UP  1  U  ,  1  ..Ip  fcl. 

END 

PPOGPhH  Mti'LC- 
r  Nfiaul.HP  ['IMEtr'-lONL. 

COHPLE:'  'iBilOt  ,  01,  '02 

DOUiOLE  PPE'.  J  "■  I  Oil  PE  f  >■! ,  2rU.i ,  F'  I  ,  C'  r  ^  R  1  ,  i..  UN  .  H';-T 


REwL  N.':;t  LIT  -  I  1  ,,  I  IN  rEN 

D  MIEN'-..  ION  ir-iIBi  *  M  •  I'pNt-'' T  1  .  NkNE-;  7  >  _ 

rOMMuN,-'TNri!iT,>  I  2  '  •  I'Ni.'K  . 1  •;  f-l  2  >,  S2'.  0  1  2  '  PE  Tn  ,  OMl.i'.  0  1  2  ' ,  tUip  , 

nunmng,nhpi  ,nn.>'E:;t,eliolit  icnmd,  ilen  ioclt  h&u,  luPn-MuP^u 

EOU  1  VhLEMi'  E  ■  P  cNF'.  I  i  ■  NNGLE  ' ,  ■:  PPl.lFv  2  ■  ,  I  1  •  ,  ■  Pl'Llp  i  .  I  2  ' 
dntm  nr, ‘I  i.i.'.  rrv'pE/?/,  i o-c o. ' ,  iuPTt),'';i.',  il.'O.- 
r.  il=din’ Pci.iF '*•2 

ELI  Inal. UOLN'  .1 C-  .> 

WRnE<U.lIN,2000  >  _  __  _  , 

i'OOO  FORMmT':  "I  !•' I'lrOLi'I  .  ■■■  ri'JPCOUT  .  LOMPUTE  PHh'-E'^'  '  (  ■s  UR  N  f  ,  . 

R.Ev4t>  '  LUlN.2'.iOI  >  IP,  fO  r,  tPH 
2001  FORNhT  'I  3h  1  ' 

IF  <  1  R  .  EO  .  1  MY  ■'  LLIOUT  = 

IF  (IPJlE  lH,  '  LLIuNT  =  ENIN 
N2CUT  =  0, 

PI  =  3  ,  1 -1  1  0 
CON  =  P  I .  1  ;;'C'3 
OCON  =  'ONOl.'  ilON  ' 

•in', 2  rnJ'HMT^'^niTEr''’cNPE':  PEP  IND'-!,  BETN,DBFT,  m'TT,  NFMI,  NUNNNG"  > 

PFmD'  1  "IN,  t-  ■  r-t.  t  '>.RI.  .:"-.BETN.DDET,M'OT,DrNT  .NLirVtHi, 

i;  M  LN';.  I'OR  NUHE.p  IChI.  INTEGRPTIO.Ei  O'^ER  'iuElD>  hNuLE 

SM'JT  -  CIH'  :.MnE'  n  -L  ON  ) 

t'oTER  =  '  f'C  M  1  t-'.-CON  ■  f  *2 
IFL  I  M  ,  ML  ,  1  (  ,  i.OTO  2 

NRITE'.  ENIN  .  EOnP  '• 

COi'iT  FOPI'ImO.  "PUT  :i-  LH-:.  •.  CR  NO, >'2  CH..  - 

Pf-,,,'  Ln.l  IN  ,  .■-'04  •  NnME'  I  '.tiHlTE'  2  >  NmNE'  3  >  •  I P 

;  ( M 1  4  0 1  ■  r!  r '  n  2  ■  i  . .  ■' 

I  i-iEL.  'j; 'CM'  U'l'f'-.  ,  lEP  NhNE  ,  lOPTN  ,  I'iC  ,  In"  ' 

It  ■  I i-t,':'  , L I’  11  '  O'O  ro 
pr.iiF'  I  ■  -  E4--rn 
R&i.iE'  I'  I 

PEUr :  i  p  i  _  ' 

CnLL  'l|.  I  IT'  if"  I.  .  U'u:'  ,  ODUr,  IL  > 

HI  ll'll  ll.iH  -  '  N'  2,4  !i'.,  '■?  I  2  •■  •  ( 

r  ir  lii.ir.l’nL  .  'nil'.’il  EE  2 1  PpL- ,  PIUlT  PN  I' I  ,  :■  B  •  f  <•  C. 


N'  M’lw  =  O'  l ' '  .  (i.il  lnn:  ,  ,  '. 

Cffi  14  Ii'  =  1  .  'lUO’i:  O’t 


HMPt 

KH  =  HMP  I 
DO  4  I  a),  HI  I 

4  KMU'v  I  *  =  C*CO-''  '  n’;.r  +  '  I~1  itOF'Ht  >''CoH 


t  call  mievO'  n:  - 

IF  '' NUMPHG  ,  EO  .  0  GOTO  t -i 
TINT  ==  Cl. 

DO  111  I  =t,tUI 
RK  =  Cl  . 

hMGLE  =  h'oT  +  '  1-1  'i-r'PHI 

PHI,  PMO  t^Pr  THF  I'tli'-LE'i  IH  cCOPECO  ;  PHL'II'  1 0  MFFEPEHCE 

IF  I  IPH  .  ME  I  il'.'  '  GOTi.i  ; 

PHI  --57  ,  Ow‘v}  •••■  M  I  HM'  hH-IhG'  •.-■1'  I  ■  '••'PEhL'  'il'.  I  ‘  )  C 
PHO --57  L'.I'.' ?  '■>  mTHI-K  Htl'U-iC.’  I  !  '-'PEwL'  1  '< 

IF',  REmL'  5  t '  I  'I.LT.I.I.  I  FH  1  •■sPH  n  I  PCI . 

IF'', PHI  .QT.I'OO,  >  PM1 -Phi  -  .•'CCi  . 

IF',  REhL<  Si'  t  '  ■> .  L,  r  .  Cl .  >  PHE=PHO  I- 1  00  . 

IF''  PH2  ,  GT  ,  I  ,  '/  PH.:  -PH2~76u  . 

PHC'IF  =  PH2-PH1 

IF'.' PHDIF  ,  GT  ,  1  Oi.'i .  >  PHC,'ir=PHMF-7i:M;i. 

IF':  PHDIF  .  LT  ,  -  1  SO  .  '  (MIL’IF^PMC'IF  +  ^it'O  . 

7  H  =  ■:  PEAL'' 51'  I  <  ■  o* '*'2  C'  A  IHaG>  5  1  ’  » 

11  =  '"  n  +1:1  i' ' 

12  =  1' PEaL':  i  l'f  I  )  •  >'♦  mINaG>  ■■2';  1  >  >  '  ♦  '  2 
12  3  ':I2<'C(2.' 

IHTEM  =»  ALOCT-:  .  5*..:  I  I  12  >  ■' 

DEGPOL  =*  '  12-11  I2+I1  ■' 

IF';  AHOC".' hMi'iL.E,  2  ,  ,  EO  .  0  .  H.iP'i  TE',  Lt'OMT  ,  lOOl  j  h1''2i_E,  I  1  ,  I  .i 

IF':  I'i  r ,  El) .!  H  C  >  chll  I'Ih I tF'.:  I t)C e ,  iepp , Pfei.iF  ic  ■ 

DO  Ci'iNPLEMEMThP'C  AHGLE 

IF  COHC  ,  mHGLE  nni''- 1 01 IiE'OIP.ED,  iHE  FIMmL  ['O-LOOP  in 
SUBPO'JT  IMc  HIEVO  iIhMMOT  L;  IC L.CiDE  5'  MNP  1  -  J  '  FM  5 
ANGLE  =  100  -mMGLF 

IF';  MONhMG  .  Ec  .  I  1  GO  i'C  I  2 

lie  --  '.PEhU'  OI ilMPI- .1  '  I  i>  v2+'.  H  '  MAG'.  0  1 '  MMF- 1  •  T  '''♦<'2 

11  =  'll  '-0  2  ' 

I2C  =  '.PEhL  MMFl-'.  .  :  .)  *  t  2  +  '  H  IMAG'.  ;'2'.  Nl'ir  1  -  1  ■  '♦♦'2 

12  =  ''.IE'-*);)':'* 

IHTEM  =  mLOGT'  ,  G".  I  1  +12  '  ■ 

Ti-iFOE  LIME::.  Pt-PF'OPN  ''•Ol  It'  wNi.LE  IMT^ijF'-'  101!  mT  PH1  =  PI/'4 
S  2INT  -  ''CO'C'  PI>'GCijii  I  <i'*2.i''  I  I  +  I  1 1-  '  *■  '  ''..‘'jM  ■'  '♦  t.2*''  le,'-)  I  2i-  > 

TINT  =  TIMT*,^ 'EIMTi-fvTEI- 
RK  =  Pi;  t-DPHI 

IF  ’  I  M',  hMGI.E  I  C'C  i-'M  '  t'Cn-C',  I'.l'  lOCON  •  .  CiE  , 1-'';  r  •  GOCi.i  S 
I  p.' I 'l'■l''’t' '  '0  O’ L  E  ,  2'  '.El''  I'l  ,  '  MP  I  T'C  ■  LUOI.* ’’  lO+'l  E  ,  i  I  ,  T  2 

Id  cnnTIMME 

RLIMT  -  mLOGT'  TIHT  ' 
angle  =  EETh 
n  =  PLINT 
12  -  FLINT 

IF':  I  '  .E'-'.  1  H  '  CMl.i_  I'lPlTF''.  IC-'2£  .  lEFf  ,C'iilF  i  L. 


M-LC'  i-':f  'OOlic-  hngle  It. "LGrt.2r ic.it  ' 


IF':  tv", 

El-' 

.  1  H 

'  1 

I'lPHE'  1. 

.l.O'.iMT  1 

I's  I) 

Ic  .I’CT  '  1 

'3 

■' j  L 1 

•FTEP 

T  .  1' 

f.r 

12  PET.. 

PC. 

1  H  + 

i-L 

ChLL  E. 

■xc 

' ,  5  ,  - 

1  1 

IF  '  ^ 

r  A 

l.F.  .  ' 

14  i.'mNTIMi 

-!C 

ijHir*  I-'. 


2  0  07 

L 

1  0  0  0 


1  C  ij  1 
1  0  02 


C 

C 


C: 

C 

c 

c 


C' 


r. 


IF':.  I'iT  ,  EO  ,  I  MS'  )  ChLL  CLOSE';  IDCS  > 

WRITE'-  LUIN.,  2  007  .  lEOR 
FOkMR  I'  "  IERR="  ,  IS  ) 

USE  this  FORMuT  for  MORE  DEThILEC  OUTi-|.!T 
FORNhT',  1  H)  SEw  ■  HIE  SI2E  PuPhMGTEP  =  .  FS  .  2  ,  1  5: : PEFPpC  T  I  VE  , 

'  I  hue:;  = ' .  d  1 4  .  'S  .  e  1 2 .  s//  ■  hhcle  '  <  1 1  ' ,  s  -sup-  i  '  ^  2  ix ,  -  s-sue-2 ' , 
t  IS'S  IHTEHSITV  ■  .2X,  DEGPOL  '  ,8,  ,  'II  '  12  '  .' 

|■ORr!HT';  F  7 , 2  ,  SE  i  4  ,  ,  FS  .  4 , 2E  I  4 . 6  .  Eo  ,  3  .' 

FORMhT^  F7 , 2,  P  1  2 , 3  .> 

EN[i 


SUBPOUTIHE  MIEVO'..  I  C 
CCiNPLE!':  SBAlK,  St..  S2 

double  precis  I  oh  BETH,.'--:  mu,  pi 
REAL  H2i:.UT 

CONMOH/IttOUT.-PI';  2  •  .-.G-hCK,  SI';  St  2  S2'  St  2  > tSiTH.,  HNU';  512  >,  N2CUT  , 

HUHHHG.tniPt  ..:iii,C;I-:;:T,LUOUT,  ICHUD,  ILEH,  ISI-.:.  T  ,,  MPO,  IOP1  0,  10720 
LOGICAL  HO  I, HAG,  HOhMGS 
COMPLEX  AHM 1 . BHM 1 

C>  I MEHS I  OH  PSP';  25S  :> .,  I  SP';  258  i ,  PSM.;  258  ,  I  SH>  SSC.  .■  .  PSPS';:  256  ,'> ,  I  SP:S';  25' 

♦  ■■?SMS';  256  > ,  I  SMS'  256  < ,  P  I H'.  256  .) .  P  1  HM  1  C  25':- .'  .  1  l  iP'  256  .' ,  I AUH';  2:56 

DOUBLE  PPECISIOM  PBIGm.  I  B  I  Gf^i ,  P  I  H  ,  P  1  NM  1  ..HPIDU,  COEi- F  ,  VTEM';  256  , 
PSP.  I:;.P,PSM,  ISM, PSPS,  ISPS.PSMS.  ISM.S.Ht'H  .?  ■  .BH'.  2  :).PANP,  lAHP, 

"■  RtOIP,  IBHP,  RhMPM.  IhHPM  ,  PE'HPM  ,  IBHPt!,.  R2E’  ,  I.CET  ,  P2ETH,  I2ETH, 
R2THP  1  ,  I  ZTtiP  1  .,  TMP  ,  TmUH  ,  PEZ  ,  REZ  I  HO  ,  FH  ,  PH  ,  i-HP  1  ,.  TWCiHP  1  ,  RNP  1  , 

+.  z  I HV';  2  ;> .  FF''  2  ,  ap'.  2  > .  deH';  2  .> ,  m.iM';  2  ,  eps  :  ,  epsz  ,  temi'  2  :> ,  tma':  2  :>  - 

t  HTH'..  2  ),  DTD'.:  2  .  TI'.  2  )  .  PIOPI  V  .  FS.IH,  P’.-:INP  t  ,.  CHIU.CHIHPt  ,8UPF':  .32  > 

Eour'/HLEHCE  ';si''  1  .-'..PiH'.  1  s;.'-,  i  :>.TAt.tH>  i  < 

+  S  t  '■  257  .> ,  P  I HM 1  •:  1  :> , ''  82':  257  ) TMP'  1  '  ' 

DATA  EPSt.''1  .D-2/,  EPS2.-' 1  .  D--0/ .  MAKI T.  -  1  0  0  On. ' 


IF't  BETA  .  LT  ,  0  .  .OR.  P  I I  :< .  L.E  .  0  .  .OR.  R.  I 2  ,  GT  .  0  ,  STOP  10  00 
HOAHGS  =  MUMAi-G  .FO,  I.' 

NOIMh',  =  RI'  2  ,H2C.Ur 


ChLCULmTE  HUMf'EP  OF  TERMS  It)  MIE  SEp-ICS  '  A  LEAST  UPPER  BOUHL 
US.IHG  EMFIPIChL  FOPMULhE  FITTED  FiT-  SIZE  FhRAMETEPS  UP  TO 
2  0, 0  0  u 


IF'.  BCTA  ,  L  T  .  JlOO  .  )  fiT  =  D'3T^-'^'l  05.tPFTHf".1 
I  pi;  BE  I  A  .  GE’ . -*2  U  0  .  .!  HI  =  BE  Th  ^  4  l-ET'h  )  >  r  t  , 
HTPI  =  HT^1 


COMPUTE  BIGh 


IF 

1C 

••  1 

MpRh;;Y  DIGA  HhS  hLPEADV  BEEH  GEM 

IF  'I 

c  .  c 

T  .  1  " 

■  GO  TO  Z’l'O 

ITRhP 

- 

1  4 

p  c: 

r-  ■- 

•":C.  ijH  .'S-TEil  ri'IC  FOR  aPRA  /  6! 

C  n  L  L 

CC’C 

C  ■  4  , 

ITAmI  I  ;  1  f'C  ,  IL'Isl  ,  L.-.'E'  T  :• 

I  CHAD 

-T 

I  D  I  C 

c 

ILFH 

r2  1 

2  C. 

T  r.r.  1 

i  ■  ' 

[ 

.IP' 

;  I  T’ .  ; . ;  '  t 

I  or  2 

•-  1 

'.  r  c  ■ 

ChLL 

ch:[' 

'•/<  i  t , 

1  [iO  M  ,  .  M  ,  tCM  ■ 
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2INV':  1  ,>  =  TErK  t  i.'EETn 
ZIHV'.  I'  )  =  TEM'::  2  >/BETh 

PREPARE  FOR  t.iOi,iM-RECiJRPEt<EG - 

COMPUTE  IHITIhL  HiCH-ORC'ER  BIGh-U  USING  LENTZ  METHOO 

15  fF<  I  i  =  NTPt  -  ZINV<  1  .) 

FF>.'2  ')  =  HTPl  vZiNV’:  2  ;> 

MM  -1 
KK  =  2*NT+3 

AK<.  1  '  =  <:  NM  *MK  .i:*^ZINV';  t  .> 

AK'<  2  )  =  <:  MM  rKK  i  +  ZI  NV''  2  > 

C'EN'-  1  )  =  AK';  1  ,) 

BEIN'  2  .1  =  AK<  2  '> 

CALL  C>i:  [>‘vN'  Cl ,  1  C> !)  ..  0  ,  .  Ff  ,  TEM  > 

Nur;<  t  >  =  den<  1  ■)  +  TEM'  I  :> 

NUM<  2  >  =  C'EN'C  2  >  +  TENi.  2> 

KOUHT  =  I 


2  0  KviUNT  =  KOUNT+I 

IF','  h;OUNT  .  GT  ,  MA;-'  1  T  )  GO  TO  40 

CALL  [iCC'VC  -  1  ,  NUr-K  1  ) ,  NUI’.'-'  2  '> ,  hK  ,  TEM  > 

CALL  DCUV'-: -1  .  BEN';  1  ) ,  DEN-.  2  > ,  .4; .  TMw  > 

I  F  ■:  TEMf  1  )  ,  GT  .  EF S  )  ,  riNB  .  TMA';  1  )  ,  GT  .  EPS  1  )  GO  TO  '1 0 

ILL-COMBIT lONED  CASE- -STR I  BE  TWO  TERMS  INSTEAD  OF  ONE 
mm  =  “MM 
KK  =  KK'  +  2 

AKiv  !:)  =  '■  MI'i  +  KK  i^ZINV’,  I  :> 

A  K I  2  :>  =  M  M  If  K  K  )  f  Z  I  N  V  ■:  2  > 

NTN’,  I  )  =  AK''  1  >fNUM''  I  '  -  hK<  2 '>i"NUM':  2  '  +  t  D  0 

NTW''2y  -  AK''  I  NUM-:  2  ■>  +  hK >  2  >|<'MUM'.  I  ) 

DTCu  1  =  aF  1  iifDENi-  1  )  -  HK(2>fDEtK2>  +  100 

DTD<  2 ''  =  A1  ,  1  ifC'EN'-.  Z  '  AK»  2  FfOEN';  1  > 

CALL  DC0V''U,NTN'  t  >  ,  HTtr  2  > ,  DTD  ,  TEN '> 

TMAC  1  )  =  TEM'  1  "i-F';  t  '  -  TEM'.  2  "FF'IZ  .) 

FF<2.»  =  lEM':  1  <  f '' f- ■  2  )  +  TEM'"  2  ('•:  1  > 

FF<  1  :>  -  TMh'.'  1  1 
MM  =  -MM 
KK  =  KKf2 

AK'i:,  1  ,)  =  '!  MM  H  h  ;  '  JIN'v''.  t  .> 
hK';  2  =  ':  MMvkF  >  fZINV'  2  ' 

CALL  DC[>V':  0  NMM'  I  Ili.lM'  2  ),NTN.TEN.) 

NUM''  1  ■'  =  Hi  1  >  f  TEM'  1  ) 

NUIT'  2  ■  =  i-F  '  2  '  f  TEM'  2  ' 

ChLL  OCOV'.  u  ,  L  E(H..  !  >,,  DL'N'  2  '  ,  DTD,  TEM  > 

C'l  N'  !  .'  =  At.’:  I  ’  +  TEM'  t  / 

I'lN'  2  '  =  hK  2  '  +  fEM’  2  ' 

K  GLINT  =  lOUNT+l 
GO  ■"0  2  0 

."0  CALI  OlC'V'.'  n  NUM',  I  '  .  Nl  IM’  2  ’  .  DEN  TT  t 

TEN'.  1.'  =  TT’  1  rfFF’,  t  ..  -  TT'  2  >*.FK’.2) 

FF'  2  '"  =  TT'.  1  'f  FF'  2  '  +  TT'  2  ’  "FF-.,  1  ) 

FF'.;  I  ’  .=  TFi;'  •  ! 

■  "f  c  r  F  C'  ■  r  on','l'F CFf iC  t 

IF',  CihCi".  r  I  '  ;  t  Dll  ’  .  L  r  .  EF22  .hND.  DhC:  T  I- •  2  '  l.T  EPS2  i  go  to  50 

MM  =  -MM 
KK  =  I  K-»2 

^1-.  ■:  I  j  ^  ',  MM  r  J.  I  ;  t.“  1  tl’,.  I  ) 

MK'  2  '  -  '  Mi  l  '  !  t  :>  t  2  IN  .  ■  2  ' 
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C  uLL  [>CCiVi,  M  1  [>0  ,  I'l .  ,  HUM  ,  TEM  > 

HUM'  1  >  =  h(  •  1  '  TEM'  1  ■ 

HUM'',  i'  '  =  Ht  i  E  TEM'!  i  !' 

CULL  [>C  ['7'!  0 ;  1  C>  0 , 0  ,  ,  DEM  ,  TEH  > 

DEN',  1  !'  =  nT',;  1  >  TEH-,',  1  ,> 

DEN'-!  2  >  =  ut  I  2  !'  TEMi:!  2  ■' 

CO  TO  20 

40  I'lPITE'  LUOUT  ,  :?.Oi:i)  >  HT,  BET  A ,  P 1  ,  UK  ,  NUM  ,  DEN .  TT  ,  FF 

JOt  FORMhTv/,'/  coin  INUED  fr..-!CT10i!  for  M-'f.UD-HT  FnILECi  ro  CONVERCl 
+  ■  HT=  ,1..-,''  '  =  '  ,  E.-:u  .  0/  REFR  1  H['E::=  '  ,  2E2  0 , 3/  aK=  '  ,  2E2 0 , 3/ 

*  '  HUM-  ,  vE2fi  ,  '  DE  H-  '  ,  2E  2  0 , 0/ '  TT-',2E2  0  3/  FF=  '  2E2  0 , 8  !) 

'STOP  10  02 

5  0  RBIQh  =  FF'  1  .> 

IB  I  GO  =  FFi 2  : 

BUFF':  El!)  =  FBI  Cm 
BUFF'  22,)  =  IBiGn 
MB  -  31 

DOMNWwRD  RECURREHCE  FOR  BIGa-H 

DO  "0  H  =  NT, 2,-1 
TEM':  I  '  =  N»2INV''  1  ’  +  RBIGh 
I  EM*,  2  1  =  M  t  ZIHV'  2,'  IBIGA 
CALL  C'C  C)V';  0  .  1  T'  0  ,  0  .  .,  TET1 ,  TMA  ) 
rbiGa  -  H+ziHW',:  i  :>  -  TMA',:  i  :> 
lE-iIG'-t  =  H^Z1H’■•'^2>  -  TMA'  E',' 

MB  =  MB -2 
BUFF':  MB  >  =  PFIGn 
BUFF'!  MB +t  >  =  IBJGh 
IF  '"  MB  .  HE  ,  1  ,'  GO  TO  7  0 

ChLI.  ECEC'  2  ICHi.iC,',BUFF.  JLEH,  lOF’l  ,  I0P2  ' 

IF  '  H  E ,  2  ,'  'TC'  T 0  ,S 0 
ME  - 


;ir|  '  =  I'l'Pl-l 

I0P2  ^  ('SECT -2 
?0  COHTIlii.'C, 

C  t-'LL  ELFC'  2,  I Hi'T' .  c'L'FF  .  ILL'  IC'Fl  ,  I  OF 2  > 
SO  MB"  ~  T:C: 

I  OP  10  •=  I  DPI 
lOP'  '  =  i;oP.i 


;ui.'  CulM 'MUi: 

Ml:  MLU 

lOPI  -  liiPIO 
I  OP 2  =  [OP2''i 

Call  FHEC'  l  ,  IC  tU.IC' ,  luff  ,  ILEH,  IGPIO,  IOP20  ' 

R  I  OP  IV  -  ID'!. 'PI'  I  ■ 

ItUI  IwL  T  '7  '.'"hUT  I  T  I  L;,  UC-LD  FOP  EF  F  i 'I  I  I'.f  I T  r  hLL  ULaT  I  ijH  OF 
HOIii-fiCmL  C'O  7FT0  IEHTC  IH  Mlt  •OERrE'.-- 
FH  = 

PH  -  ID" 

MN  =  t 

IHlTlA'LlCu  P  ! '.  v  I  T  I -EE'S'OEL  FUHLTIijtl  t  Th  i’l'ip  iJPMaPD  PECUPPEHC 
FoIH  =  I'  ;L'  lLTh  ' 

CHIM  =  !!"''■■'  '  FE  T.  > 


P':.  INPl  P'S IH, 'BP, Ti-i  -  PhIM 

CHINP!  ^  SHIH.-'e!,- Tm  t-  P'}.  in 

P  Z  C  T  f  i  =  P  I  { t 

izcrn  =  cHifj 

PZTMP!  P'l.  IHPI 
IZTMPI  =  i.HIMP) 

imitiplitp  ppppiup-:-  cuPrPTC I irnTS  .  m --ji'p-m - 
POP  USE  It)  M'i.  ,l!METr',  PmOOP  SETtlE''^ 
hNU!  =  '  O  .  0  1.1.  I-  t 
l.'tjl) )  =  I  0  ,  0  ,  u  ,  It  .1 

UtXT  =  -V  . 

It)  I  T  1  III,  T  .■"£  Htir,l..ll,i;i-  r-')i  Pin  Htlf'  'SUH'.S  ri..|-  .+ 


b  ULL 

0  . 

C>  0  . 

n  ^ 

! 

MM  ) 

CALL 

f.l''U'lMV-, 

Ij  , 

Dll , 

u  , 

I  fP,  ! 

.MM  ^ 

CALL 

[''/I'l 

0  , 

,  t'  O  , 

0  / 

b  'bfl,  ) 

,  fUi  J 

CALL 

F'Vfl.CVb 

!.! 

.  0 , 

0 , 

; '.  fi .  1 

,  M'M  > 

CAlL 

Lr.tpIflV' 

1^1 

cm:i 

n  , 

.  l'"  :.|'  S  , 

1  MM 

C  »l,.L 

C'vrtoV' 

fi 

,  D'l 

0  , 

■  * 

1  MM 

CALL 

DVIIO'.,'.: 

M 

DO 

‘1 

P- .bV; 

1  MM 

CALL 

DVUS'Vt 

.  D  0  . 

n  ^ 

,  i  :."i 

?  rif» 

CALL 

DVfl.'SU' 

.  [  0  . 

,  , 

,  I  UP! 

.  •  M 

Call 

C'VUOV' 

1 

,  DO  , 

•J  J 

,  f  ■'  1 1  ,  1 

,  MM 

!0R  I  TE'' LUOUT  ..SO'-..'?  ■  i)r,l'UUtlT.  1 1.  ,  t  II  IMhI  IS 

formp F'  s:, , It;. ,e  I z  .  7  ' 


DO  5  00  li  =  I  ,UT 
PE'.IGh  =  B'.'Pt-.  Ht-.; 
isis^j  =  rurP'  Pb' 1  ' 

I.IM'IPI.IIE  T''E  './-iBIOU':  li  'I-'IEF  M  hL  COP.IT  i  F  I  pH '■ 
run  r  ft 
i  l,,il.'!l  if-  I  -  !■  tl  ■  tlf'  t 

;  I !p  !  =  '.  [ •  I •  If  1 

..  L  rr'  =  f  1 1  t  :  I  ; 

lirU.'U  =  I  ['F- 1 


ChLuI.U.hTE  the  MIE  '-.EPiE'-.  b '0  E)- P  i  C  I  i..tl  T'S  lITT 


IF',  UOIMpG  •  G'li 

TO  ~  1 ' 0 

GEUFPhL  fM 

SE 

CALL  r.:C  C' I.I  . 

DSh,  IBlGb  f 

I  Tf,-'.  • 

TEt'!'  1  :■  ^  TCf'.' 

1  I  I-  '.'^U.'CETi-) 

iriA'',  1  :>  =  TKn 

t  ■♦!-r:T;ii'i  - 

TEU'  Z  '  1  I~ 

TI-iA'  Z  '  -  TE'-'' 

1  )  1  1  CTHP  1  *■ 

TE'!'  Z  '“PZ 

H,  !■:  1  '  =  '  TE'.'i 

.  1  ■  tt- : I !!,-  • 

.  f  D'h 

'  TMhv I  ■ 

‘  1  r  i 

(  V  J  ' 

Mi)',  z  1  =  •  ’  TEI'I 

',  2  •♦P  Itf'  l  '* 

rj'.i-t  !  )  - 

Ti'1'r-i'  !  • 

*  y  y,  »  ’•  V-;  • 

^  t  Z  ' 

FEM'  1  '  =  PI'  1 

'l‘■•r■GM  -  PI 

■  .0  ■■  b  F  p'  i  biH 

TEI'D  Z  1  =  PI-  1 

•  b  !  ...  l  t’j  M  ^  'P  I 

'  ■  '  i  i'-'  I  C.H 

Ti'U-i'  i  ■  -  i'Ci-' 

;  '  t  - 1 T  t :  I- ;  - 

1  L '  i  E 

-  TT 

I  ■  1  i  z  ■  t"  1  ^ 

TC'-  ;  .P" 

Lli'  1 


j  r  I  ■  1-1  ■  I  i  '  ■  I  III'.  I  '  ' 


I  I  F’  ; 


BE  Th 


1,  B-'JUB-N- I  .) 

.  ;  -  PT  i-)LL  hUGI.  F  ;■ 


tiEEr>ED 


Lt-H  hHI'  L1TTlE--B 

1 : 
t! 

P'MiiP  I  i  *  iMu'  Z  '  ' 
!  ir  •  -P'S  I  I)  '■♦'Tl'ln'  Z  '  ' 


I  b  b . 


h 

.1  F"'- 1  Ur  I  Z 


7110  COMTiHilC; 

H  0  -  M  B  Cl  R  P  T  I  0  M  C  H  t 

TEN'.  1  >  ■  PI  OP  I  /♦PRIi,H<rtl.'BETM  (♦■PZTHF  1  -P.- ETM 

TEH''  C'  '  =  PlOPlVt  JUHfFH. 'CL  !«  >•♦  irTHP  1  M  2LTM 
TMA';  1  1  =  '.RIOPIVtP|iIGP»-*rn/&ETA  )i*PSIHP1 -PilH 
ChLl  Dcr.''  ■;  0 .  rii.!'  t  .  1 1 ,  ,  teh  .  hh  ■ 

TEM'  1  I  =  <■  I  .  1  ')*RZTHP1 -PZETN 

TEH''i  )  =  V  PI'..  1  -♦<P.Eiir.M<-Ftl,'eETP  >*l:'THPI -lE'ETtl 

TMO'I  1  I  ~  ■  PI'  1  f.pRIGMfPH/BETi-i  )*PSIMP1 -PC  IN 

C^ILL  [lODVi.  U  ,  Tl'Ur  1  > ,  0  ,  ,  TEN  ,  BN  < 

.750  0E;-;T  OP. .  r  •-Tl'.ii.ulP' 1  '-.NGL'  nN',  1  >  .'♦CNC.L'.  C-iN'  t  ' 

IF'I  NCmNQS  O'TiTO  4S0 

PUT  NIE  CFPIEC  COEFFICIENTS  IN  FORM  NEEDED  FOR  COMPUTING  S+, 

PmNP  =  nOFFP  <"■  hN'  1  ,'t-p.N'  1  )) 

InNP  =  C I'lErr  f- HiP  ' t-E'i'!'.  .2  >  > 

RDNP  =  Ci;iEFr+'  Ml!'  1  '-  DN';  t  '  • 

I9NP  =  EuCFF  hN<  2  •■•DN'.:  2 
RPHPM  NM'ipuNP 
IMNPN  =  I'll'l  f  IhMP 
ReNPM  =  hl'Ui'T'IIP 
IBNPII  =  Mi'WICNP 

rtC'D  UP  CUi'i',,  (.iNll,  E  UPWARD  PECUPCING  ON'.JLmP  FUNCTIONS.  LITTLE  PI 
NND  LITTlE  ThU 

CALL  DVMPV'  ,;HU,  !  .  P  I  N  1  .7TEN,  1  .  NN  ) 

CALL  DVSNB'  VTEM.  I  .PINNl  ,  1  ,TNP,  1  , NN > 

CALL  DVSnT’'  FN,TflP,  1  .T.-'.UN,  1  NN  ' 

CALL  DV'jl.iLi'.  TmUN  I.PIt.Nt.!  rMUN,l,NN> 
i.'wLL  C'VMOP  I  IN,  1  .PItiNl  ,  1  ,HN  ' 

Chli,  [ivr  I V' m  I  I'N  ,  tnp  i  .  vten,  i  .  pin  i.nN' 

CmLL  I'v.ipi"  r  INNl  ,  1  ,  ThUN,  1  .  /IE,I,  I  ,NN 
LhU.  DVPI'.'  ,I"Nr  , I'Lii ,  t  P  P  1  ,  I  CP  ,  I  ,  NN  ' 

CALL.  DVPI"'  ImUP  , '.'TEM,  !  .  1  : 1' ,  t  ,  ICP,  I  .Nil  ' 

CALL  DVPIV'.  N'  11 , 7  TEI't ,  I  ,F  iNC  ,  \  ..  PCMS,  1  .  NN  ' 

CALL  C'VPr,''  !i;.Nri!,'/TLi'l,  '  ,  ICIIC  ,  1  ,  IONS,  t  ,NN  ' 

CALL  L'VCUP'  MIUNt  I  ,  rnlJN,  I  ,VTEN,  1  ,  NN  ' 

CmLL  OVL  17'  PPMf- ,  V  I  l.l'l ,  I  .R-.l'l,  I  .P'CM,  1  ,NN  > 

i:m1  I  nvr-!V'  !i;:UP  vifN.i  i  ;.ii.i  I  nn  .• 

CmLL  DVI'I','  fC-U!-:!  7  I  Li-!  .  t  .P'CP'’. ,  i  kCP'-.  ,  '  ,  NH  ' 
i.'mCL  U'.'f-  I  7’'  i'.UPN  ,  if  N  I.ILP',!  JCPo.t.NN' 

UPI'hTE  RtLFVMNI  OI.Im'NT  1  T  1 E C  PmR  NE;-;!  PmCC  IUFOUGH  LOOP 
45 U  I'tM  "  -  MN 

FN  ---  FHfM 
RN  Rtif-I 

hNI'I  I  -  I  I'lC'l  '.N'.l  '  Ml!'  I  '  '  ,  iNiM,'  hN>  2  '  '  > 

ENMI  =  I  i;'-  '.  .111,1  '  UN-  I  -  I  'I'NC.L'  PNv  2  >  '  ' 

imIi'iImI.  :  i'.TII  '.'"■.'•.Em  n.INlTlOl!.  PC  Ilf  i^'  nr  Ff.i.UPPEUCE 
PZt:!  :  '.  Ti:i:'t:;  !  .  ;.E  ’  ■  1 1',' :(,ri  ~p:x  IN 

lEFT  '  Ii'i'.:ii:'' t.-i-l- Fh  'MZTNI'I -ICETN 
PEC  ru  -  (■  _  V"'  1 

IZLIN  -  PLT'lMt 
PC  NIC  1  '■  I'.'Ci 

I 1 1 '  ■  1  ■  ! .  r 

py III  =  I  '.I III  ' 
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o  o  o 


PS  IMP  1  P2THP1 
MB  =  MBi-i: 

IF  '.MB.LT.li)  i:0  TO  500 
MB  =  1 

I0P2  =  I0P2+2 

IF  <  I0P2 . LT , ISECT  )  GO  TO  490 
lOPt  =  10P1+) 

I0P2  *  0 

490  CmLL  EXEC';  1  ,  ICNI.IO,  BUFF,  ILEN,  10P1  ,  10P2  ) 
500  COM r I HUE 


QEXT  -  2  .  i'0EXTXeETM'*+2 

IF';MOrtMGS  )  PETIJRH 

RECOVER  SI  MtU'  '52  FROM  S-*-,  S- 

00  300  .1  =  1  ,NH 

SI';  J  )  =•  0,5'fCMPLX'  SUGL'-.  RSP-;  J>  +  PSI1<  :',SMGL';  ISP',  .1 '>+l'rt-P.  0  )  >  > 

S2<  J  I  »  0 , 5>*CMPLX'.  SHIjL'  RSP'.  ..I  >  -PSM<  vl  :>  > ,  SMGL'  I  SP'  J  >-  I SM'',  J  >  .) !) 

3 1 ';  HMP I  -  0  i»o ,  s^cMiH.:-:''  sngl':  rsps',  j  .>+rsmS',  j  :<  :> ,  stJGL':  ISPS'-;  .,i  >+  I3ms<  j  > 
SOO  S2<NMP)-J  )*  0,5  ♦•rMPLX'vSMGU':  PSP'S'  0)-RSM:.'  .J  )  ' ,  •..MGL<  ISPS':  .J  :>- I 'SM'iX  0  :> 

c 

RETURM 

ENO 

SUBROUTINE  OCOVOI.PI  ,R2,S,Ci:> 

C  PERFORMS  COMPiE;;  NIILTIPLICRTIOM  ip  DIVISION  IN  DOUBLE  PREC 

DOUBLE  PREC  1  3 1  ON  R I  ,  R2  ,  'S':  2  > ,  O'  2  ,  D I V 

IF  <R2.EO,0.  .AND.  'S':  2  !> .  EO  0  .  >  GO  TO  10 

DIV  =  S';  1  «'2+'S<  2  > *  +  2 

IF  <  N  1  1,2 

1  0";  1  .>  ~  DSOPTi  '  PI  f  *  2  ♦■P2  I*  <"2  >.  D I ' 

RETUP'N 

2  O'  1  .>  “  ';Pi  >:.’  1  ' •'p;.' »i.',  2  '  >. 'D I V 

I"'''  2  '  "  *;  p«.' ♦  i-’  1  '  ■  P 1  *'S'  2  ' ',  (' ;  V 

RETUPII 

10  O';  1  )  =  R 1  /S',  1  ,) 

Q<  2  I  =  0  . 

RETURN 

END 
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T  =  00n04  1C  OM  CP00023  UUNil  00019  BLK9  P  =  0On0 


C 


pPQur'BM  'IB  IT  ,.e•^^ 

C'IMEM'-li'"l  BP.UF.r  "  >,  !C-Ce<  144  .,HmME*.  - 
PEmL  M-:.  niTEK 

OniJBUE  PF^EClilOH  H.'.r.TH'. 

^  EoS  r-  I  ■  I 

>  '  PBUF'"  4  > .  I  '  '  Bt'Uf  V  5  ■ .  C.2 1  > 

CiiUTP  I 'I 'I  '  iOPTtl.  M.' .  lE.'l 
LU  [N  =  L0uLU‘  IC'iJMIT  •' 

UipiTE.  LUIM„Em,mm-,  .  _  _  _  HtiCLE"  > 

FiIRMm!  ’'EHTEr-  B  1  ■  >  h  ,  F  i  HT  .Lt4 

PEFD'  LUltl,  *•  ■  FI  .BETH.FTS.'r.TRT,  -TwHu 

WRITE'- LU  IN.  con  i  ■  „  ., 

:.,.R 

FORHWT'.  7h:;  ,  IC  ' 

PI  a  T  ,  1  4  I  '5  lO F 5 -3 '0 E’ 0 
p  4  a  F  I  (•'  4  , 

WQ2  =>  '-.QRTi,  C.  ' 

H  »  1  ,C''’.'.'R(  ,  - 

THC  =  OaTBtt':  t  .C'''.'C'30F:TOI*+2-1  Du.' > 

TH2  ~  TMT, 
mi;  =  THC-O  .  WE'C 

th4  »  thc-i  .■?e-7 

PHC  =*  P  I  -  2  .  ♦•TUC. 

H'J  *  M'tl 

=  BErR/C  2  DO' FI  » 
hF;  s  luPT'  H  "'OE''  ■RHi'jL''  THC  ■  I  ’ 

IW  a  'vlPYtTHC 
E  •=  •  TMC-TH  i,  PT  ;. 


2000 

2  0  0  1 
2  OOi 


E4 


F./ 


‘""ofeh^C'mTv.  file  dcsimgated  b ■ 

rpLL  OPEfI'..  IDlO-  Ibl-F.HHllE,  I'-IH  TH.  - - 

IF'-  lEFP -LT  . )  WO  Til  OOO 

PB'.'F'.  I  '  =  BFTh 
ppi.'F'-'-  3  •'  =  '  ■  •’0 

‘^^'.TOTE  FlJMTTiiV'  (-r-Mt'.'TFL)  I'l  FEl'F  ETOiV, 
Ci4L  L  ''IF  ITT"  I  i  'C  B  I EV-  P  ,  B  P'.'l  ,  I  L  ■ 

1  0  '-;T  ^  l-OI"'  C'b-LE-  TW  •  ’  . 

■  ■  ,  T  H  (I  ■  r  ■•■■'■  ’'■T-  \  Pm-'  H  +  '-T  ;  ♦  ♦C  ' 

R  a  Dm  Ti-iH'  is  I  -  .  . 

0  =  P-TH 
?  ,  H' 

_  i;ijC'  TH  ■ 

pe,  r.  PH  iH  Df::.B'E£E 
r-i;*ET  ;  ■■■'S'? 

^  o:,wi'uTE''vHH'e.E  HHLST  DlVpWDttLE  Fhl  IuP 

-  T  f.  1 ;  HH  '  • 

=  ptTH 

TmM'  D  n-'TM'F  ;  -  - 

■  "-UF  r.  ■  ■;  ■  __  , 


Fll 

lT 

F 

pri 


II 


FP 


u  PHCE  'DTwTrMrHT 
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i 


i 


I 

i 


I 


IP  =  G  1  *4  .  ♦ '  t  .  -r-P  ;  t 
I  :■  =  Gt  k4  .  I  .  -FS  > 

ETA  =  PHC-PH 
W  =  S  I  Ai;  ETh  y*AK 
THe  THC>E7A,-'2 

DELI  =  ATmH>,  SC'PT'  SIH'.  THe>*t-2-1  . /NS  >/C  U':':  TAB  .V) 

DEL2  «  2,-*hT^^N<NS>TAH-;DEL1  >  ) 

DELI  =  2.*DEL1 

ETAC  =  2  ,  >£  ETh  COi^  THD  >-C  T  +CO'£';  R  >/M  > 

Ei1=EThC*-P4«-C  EL  1 
B2“ETmi:  +  P4  +  C'FL2 
3I1=G.0RT'  I'  • 

GI2*£.0RT<  IP  :> 

CBI=i:nsi:  B1  - 
CB2=C0£.'.,  B2  • 

'5  B  I  -  S  I N B  I  ,) 

SB 2=  SIM’:  B2 
GOTO  40 

THIS  CODE  USED  MriEM  PH;.-pHi; 

30  ETAD  =  ETAD-CETh 
ETA  =  ETAD, '5?  ,  2053 
PH  =  PHC-ETh 
PD  =  PH’+^F.OOSS 
W  *  SIMC  ETA.)»AF 
S  I  I  *  0  , 

SI2=0, 

cei  =  i . 

CB2= 1  . 

SB1 sO , 

SB2==h  , 

4  0  Call  freS'.  oi, fs,  fc  > 

FSapS/'Sljt 

FC-FC,-'S02 

s  I  P=Fc  tcp  I  +r-'s  ♦sp  1 1  n  I 

S  2  P  -  F  C.  +  C  B  2  -t  P  S .« I  2 
SI  i  =fsh:bi  -(*■:  fSB i 
S2I=FS  m:  B2-Fi:  cSE  :■ 

COMPUTE  FHw-,e  DIPFEPENCE 
PH)  =57 , 2A50+H  I  hM’.  S!  1/'S1R  '/ 

PH2  =  57 . 2050  *  hTmMi:  :S2I  .•‘S2P  ) 

I F'.',  S  I  P  ,  L  T  ,  ij  ,  F  H  I  =  t ij  ,  +F  F|  1 
IF'  FH  t  ,  G  7  I  ■?  C  .  '  F  H  I  =PH  I  —  'It-  0  . 

I F','  SPP  ,  L  T  ,  0  ,  >  PHr;=  t  3  0  .  +PH2 
IF':FH2.GT.  I  C’O.  F'H2=PH2-7'7i'i , 

DIFh=PH2-PM I 

IF':  DIFh  ,  GT  ,  I  00  .  DIFm-PIFa-SE  0  . 

IF'^C'IFk.LT,  -  f  ij  ,  t  L  I  (-  H  =  C'  I F  A  <-7';.  ('  , 

COMPUTE  total  IHTEli-.ITIEC 
T  I  S  = '  .  S  I  P  +  t'  2  +  S  1  I  *  2  > 

TIP  =  '  S2F'-*2  +  S;i  •  *2  : 

IMTEM  =  ALOGr’  Sf,  T  ISaTIP  .>  > 

IF'' MroOD':  PD,  1  E -4  < ,  EQ  ,  n  ,  .  i.'P  I  TE',' LUI M  ,  1  0  0  n  ;)  F  D  i' I S  ,  T I P  ,  D  I  Fa 
0  Fi'.iRMhT'  '=’8 . 4  2E  I  4 . 0  .  F  I  0 , 4  ) 

CALL  I  IF  I  TF':  I  DC  B  ,  I  EPP  ,  PBUT  ,  I L  > 

IFOH.LT  THC  )  GOTO  SO 
IF':  PD  ,  LT  ,  S  I  aMG  >  C'DTO  S'l 
GOTO  O'V? 

5  0  IF';  TH  .  GE  ,  Tri:  ■  E^E2 

]  F'  10  ,  GE  ,  rU  T  •  L-  E  7 
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ETrtOeO, 
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GOTO  30 
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RETURN 
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END 
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Paper  No.  7 


Glory  in  tlie  optical  backscatteriru)  from  air  bubbles  (U.  S.  Lanyley 
and  P.  L.  Marston)  Accepted  for  Publication  in  Riysical  Review 
Letters 


Abstract 

Obfiorvat -tons  of  1  tglit  backaca l  toroJ  from  air  bubble's  In  viscous 
liquid  demonstrate'  an  onhauccmenit  due  to  axial  fociuiio,;.  A  physical-optics 
approximation  fo’-  the  cross-polar;*  xe'd  scatLcrinq  correctly  dc'seribe'S  the 
spacing  of  regular  features  observed,  Tfie  non-cruss-1'’f'larize’'d  scaltering  is 
not  adequately  elescrlhed  by  a  single'  e-lass  e>f  I'ays. 
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The  Mle  solution  for  electromagnetic  scattering  by  a  spliere  frequently 
does  not  lead  to  direct  interpretation  of  the  angu.’ar  scattering  pattern.  Con¬ 
sequently,  models  have  been  developed  to  facilitate  an  understanding  of  the 
structure  in  the  scattered  intensity  present  where  intensity  is  plotted  as  a 
function  of  the  scattering  angle  (j)  or  the  size  parameter  x  =  ket  (k  =  wave- 

number;  a  “  sphere  radius).  These  models  have  emphasized  tlie  angular  regions 

2  3 

where  diffraction  is  important  lor  a  drop  of  water  ii\  air:  the  rainbow,  ’ 

(fi  =  180“,^  ^  and  (p  ^  0°.^’^  In  the  scattering  of  light  by  a  spherical  air 

bubble  in  a  liquid  or  in  glass,  the  real  part  of  the  refractive  index  of  the 

sphere  is  less  than  that  of  the  surroundings  and  the  models  must  be  signifi- 

7  8 

cantly  modified.  New  phenomena  appear,  sucli  as  diffraction  ’  in  the 

region  of  the  critical  scattering  angle  '|)^.  Here  we  report  the  first  detailed 

observations  of  backscattering  by  air  bubbles  in  liquids  and  give  a  model 

which  describes  some  of  the  observed  features.  Wc?  refer  to  this  as  glory 

3-5 

because,  as  in  tlie  case  of  drops,  the  (fi  180"  scattering  is  enhanced  when 
X  is  large. 

3  4 

Van  de  Hulst  ’  gave  a  partial  explanation  of  the  enhancement  for  crops 
by  noting  the  axial  focussing  of  those  backsca ttered  rays  which  have  a  non¬ 
zero  impact  parameter.  Wlicn  modeling  this  focussing  in  the  far  field,  dif¬ 
fraction  provides  an  e.ssential  correction  to  ray  optics  because  the  factor  in 
the  scattered  intensity  v/hich  accounts  for  geometrical  divergence  of  the  rays 

goes  to  «>  as  4^  ’■  180".  Examination  of  this  factor  in  ray-optics  models  of 

g 

scattering  by  bubbles  sliows  that  this  Is  not  restricted  to  drops.  We  have 
modeled  the  backscattering  with  a  physical-optics  approximation.  The  proce¬ 
dure  is  to  (a)  compute  amplitudes  in  an  exit  plane  in  contact  with  the  bubble 
via  ray  optics,  and  (b)  allow  i.his  wave  to  diffract  to  the  far  field  where  the 
distance  from  the  bubble's  center  R  >>  kfl*'. 
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Fig.  1  lllustratc.s  several  rays  which  lead  to  hackscattering.  The  patlis 

are  determined  by  the  number  of  chords  p  and  ni  ~  m./m  whi-re  the  refractive 

1  o 

indices  of  tlie  inner  and  outer  media,  m,  and  ni  ,  are  taken  to  he  real.  Fig.  1 

1  ( 1 

is  drawn  with  m  ^  =  1. A03  which  corresponds  to  an  air  bubble  in  the  dimethyl- 
siloxanc-polymer  liquid  used  in  the  experiment.  All  rays  .satisfy  .sln0  =  nisliip. 
For  ({)  “  180°,  tb.e  off-axis  (or  £l_or^)  rays  h.ive  o  =  t)  and  f)  =  p  where  0  = 
pp  +  (2g+2-p)90",  g  i.s  a  non-ni-)',at  Ive  integer  (g  =  0  for  rays  in  Fig.  1)  and 
m  <  1  requires  tliat  p  >  3.  The  exit  plane  ((la".lu'd  line  in  IMg.  1)  touches  C' 
with  its  normal  parallel  to  the  propa^’.ai  Imi  dlrei  tion  o1  the  incident  wave. 

Our  description  ol  the  fii-ld  in  (he  t'x  i  t  plane  is  lacilitated  by  consid¬ 
ering  the  propagation  ol  a  wavelet  de  which  lies  i  Inso  to  the  backsc.attered 
path.  Fig.  1  show.s  de  for  p  “  .1;  if  enu'ig.e.s  as  curve  d'e'.  J'hi.s  curve  appear.s 
to  come  from  a  ring, -like  sottree  at  F  known  as  the  loc;.l  circle  in  the  analo¬ 
gous  p  =  2  scattering  from  droju-;*  with  /?  •  m  •  Z.  The  source  is  ring-like 
because  tlie  figure  may  be  rotated  around  the  ('.f,'  axis.  The  radiu.s  of  the  ring 
is  b  ■=  dslnO.  After  the  iueldent  lay  cross.es  the  dashed  vertical  plane  (the 
entrance  plane),  the  ptej'.ijjatJ^tni  phase  delay  for  if'acli  i  ng,  the  exit  plane  is 
n  ••  ka[l  -  cost)  +  (1  -  cosB) sec  (O-B)  +  2mpeosiil.  The  ray  crosses  the  exit 
plane  at  a  radius  s  from  C’  with  s/a  =■  sint'  -  (1  -  cosB)  tan  (')-B)  .  The 
radivis  a  of  arc  tl'e'  follows  from  the  curv,(i\it('  at  s  •  b:  a.  =  k(d  t|/ds^)  - 

a[l  +  ''.i(pT-l)  ^cosb)  where  i  -  tani'/tanO.  The  spreading  of  the  wavelet  is 


characterized  by  q  -  11m  d'eT/de  as  do  ►  0  where  the  bar  denotes  the  arc  leng.th. 

An  equivalent  expression  for  q  is  |llm[h  -  s(0)l/(b  -  a.sln(J)|  as  0  a-  0; 

its  value  from  1, ' Hospi lai ' s  rule  is  a/ (a  -  a).  Vectots  e^  {I  =  1,2)  denote 

orthogonal  basis  vectors  In  both  the  enti.ince  and  exit  planes;  is  chosen 

parallel  to  the  pol  .arlzat  1  on  of  the  incident  w.ive  '  s  ei((ctrtr  field  E^^exp  (-iuit)  . 

0 

Tn  the  exit  pl.ine,  the  fiel<l  of  t  ue  outgoing,  pth  glory  wave,  is 

computed  by  applying  Van  de  Hulst’s  method  of  first  decomposing  the  fields 


4 


perpendicular  and  parallel  to  the  scattering  plane.  ’  Exit-plane  polar 
coordinates  centered  on  C'  are  where  is  the  angle  relative  to 

and  8  and  i(j  denote  local  basis  vectors.  We  assume  x  >>  1  and  use 
Fresnel's  coefficients  r^  for  the  internal  reflections  where  j  ®  1,2  for 
fields  parallel  to  ij'  and  s,  respectively.  If  |s-b|  <<  CL,  the  multiple 
internal  reflections  give: 

Ep  -  F^exp[in  +  ik(s  -  b)^/2al  (1) 

where  n  “  M  +  n(p  =  P),  F  (>10  =  c^sin'4/  +  c^cos  'p,  F  (i|;)  =  -  c^)sin2t|^, 

and  Cj  ••  (-1)^^^  phase  term  p  accounts  for  the 

crossing  of  caustics  or  "focal  lines";  its  value  is^’^^  -■'f(p+8)/2.  The  r^ 
are  evaluated  at  9:  r^  =  sin(0  -  p)/sin(9  +  p) ,  “  tan(0  -  p)/tan(9  +  p) . 

The  sign  factor  in  Cj  acconnt.s  for  a  geometrical  inversion  (present  when 

j  =  2  and  p  Is  odd)  which  is  not  evident;  in  descriptions  of  p  =  2  glory 

.  j  3,4 

in  drops. 

The  field  at  a  distant  point  Q  is  computed  as  follows.  The  left  exten¬ 
sion  of  the  CC  axis  makes  an  .ingle  y  with  C't).  When  y  is  small  and  C'Q  = 

2  12 
R'>>ka  ,  scalar  diffraction  tlieory  and  the  Fraunliofer  .approximation  give: 


kBj  ^ 

27tiR'  q  ‘ 


ik(s-b)  /2a^ 
iW  e  ds 


^-iks  sinYcos(t|J-f;)^^ 


where  ^  is  the  angle  between  and  the  projection  of  C'Q  on  tlie  exit  pl.ine. 
In  Eq.  (2),  the  approximation  given  by  Eq.  (1)  has  been  extended  beyond  its 
useful  domain  in  anticipation  of  the  stationary  phase  approximation  (Si’A)  of 
the  integral.  Direct  evaluation  of  Kq.  (3)  gives  W  (y,0  W  (Y,C»s  =  b) 
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“  +  C2)Jq(u)  +  (c^  -  C2)J,,(>i)  cos2f  J  and  W*"  -  ii  (c^  -  c^)J2(u)  .‘•^u2C  where 

u  =  kbslny*  The  SPA  ol  Eci.  (2)  Rives  tlie  prh  Rlury  contribution  to  the  scat- 

O 

tered  field  when  kb  /u,  and  thus  x,  are  larRe.  In  tlie  experiments  to  be 
described  x  ,>  AOOO  and  the  SPA  Is  applicable. 

The  total  field  may  be  approximated  bv  summlnR  the  from  Kq .  (2) 
with  the  fields  due  to  axial  reflections  aird  fjurl  ice  waves.  Surface  wave 
contributions  should  bo  small  for  tire  observed  bubbles  diu;  to  the  largeness 
of  X.  To  determine  which  glory  and  axial  terms  .art'  important  to  tlie  total 
field,  and  for  othei  heuristic  reasons,  is)irsid<‘r  the  x,-polarlzed  intensity 

y, 

1  of  the  pth  field  ta_keji  alone.  The  SPA  ol  Kq .  (2)  j’.lves: 

P 

-  (2./n)xl  f  IW^  ,  (A) 

P  k  P . 

2  2 

where  1  =  la  /AR  is  the  total  intensitv  at.  a  d  i  si  aucc'  K  =  Ctj  from  a 

K  i 

9 

perfectly  reflecting,  splnu'e  of  r.tdius  a  predicted  by  lay  optics,'  1^  is 

2  1  2  i 

the  incident  intensity,  and  f  b'  a/a  n  ^  b  (a  -  a)/a  .  in  Eq.  (A),  R 

P » b 

has  replaced  R'  from  (2)  and  t  becomes  ISv)"  -  [j  because  R  •  <t. 

3,9  ■'? 

Geometrical  optics  I’.Lvi’f,  the  intensities  I  oi  si'parato  axial  reflections 

P 

(e.g.  p  =  0  and  2  in  IMg,.  1)  wb  i  h  aic  propt)i  t  i I'lia  1  to  a".  The-  strong, e.st 

reflectifin  has  p  ~  0  and  .t'.  •=  1;  for  >  -•  0,  l'  ^  1  (m- 1  I '  /(mil)'  wliile 

V/  K 

"2  ii.  ( 

1„  =  0.  Since  f  does  iu)t  depend  on  a,  1  ■'  ka  .uid  glorv  terms  dominati;! 

0  p ,  g  p 

the  'nackscal  t  er  1  ng  when  a  is  lai  g.e. 

Consldei  a  buld)le  with  x  AOOO  and  m  1 .  At)  !  The  strongest  glory 

terms  have  g  -  0  and  ]'  3, A,  and  S;  the  ’*/’(>  1  '  respectively 

1.03,  0.A3,  and  0.16.  The  1 decrease  with  Increasing  p  due  to  the  partial 

reflections  in  the  bubble.  The  strongest  axial  r.iy  gives  “  0.028.  The 

(J  R 

interference  ol  the  flelrls  depends  on  a  anil  our  Mli'  computations  verify  that 


the  backscattered  Intensity  Is  not  simply  propoi t Iona  I  to  a  even  for  this 
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large  value  of  x.  The  S,  2  (cross-polarized)  scattering  is,  however, 

nearly  dominated  by  the  p  =  3  glory  term.  Due  to  symmetry,  1-2  scattering 

2 

vanishes  as  Y  ^  •  The  I  (y  ^  0,f;)  have  maxima  at  =  *^45'’  and  4135“  and 

P 

they  vanish  at  f,  =’  0“,  +90“,  and  180°.  Let  y  =  j  locate  the  first 

P 

2  2 
maxima  of  “  45°).  The  largest  I  ~  2  terms  have  “ 

0.53  and  0.10  for  p  =  3  and  4.  To  the  extent  that  ^  i  scattering  may  be 

neglected,  the  2,  =  2  Intensity  will  be  quasl-per indie  in  y. 

We  have  numerically  verified  the  validity  of  Eq .  (4)  by  using  Debye's 
3  4 

localization  principle  ’  to  modify  Mle  theory  so  that  only  partial  waves 

associated  with  p  =  3  rays  were  included  in  the  Mie  series.  Furthermore,  when 

Eq.  (4)  is  applied  to  spheres  with  certain  m  >  1,  the  resulting  1^(Y  =  0) 

P 

agree  with  the  glory  "analog"  tabulated  in  Ref.  11.  This  analog  was  derived 

by  applying  the  Watson  transformation  to  the  Y  “  0  Mle  series. 

Fig.  2  diagrams  the  experiment.  A  syringe  injected  bubbles  into  the 

liquid.  The  liquid  had  a  higli  kinematic  viscosity  (=600  000  cS)  and  a  single 

bubble  could  be  observed  for  hours.  The  laser '.s  power  output  was  5  mW  and  the 

beam  diam  was  5  mm.  The  wavelength  in  the  liquid  2Ti/k  was  632.8  nm/l.403; 

e^^  lay  in  the  splitter's  plane  of  incidence.  The  camera  was  iocuse.d  on  ™  so 

7  ]  2 

the  photographs  recorded  the  far-fleld  intensity  pattern.  ’  Photographs 
were  made  with  a  =  0. 3-0.8  mm  corresponding  to  x  -  4000-11000.  Exposure  times 
were  typically  5s  for  TrlX  film  and  a  200  min  focal  length  camera  lens. 

Fig.  3  demonstrates  that  the  scattering  has  roughly  the  dependence  on 
C  predicted  by  Eq.  (4);  f,  =  0°  corresponds  to  scattering  toward  the  top  of 
the  photographs  and  y  =  0°  corresponds  to  the  center  of  the  symmetry.  Fig.  3(b) 
shows  that  the  ?  =  1  scattering  for  Y  ^  0.2°  is  significantly  stronger  for 
^  =  ±90°  than  it  is  foi  f,  =  0  .  This  agrees  with  the  following  model  results: 
(1)  >>  1  (for  p  =  3  we  predict  c^^/c^  -  -5.2);  and  (ill  for  this  x. 
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1 

the  depend  only  weakly  on  £,  and  are  dominated  by  the  .  One  prediction 

of  Eq.  (4)  could  be  quantitatively  checked;  when  both  siny  -  Y  and  u  >>  1, 

2 

the  minima  in  should  be  spaced  by  Ay  rad  such  that  kbAy  =  tt  where  for 

p  =■  3,  b/a  =  0.447.  Fig.  4  compares  this  with  the  mean  spacing  of  ~40  dark 

rings  lying  outside  the  9th  ring  from  the  center.  The  error  bars  combine 

uncertainties  in  measured  a  and  Ay  with  those  of  corrections  due  to  refraction 

at  the  cell-air  interface^  and  the  tilt  of  the  cell.  Fig.  4  shows  that  p  =  3 

rays  dominate  the  1=2  scattering.  Tlic  modulations  of  the  intensity  along 

^  =  ±45“  in  Fig.  3(b)  show  that  other  rays  contribute  to  X,  =  1  scattering 

1  2 

since  the  predicted  «  [dQ(o))  . 

In  conclusion,  backseat tering  from  bubbles  can  be  enhanced  by  axial 

focusing.  The  number  of  significant  glory  terms  depends  on  m.  The  main 

3 

contributions  differ  from  those  for  water  drops  where  surface  waves  and 

other  diffraction  related  tt'rras^  play  an  essential  role.  If  focusing  were 

7-9 

not  present,  scattering  by  large  bubbles  would  be  '^<  1  in  the  region 

R 

(<()^  +  10°)  <  <{)  £  180°  where  =  2  cos  ^m  =  89°  for  m  ^  =  1.403.  We  also 
find  evidence  of  p  =  3  glory  In  Mie  computations  for  bubbles  in  water . 

This  work  was  supported  by  the  Office  of  Naval  Research.  P.  L.  Mnrston 
is  an  Alfred  P.  Sloan  Research  Fellow. 
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Fig.  Captions 

Fig.  1.  Rays  which  contribute  to  hackscattering.  Tlie  local  angle  of  inci¬ 
dence  is  0  and  C  is  the  bubble's  center. 

Fig.  2.  Apparatus  for  observing  backscattering  from  bubbles. 

Fig.  3.  Photographs  for;  (a)  crossed  polarizer  (£.  =  2  scattering); 

(b)  uncrossed  polarize.  (i.  =  1);  and  (c)  no  polarizer.  The  incident  polar¬ 
ization  was  vertical.  0.  =  0.49  nun  and  x  =  6830. 

Fig.  4.  Measurement  and  model  for  the  angul..r  separation  of  the  dark  rings 
in  the  S,  =  2  scattering. 


84 


ig-  2 


86 


Fig. 


I 


i 


[ 

I 


Vi 


CM 

in 


CM 

C7N 


L.  r-  O 
0  3  • 

r—  0.-«  »T3 

>)  O  05  c 

ro  u  o  a> 

i—  Ck.  TJ 

o  m 

•  W  O  i'* 
1.  O)  CO  0) 
£3  O  ot 


K  C 

o  v_  c  o 

O  >  C7>  w 

V.  O  O  ** 

•  O  C  05 

00  +J  fO-r-*— 

o5aj^f— 

•  u  u  m  o 

o  o  o 

.o 

V.  m  •  *4-  05 

o  _J  +J  o 

in 


>y  >^  C  > 

d;  ♦J  I— t  4 


4->  .—  0) 

C  *.  •»- 

3  ^  .  dJ  C3 

Q  > 

lA  u  c 
•  u  c  05 
cr^  i/>  =5  ui 


o  >> 

O  -c 


C  <*.  ‘r-  (tj  ^ 


c  w  ^ 

x:  c  05  m 

o  05  'r-  - 

>r>  i-[  V*  X 

•P  3  m  “t 

•  0.  o  ^  o 

Q  0»  X:  -r-  <t 

Q.»~  r- 

•  dJ  03  01  <C 

CD  O  Q  X  O 


V.  O  l/>  cn 


w  E  r- 
O  (U  CM 
•!-»  05 


t/5  rtj  O 

.-I  r-  tC 
t4j  >»r5 


•  r>  O  (/>  dl 

^  U  r- 

<c  •  45  w 

*J  >  4^ 

•  a /a 

U  45  C  45 
o  o  r)  i/i 


O  45  o 

O  U  05 
3:  O  45 


45  Vi  • 
E  ^  c^ 
■M  45  • 

c  >  o 

I  03  •»“ 

Q.  tr  • 

)  45  ro  c: 

I  o  o 


X  O 
io  x:  x: 
.  to 

U  -C  05  05 
O  O  3 


CL  ^ 

O)  CM 
U  O 
CT*  • 

O  ^ 

c 

03  •CO 
45  •—  CT\ 
1>  O  •—  O 
45  O  05  Cs5 
t—  X  o> 
>>  • 

5—  4-»  CL 
VI  3  «t 
•  C  I-  L_> 

UJ  *-•  -O 


C  •  45  03 

x:  v)  >  f“ 
o  4J  in  ^ 

'-3  Q.  O 
•»-  CM  *3 
•  k.  ^ 
k.  u  CM  to 
O  C/5  CM  -J 


>. 

(■'i  U 
r—  o  1/5 
U5  M  C^ 
03  CO 
45  i.  O 
T3  O  CM 
O  X) 

CJ  05 

—5  * 

•  1.) 


u 

45 


O 

5: 


3  45  O  Cn 
05  (J  C 
CO  05  45  ••• 


.  U  -  CL 
(_>  J  X  lO 

to  3 

L. 

Or-  *45 
05  •—  > 

•  >  m  »— 
u  lO  t  •»“ 

O  sc  QC  to 


•r-  4)  O  05 
2  O  •  45 
(5  CD  O.  to 


5'S 


>»ro 

<-  O 

t51  O 

U  r- 

45 


.  I—  > 

•3  ca.**-  >• 

CL  c  z 
C  ct  3 
QJ  • 

i.  M 

i>  o  t.  u 
0)  o  o 


L.‘  ^  {  S 

Q  O  Z  Z 


41 

O 

fs. 

«o  O 

-- 

4-»  ro 
45  iO 


a> 


-to  •* 
X  *-*  >» 
03  1/1  LJ 
r-  to  'i~ 
U.  O  tJ 
tj 


o 

•r* 

•4-5 

3 


00 
<1.  :7i 
45 
1/1 


•4-5  O  =3  -M 

O  k- 


41 
C  -o 
-  o 


s 


Q. 

>> 

41 


<>! 


45 

4-5 

01  O 
4-5  3  Of 

c  *o  a> 

45  05  m 

e  w  cr» 
P  OI 
C  5-  4-5 

•r-  03  to  < 

3  CL  O  O 
X  •-  >» 

i/i  <e  45 

•  (J  >  L. 

X  -f-  41 

to  Z  -t-i 

•  ^  C 

k-  £  to  o 
a  D.  =3  X 


n.  >> 

<  V.  : 

o  * 


•.J  O  45  CM  • 


45 

n 

Q 


»  •»“  X  45 
5  C  O  > 
3  00  05 
X 

>  45  • 

1  ^  O  » 
J  03  •  ^ 

:  >•  a.  z 


o 

U  to 


J/  ti  K-  r-. 

u  1-  CO 
05  05 
r—  45  O 


XJ 

•  01  -  — 

CO  •»-  45  ••- 


>  4.5 

a.  io 

CL  C 


r-1  •r^  *#■ 


$ 

C 


03  (/I  ^ 
U  O  O 
^  .r~  O  O 

o  cr  X 

O  to  r-K 

O  -C  »4- 

O  O 

.  45  X 

Q  X  >»♦— 

4-5 

^  *4-  •»-  • 

^  o  V)  C 


rt**-  3 
o,  c  O 

O  3  X 


J  >>  VI  I 

x:  •- 

•  o.  X 


V  45  i/l  *»- 
3  S  I-  lo 
03  P  k.  U 

-J  W  41  05 


05  vi  M 
^  k. 

L.  41 

.  > 

01  4J  •»-  «£ 

U.  Qu  cr  O 

01  :--j 

.  o  • 

X  T?  t3 

VI  k-  u 

•  u  o  o 

ac  .F-  V-  w- 

VI  r  c 

•  >S  03  03 
L.  X  4->  4^ 

CO  a.  VO  wO 


y-  VI  L. 

■6  -r-  o  X 
.t;  >  *-*  z 

05  o 


•j:  03  -J  u 


ct  ^  cr 

(J  *0  3 
•  O  C  43 

W  C  0)  r- 

o  to  VO  «t: 


5  4)  >tm 

f—  x:  4-5  o 

. —  *J  -r-  ro 

45  05  V>  ^ 

Lt  ri;  k_  o^ 
45 

•  1*-  > 

CD  O 

c  U) 
x:  ■*->  :3 
Q.  C  • 

-  .  o  o 

"O  U.  ‘4-  *4- 
03  C  C 
•  to  0) 

c.  4(  4-i 


a.  c 

w 


c  ^ 

U  CM 

o  o 

*4-  O 

4-5  T-  Ol 

JSg 

U  V-  v 
•P*  to  o  • 
C  CL  VI 
03  05  >>  05 
3  O  4-5  p— 

oc  -tv 

VI  V>  Cfi 
•  U)  t  5 
t-x  •*-  05  4 
VI  > 

.  >>*r-  VI 

c  z  c  o 

Q  t3L  3  -J 


i 


1 


i 

i 


s 


